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SECTION  1 


INTRODUCTION,  OBJECTIVE,  AND  SUMMARY 


The  role  of  the  Charles  Stark  Draper  laboratory,  Inc.  (CSDL)  in  this 
study  has  been  to  assist  the  government  (USAE/RADC)  in  assessing  the  via¬ 
bility  of  a  space  based  millimeter  wave  (MMW)  antenna.  While  several  other 
contractors  were  involved  with  the  study  of  the  various  aspects  of  the  de¬ 
sign  of  a  MMW  antenna  system,  CSDL's  responsibility  was  to  provide  tech¬ 
nical  support  in  the  areas  of  structures,  dynamics,  and  control  for  a 
specific  antenna  concept.  The  results  of  these  studies  provide  information 
for  use  in  comparing  a  MMW  surveillance  satellite  with  other  surveillance 
satellites  using  visual  or  infrared  wavelengths. 

The  objective  of  the  CSDL  work  is  therefore  to  design,  analyze,  and 
evaluate  the  performance  of  a  MMW  antenna  system  concept.  The  base  line  was 
selected  from  the  Aerojet/NTL  study.  The  design  should  rely  on  a  traditional 
approach  using  state-of-the-art  technology  for  materials  and  spacecraft  de¬ 
sign  techniques.  More  specifically,  this  -tudy  is  limited  in  scope  to  the 
generic  truss  constructions  that  are  generally  proposed  for  large  satellites 
requiring  high  precision  for  their  optical  components.  The  truss  structure 
used  in  the  design  of  the  MMW  antenna  can  be  compared  to  an  alternate  struc¬ 
tural  concept  employing  an  electrostatically  shaped  reflector.  Figures  of 
merit  for  both  systems  will  be  compared  in  order  to  assess  their  relative 
desirability. 

In  designing  the  structure  of  the  satellite,  all  avenues  are  explored 
in  order  to  obtain  an  adequate  performance  while  using  passive  means  of 
control  only.  An  active  system  could  be  advocated  lor  instance,  on  the  basis 
of  comparison  between  the  a -vantages  of  weight  saving  and  problems  associated 
with  the  added  complexity.  The  decision  between  active  and  passive  approaches 
depends  mainly  on  reflector  size  and  its  associated  deformational  tolerance  re 
quirements.  Among  the  passive  means  of  controlling  the  antenna  deformations, 
the  judicious  design  of  thermal  blankets  appears  to  be  the  most  promising  ap¬ 
proach.  This  i3  a  direct  consequence  of  the  fact  that  the  thermal  effect  of 
the  solar  radiation  constitutes  the  more  severe  disturbance. 

In  tae  present  study,  the  numerous  issues  involving  packaging  and  deploy¬ 
ment  are  also  considered  as  they  impact  the  structural  design.  For  a  reflecto 

■) 

of  the  size  considered  (-1600  m”)  a  relatively  high  packaging  ratio  is  also 
needed  if  the  satellite  is  to  be  stowed  into  the  cargo  of  a  single  shuttle. 
This  very  desirable  requirement  restricts  the  freedom  of  the  structural  de¬ 
signer. 

In  earlier  studies,  done  by  other  contractors,  the  most  promising  candi¬ 
date  sensor  configurations  have  been  identified.  The  present  satellite 
design  is  based  on  the  single  surface  antenna  concept  described  in  the  report 
by  Aerojet  Electrosystems  Company,  (AESC)  and  NTL,  Inc.*  The  elliptical  torus 
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reflector  is  rectangular  in  projection,  with  side  dimensions  of  approximately 
60  m  and  30  m.  The  feed  surface,  about  30  m  above  the  reflector  is  supported 
by  a  tower  designed  in  a  way  that  its  members  do  not  interfere  with  the  in¬ 
coming  rays. 

This  report  presents  the  approach  to  the  design  and  the  analysis  of  a 
structural  system  capable  of  supporting  the  millimeter  wave  antenna  concept 
suggested  by  AESC/NTL.  The  structural  distortions  resulting  from  environ¬ 
mental  disturbances  and  on-board  dynamic  loadings  have  been  translated  into 
antenna  performance  indices.  Line-of-sight  errors  and  loss  of  gain  have  been 
calculated  for  various  beam-pointing  directions  which  use  different  portions 
of  the  antenna  surface.  It  appears  that  it  is  possible  for  the  performance  to 
fall  within  acceptable  limits,  provided  that  the  surface  panels  are  fabricated 
with  graphite-epoxy  composites  because  of  their  high  modulus  of  elasticity 
and  iow  coefficient  of  thermal  expansion.  The  design  of  the  thermal  blanket 
insured  small  temperature  excursions  and  negligible  gradient  such  that  the 
surface  of  the  panels  was  oniy  slightly  distorted.  The  moti  *s  of  the  re¬ 
ceiver  supported  at  the  tip  of  the  tower  were  more  signifies  it.  Some  active 
positioning  scheme  as  a  function  of  time  along  the  orbit  path  might  be  re¬ 
quired,  unless  some  judicious  compensation  technique  is  designed  to  provide 
thermal  insensitivity  for  the  30  m  open  tower.  Because  of  the  500  km  orbit 
that  was  tentatively  selected,  the  atmospheric  drag  causes  a  substantial 
orbital  decay  rate,  thus  creating  a  need  for  additional  propellant.  Preliminary 
estimates  put  the  needed  propellant  at  2000  kg.  The  structural  design  concept, 
including  124  surface  panels,  the  back-up  truss  assembly  and  the  receiver 
tower  can  be  packaged  into  a  single  shuttle  cargo  bay.  While  the  packaging 
scheme  presented  is  still  preliminary,  the  use  of  open  cross-sections  for  the 
truss  construction  greatly  increases  the  packaging  ratio  without  endangering 
the  stability  of  the  individual  members.  As  is  usually  the  case  for  large 
high  precision  satellites,  the  total  mass  of  8800  kg  makes  this  design  volume 
limited  with  respect  to  shuttle  capabilities. 

This  report  sequentially  covers  the  steps  involved  in  identifying  the 
problem,  designing  the  structure,  constructing  a  mathematical  model,  quanti¬ 
fying  disturbances,  obtaining  the  structural  response  and  evaluating  the 
cor-  esponding  performance.  It  explains  the  approaches  used,  with  their  assump¬ 
tions  and  limitations.  It  is  important  to  recognize  that  the  performance  of 
the  antenna  structure  has  been  evaluated  as  the  satellite  responds  to  typical 
static  and  dynamic  disturbances  expected  along  a  500  km  altitude  circular 
orbit.  The  main  objective  of  this  study  was  to  check  whether  a  satellite 
of  this  size  (60  m  x  30  m)  could  be  designed  to  perform  adequately  at  a  3  mm 
wavelength  in  a  purely  passive  way.  The  report  includes  an  introduction, 
followed  by  six  sections  and  a  conclusion.  An  outline  of  the  content  of 
each  section  is  presented  below. 

Section  1  above  defined  the  problem  and  the  objectives  of  the  study.  The 
approach  chosen  in  order  to  provide  a  useful  answer  to  the  problem  was  ex¬ 
plained  with  the  assumptions  and  the  reasons  involved.  A  summary  of  the 
important  results  is  included. 
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Section  2  deals  with  the  structural  design  for  the  millimeter  wave 
antenna  based  on  a  concept  suggested  by  AESC/NTL.  It  describes  the  approach 
followed  in  selecting  the  truss  configuration  and  the  structural  materials, 
the  stiffness  needs,  and  stability  criteria.  Some  design  optimization  under 
the  constraints  of  packaging  and  deployment  requirements  have  also  been 
attempted.  A  summary  of  the  mechanical  and  inertia  properties  of  the  satel¬ 
lite  components  is  presented. 

Section  3  addresses  the  problem  of  disturbance  identification  and  quanti¬ 
fication.  Steady-state  and  transient  thermal  loadings  are  calculated  for  a 
number  of  thermal  models  in  order  to  minimize  the  expected  responses.  Models 
of  the  reflector  assembly  and  the  receiver  tower  have  been  analyzed  separately. 
An  evaluation  of  the  atmospheric  drag  is  presented  for  two  extreme  solar 
activities.  The  effect  of  solar  radiation  pressures  were  obtained  along  the 
entire  orbit  and  expressed  in  terms  of  force  and  torque  resultants.  Because 
of  its  asymetrical  configuration,  the  satellite  is  subject  to  additional 
torquing  from  the  effects  of  gravity  gradient.  All  these  force  and  moment 
resultants  translated  into  dynamic  disturbances  of  thrusters  needed  to  compen¬ 
sate  for  loss  of  pointing  accuracy  and  for  station  keeping  requirements. 

Other  mechanical  disturbances  such  as  the  harmonic  forcing  function  of  a 
retrigeration  unit  have  been  identified. 

Section  4  discusses  the  approaches  to  structural  modeling.  While  the 
model  for  static  analysis  derives  its  simplicity  from  the  permissible  decou¬ 
pling  of  the  surface  panels  and  the  back-up  trusses,  the  dynamic  analysis  is 
less  rigorous  since  the  decoupling  is  based  on  the  frequency  separation  of 
relatively  rigid  panels  and  flexible  main  structure.  Two  distinct  models, 

(a)  the  back-up  trusses  with  the  receiver  tower,  and  (b)  the  typical  surface 
panel,  have  been  constructed,  and  modal  analysis  furnished  information  on 
their  dynamic  characteristics.  For  the  surface  panels,  several  materials  and 
construction  approaches  have  been  considered  so  as  to  suggest  the  most  effi¬ 
cient  solution. 

Section  5  explains  the  choice  of  static  and  dynamic  disturbances  believed 
to  be  critical  for  the  antenna  performance.  The  thermal  loading  due  to  the 
sun  flux  at  68°  incidence  angle  with  the  pointing  axis  was  treated  in  detail. 
Four  different  zones  on  the  reflector  were  subjected  to  four  different  tem¬ 
perature  excursions  and  gradients.  The  receiver  tower  was  alsc  analysed  for 
several  thermal  conditions.  Dynamic  responses  were  calculated  for  three 
disturbances.  The  first  two  are  derived  from  two  possible  means  of  attitude 
control  through  thrusters  impulses.  The  third  dynamic  loading  was  a  harmonic 
forcing  function  representation  of  a  cryo-cooler  machine  needed  for  the  re- 
*  frigeration  of  the  receivers.  Absolute  distortions  of  reflector  surface 
points  were  obtained  as  a  geometrical  combination  of  individual  panel  distor¬ 
tions  and  back-up  truss  deformations,  for  all  the  loading  conditions  mentioned 
above. 

Section  6  describes  the  degradation  of  point-spread-functions  (PSF)  ol 
the  30  m  x  60  m  elliptic  torus  antenna  due  to  its  deformations  produced  by  so¬ 
lar  flux.  The  antenna  is  assumed  to  be  operating  at  a  wavelength  of  3  mm  with  a 
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cosine  feed  taper,  and  the  materials  considered  are  graphite-epoxy  and 
aluminum.  It  is  shown  that  the  graphite-epoxy  antenna  performs  much  better 
than  the  aluminum  antenna.  Strehl  ratios,  line-of-sight  error  (LOS), 
power  on  a  detector,  and  profiles  of  the  PSFs  are  given  for  three  regions  of 
the  antenna.  Airmat  antennas  were  also  considered,  but  owing  to  their  large 
deformations  (~A)  under  solar  flux,  they  were  not  investigated  in  any  detail. 

Section  7  deals  with  the  issues  of  packaging  and  deployment.  The  main 
objective  was  to  limit  the  required  volume  to  one  shuttle  cargo  and  suggest 
a  deployment  scheme  that  is  as  close  as  possible  to  an  automatic  maneuver.  A 
packaging  method  presented  for  the  back-up  trusses  used  open  cross-section 
members  and  greatly  improved  the  packaging  ratio  over  the  more  traditional 
use  of  tube  members. 

Section  8  summarizes  the  findings  of  the  study  and  offers  recommendations 
based  on  these  results.  While  many  refinements  of  the  design  are  suggested 
in  order  to  improve  the  spacecraft  overall  efficiency,  it  was  shown  that  the 
MMW  antenna  structure  can  perform  satisfactorily.  This  performance  can  be 
achieved  with  a  passive  structural  system  that  can  be  launched  by  a  single 
shuttle  flight. 

Before  proceeding  with  the  detailed  presentation,  it  is  helpful  to  con¬ 
sider  the  block  diagram  of  Figure  1-1.  It  emphasizes  the  major  possible  tasks 
and  the  way  they  are  interrelated  in  a  flow  towards  the  determination  of  the 
MMW  antenna  performance.  The  blocks  involving  active  control  techniques  were 
not  needed  since  a  satisfactory  performance  could  be  achieved  by  use  of  a 
passive  system. 
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SECTION  2 


ANTENNA  DESCRIPTION 


2.1  Antenna  Selection  and  Description 

The  selection  of  a  MMW  antenna  concept  was  based  on  the  trade-off 
studies  carried  out  by  the  Aerojet/NTL  team.  This  section  is  a  summary 
of  their  findings.  It  explains  the  merits  associated  with  the  millimeter 
wave  sensors,  their  mode  of  operation,  and  suggests  as  a  candidate  a  single 
surface  antenna  working  at  a  wavelength  of  X  =  3  mm  which  corresponds  to 
a  transmission  window  frequency  centered  at  95  GHz.  The  following  text  is 
partially  borrowed  from  the  AESC/NTL  Report. 

2.1.1  MMW  Antenna  Characteristics 


Enhancement  of  satellite  surveillance  can  be  derived  from  the  cloud 
and  weather  penetration  of  microwave /mil lime ter  wave  sensors.  Initial 
analyses  quantitatively  verified  the  feasibility  of  spacebound  passive 
microwave  technology.  The  effective  satellite-based  surveillance  of  stra¬ 
tegic  aircraft  and  naval  vessels  in  cloud-cover  and  weather  situations  can 
be  achieved  where  surveillance  is  presently  denied  by  the  limitations  of 
sensors  restricted  to  visual  and  infrared  wavelengths. 

Smoke,  dust,  fog,  and  clouds  frequently  prevent  satellite-based  obser¬ 
vation  of  surface  targets  wich  visual  and  infrared  sensors.  The  ability  to 
penetrate  these  cover  conditions  and  to  perform  equally  well  by  night  or 
day  makes  microwave  passive  sensing  an  attractive  candidate  to  significantly 
extend  the  capabilities  of  satellite  based  observations  beyond  their  present 
limitations . 


Passive  microwave  sensors  respond  to  millimeter  wavelength  thermal 
radiation  emitted  by  all  objects  in  proportion  to  their  physical  temperature 
Tfo  and  surface  emissivity  £  .  Useful  sensing  of  military  targets  results 
predominantly  from  high  contrasts  of  target  emissivity  with  respect  to  back¬ 
ground  surface. 

Oxygen  and  water  vapor  are  the  only  atmospheric  gasses  with  significant 
absorptions  at  millimeter  wavelengths.  These  combine  to  establish  the  clean 
air  transmission  windows  centered  at  frequencies  of  35  GHz,  95  GHz,  140  GHz, 
and  230  GHz. 


D  =  XRd 
a  c 


Design  of  a  microwave  pushbroom  sensor,  shown  in  Figure  2-1,  is  based 
on  optical  principles.  A  satellite  altitude  R  requires  an  antenna  dianeter 

^  to  achieve  a  surface  resolution  d  .  The  sensor  wavelength,  X, 

c 

in  millimeters  is  related  to  frequency,  f,  in  GHz  by  X  =  300  f  .  A  "push- 
broom"  of  N  contiguous  surface  resolution  cells  results  in  an  instantaneous 
crosstrack  coverage  swath  width  of  dimension  Wg.  As  a  consequence  of  the 

satellite  ground  velocity  v  ,  the  sensor  has  an  earth  surface  area  search 

.  S 

rate  A  =  W  x  v  . 

s  g 
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Pushbroom  microwave  satellite  imager. 


The  filled  pushbroom  sensor  is  configured  with  a  number  M  of  feed/ 
receiver  modules  which  sample  the  entire  crosstrack  swath  at  distance  inter¬ 
vals  equal  to  the  surface  resolution  cell  diameter  dc.  The  antenna  beam 
associated  with  the  individual  feed/receiver  modules  are  situated  at  fixed 
pointing  angles  in  the  crosstrack  plane.  Thus  satellite  motion  scans  each 

feed/receiver  module  to  sweep  along  a  continuous  inline  track  of  width  d  . 

c 

The  outputs  of  all  feed/receiver  modules  are  sampled  at  intervals  t.  =  d  v 
(the  dwell  time) .  1  ° 


Conceptually  the  satellite  based  passive  microwave  imager  resembles 
quite  closely  in  form  that  of  an  infrared  imaging  system.  A  reflecting 
system  focuses  the  radiation  of  individual  surface  resolution  cells  onto 
the  feed  apertures  situated  at  the  focal  surface. 

As  in  the  case  of  an  infrared  sensor  the  detected  output  from  the  micro- 
wave  receiver  output  varies  with  time,  having  a  magnitude  which  is  directly 
proportional  to  the  incoming  ladiaticn  intensity.  An  analog  signal  proces¬ 
sor  and  digital  data  processor  serve  the  same  functions  as  in  the  case  of 
an  infrared  mosaic  imaging  system. 

2.1.2  Description  of  Selected  Concept 


Figure  2-2  shows  the  general  form  of  the  reflector  antenna  selected. 

The  microwave  antenna  features  efficient  performance  and  the  reflector  sur¬ 
face  is  an  elliptic  torus.  The  feed  horns  are  deployed  on  a  cylindrical 
surface  passing  through  the  focal  circle  of  the  ellipses.  The  other  focal 
circle  lies  on  the  surface  of  the  earth.  The  design  of  the  feed  horns  pro¬ 
vides  compensation  for  the  aberration  caused  by  the  circular  crosstrack 
reflector  profile.  Each  horn  accepts  incident  energy  from  approximately 
half  of  the  reflecting  surface  and  provides  for  side  lobe  control  through 
proper  tapering  of  the  energy  accepted. 

The  elliptical  torus  reflector  is  not  a  perfect  collimator  so  "point 
source"  feed  horns  are  not  adequate  for  maintaining  pattern  gain  and  quality. 
Oversized  feeds  with  the  ability  of  controlling  the  amplitude  and  phase  over 
an  extended  focal  region  are  used  to  maintain  antenna  performance  levels. 
These  feeds  can  be  small  arrays  of  horn  radiators,  multimode  horn  feeds 
such  as  the  scalar  feed  or  corrugated  horn  feeds  or  horn/reflector  combina¬ 
tions. 

Tolerance  limits  must  be  imposed  on  the  positions  and  orientations  of 
the  antenna  feeds  relative  to  the  reflector  surface,  RMS  limits  on  surface 
roughness  and  surface  deviations  due  to  operational  disturbances  are  pro¬ 
vided  as  preliminary  requirements  for  satisfactory  performance.  Table  2-1 
summarizes  the  various  tolerance  limits. 


Table  2-1.  Tolerances  on  reflector  surface  and  feed 


X  =  3mm 

Primary  Reflector  Surface  RMS  Surface  Deviation  0.15mm 

RMS  Surface  Quality  0.05mm 

Feed  Tolerances  Torsional  1° 

Axial  1.0mm 

Decenter  0.15mm 

Pointing  Drift  50.0  yrad 

Accuracy  1.5  yrad 


2.2  Structural  Concept  and  Design 
2.2.1  Introduction  and  Objective 

The  task  of  the  structural  designer  is  to  provide  a  support  structure 
to  the  millimeter  wave  antenna  configuration  such  that  the  structural  system 
will  perform  within  tolerable  limits  when  exposed  to  orbital  environment. 

In  designing  a  space  structure  compatible  with  the  geometry  of  the  large 
optical  system,  it  is  important  to  generate  sufficient  stiffness  with  a 
minimum  amount  of  structural  mass.  Since  these  needs  are  basically  in  con¬ 
flict,  the  structural  design  will  emerge  as  a  result  of  trade-off  studies 
and  optimization  techniques.  Many  constraints  emanating  from  optical, 
structural  and  environmental  sources  cannot  be  ignored  in  the  optimization 
process.  In  addition,  the  size  of  the  shuttle  cargo  bay  will  impose  con¬ 
straints  on  the  structural  design  associated  with  the  specific  packaging 
and  deployment  technique. 

The  objective  of  the  CSDL  work  was  to  suggest  an  antenna  design  that 
uses  conventional,  state-of-the-art  structural  concepts  and  technology. 

More  precisely,  the  overall  structural  design  approaches  have  been  limited 
to  the  solution  space  involving  truss  construction  only.  Although  this 
traditional  concept  requires  a  relatively  larger  amount  of  mass  when  compared 
to  more  advanced  systems,  it  offers  the  advantages  of  simplicity  and  reli¬ 
ability  and  thus,  provides  a  desirable  level  of  confidence.  The  figure  of 
merits  of  this  design  and  its  performance  characteristics  will  be  contrasted 
with,  and  compared  with,  those  of  an  alternate  structural  concept,  i.e., 
an  electrostatically  shaped  membrane  reflector  proposed  by  another  contractor. 
Some  parameters  such  as  total  mass  and  packaging  ratio  are  easily  assessed, 
others  involving  control  complexity  and  ultimate  reliability  do  not  lend 
themselves  to  quantitative  measures  of  comparison. 
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2.2.2  Approach  to  Antenna  Structure  Design 


Designing  a  large  structure  that  performs  to  the  required  stringent 
optical  performance  is  challenging  and  the  task  is  further  complicated  by 
the  space  shuttle  launch  constraints.  It  was  stipulated  that  the  entire 
satellite  will  be  packaged  into  the  cargo  bay  of  a  single  shuttle.  This 
requirement  limits  the  designer  in  his  freedom  to  select  the  various  struc¬ 
tural  components,  their  size,  shape  and  dimensions.  A  third  tentative 
requirement,  that  the  antenna  have  a  satisfactory  performance  without  the 
help  of  active  means  of  control,  only  compounds  the  difficulties  in  an 
already  complicated  situation. 

The  large  size  of  the  reflector  dictated  segmentation  into  panels  for 
the  "rigid"  surface  structure.  In  order  to  take  best  advantage  of  the 
shuttle  cargo  dimensions  (4.57  a  in  diameter)  and  at  the  same  time  have  an 
efficient  assemblage  pattern,  a  panel  of  hexagonal  shape  with  a  diameter 
of  4.5  m  appeared  most  advantageous.  The  124  panels  are  supported  on  kine¬ 
matic  mounts  by  an  orthogonal  system  of  graphite-epoxy  trusses.  This  inter¬ 
face  between  panel  and  back-up  trusses  allows  for  structural  decoupling 
in  response  calculation  to  static  loadings,  while  the  choice  of  a  relatively 
high  fundamental  frequency  of  the  panels  should  permit  the  assumption  of 
mild  dynamic  coupling  between  the  panel  and  the  back-up  truss  assembly. 
Another  major  structural  component  consists  of  the  receiver  support  tower 
that  extends  some  30  m  above  the  reflector  surface.  Because  of  the  undesir¬ 
ability  of  optical  blockage,  this  structure  was  designed  in  such  a  way  that 
it  does  not  obstruct  any  incoming  or  outgoing  signals.  In  order  to  achieve 
an  acceptable  bending  stiffness,  the  tower  space  frame  construction  increases 
in  depth  from  the  top  down  to  extend  about  8  m  beyond  the  reflector's  edge. 
The  design  of  a  large  passive  reflector  that  performs  within  given  tolerances 
is  limited  by  state-of-the-art  technology  in  structures  and  materials. 

With  today's  technology,  the  ratio  between  reflector  diameter  and  surface 
accuracy  can  reach  approximately  10^  for  a  passive  structure  exposed  to  a 
typical  orbital  environment.  The  present  design  goes  beyond  this  number, 
since  the  radio  equals  10^  for  a  reflector  of  50  m  (equivalent)  diameter 
requiring  a  surface  accuracy  of  0.05  mm. 

In  order  to  achieve  this  goal,  many  avenues  capable  of  optimizing  the 
structure  have  been  explored.  In  this  attempt  the  overall  structure  was 
viewed  as  a  time- invariant  mechanical  system  acted  upon  by  environmental 
disturbance  inputs.  The  distortion  of  the  system  is  the  output  from  which 
the  antenna  performance  is  evaluated.  The  wide  range  of  options  at  the 
designer's  disposal  to  optimize  the  structural  system  must  be  exercised  sub¬ 
ject  to  existing  constraints.  Figure  2-3  attempts  to  list  such  options, 
many  of  which  have  been  used  to  minimize  the  structural  distortions  of  the 
back-up  truss  assembly.  A  special  effort  was  made  to  design  an  efficient 
thermal  blanket  over  the  reflector's  back-face  since  it  was  determined  that 
the  sun  and  earth  thermal  radiation  effects  translated  into  the  most  dele¬ 
terious  distortions.  The  thermal  model  of  this  structure  and  this  thermal 
shield  was  used  for  trade-off  studies.  The  results  are  discussed  in  detail 
in  Section  3.3. 
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In  the  selection  of  structural  materials,  whether  for  truss  members 
or  panel  egg-crate  construction,  graphite-epoxy  composites  invariably  emerged 
as  the  most  desirable  candidate.  This  choice  hinges  on  the  excellent  E/p 
ratio  (modujus  of  elasticity  over  density)  which  advantageously  affects  the 
stiffness  while  keeping  the  structural  mass  to  a  minimum.  At  the  same  time, 
graphite-epoxy  exhibits  a  very  3ow  coefficient  of  thermal  expansion  a, 
usually  assumed  to  be  0.2  x  10"^  m/m°K.  This  value  can  be  decreased  in 
one  privileged  direction  of  the  material,  if  the  proper  combination  of 
plys  can  be  achieved  in  its  construction. 

2.2.3  Antenna  Structure  Description 

The  millimeter  wave  antenna  selected  for  this  study  includes  a  large 
reflector  surface  in  the  shape  of  an  elliptical  torus  and  a  feed  horn  deployed 
on  a  cylindrical  surface  passing  through  the  focal  circle  of  the  ellipses. 

In  projection  the  reflector  is  approximately  rectangular  in  shape  with  side 
dimensions  of  60  m  and  30  m.  Since  the  area  of  the  reflector  employed  by 
each  receiver  is  a  circle  of  approximately  30  m  in  diameter  the  corner  areas 
of  the  rectangular  surface  are  not  needed  and  therefore,  have  been  removed. 
Figure  2-4  shows  a  projected  view  of  the  reflector  surface  and  illustrates 
the  arrangement  of  124  of  the  hexagonal  panels.  These  panels  are  each  sup¬ 
ported  by  three  kinematic  mounts  located  at  the  intersections  of  longitudinal 
and  transverse  back-up  trusses.  At  about  30  m  above  the  reflector  a  slender 
11  m  long  box  housing  the  feed  horn  equipment  is  supported  by  space  frame 
tower.  The  design  geometry  of  the  tower  structure  was  dictated  by  the  need 
for  non-obscuration  of  the  reflector.  For  that  reason  the  entire  tower 
structure  extends  beyond  the  reflector's  edge  and  its  triangular  shape  is 
needed  in  order  to  develop  the  bending  stiffness  required.  Its  attachment 
to  the  reflector  back-up  trusses  attempts  to  avoid  localized  reactions. 

Two  main  oversized  trusses  underneath  the  reflector  insure  a  gradual  inter¬ 
action  of  tower  and  back-up  trusses.  A  perspective  view  in  Figure  2-6 
illustrates  the  main  features  of  the  antenna  structures.  Details  of  rela¬ 
tive  dimensions,  elements  shapes,  sizes  and  material  properties  are  given 
below. 


The  entire  reflector  surface  is  provided  by  124  hexagonaliy  shaped 
graphite-epoxy  panels.  The  largest  dimension  of  the  panels  was  limited  to 
4.5  m  so  that  the  panel  can  fit  into  the  shuttle  cargo-bay  diameter  of 
4.57  m.  In  order  to  preserve  good  dynamic  characteristics  it  was  important 
to  design  a  relatively  stiff  and  therefore,  deep  panel.  Limitations,  however, 
were  imposed  by  the  packaging  requirements  and  the  choice  for  the  panel 
depth  was  finally  set  at  6  cm.  In  this  way  it  is  possible  to  stack  the 
124  panels  in  about  7.8  ra  length  of  the  shuttle  cylindrical  cargo.  The 
hexagonal  panel  is  an  egg-crate  construction,  using  thin  graphite- epoxy 
face  sheets  and  an  orthogonal  grid  of  webs.  It  weighs  just  under  4  kg/m^. 

To  minimize  the  mass,  the  back  face  sheet  and  the  webs  have  cut-outs  as 
shown  in  Figure  2-5.  The  three  kinematic  mounts  are  located  along  radial 
lines,  120°  apart,  and  preserve  the  symmetry  of  the  hexagon.  In  the  present 
layout  the  back-up  trusses  offer  hard-point  supports  at  the  apex  of  the 
hexagonal  panels,  but  it  might  be  desirable  to  move  these  supports  inward 
along  the  panel  radial  lines,  in  order  to  increase  the  panel  fundamental 
frequency.  Comparative  studies  have  been  made  and  are  reported  on  later 
in  this  report. 
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The  back-up  trusses  were  designed  as  an  orthogonal  set  of  1.5  m  deep 
trusses  and  spaced  such,  t’-at  their  intersections  coincide  with  the  comers 
of  hexagonal  panels.  There  are  9  longitudinal  trusses  at  3.375  m  on  center 
and  16  transverse  trusses  at  3.9  m  spacing.  Their  length  varies  and  some 
have  outriggers  in  order  to  accommodate  the  irregular  arrangement  of  the 
peripheral  panels.  Diagonal  members  in  the  top  and  bottom  faces  provide  the 
necessary  lateral  stability  of  the  assembly.  The  individual  truss  member's 
area  has  been  selected  so  that  the  bending  stiffness  of  the  platform  would 
be  sufficient  to  result  in  a  fundamental  natural  frequency  of  1  Hz.  The 
individual  truss  members,  usually  are  designed  with  the  traditional  thin- 
wall  tube  section,  because  this  section  is  most  efficient  in  terms  of  member 
bending  stiffness,  natural  frequency,  and  stability.  In  this  design  an 
open  thin-wall  section  was  suggested  since  it  offered  the  advantage  of  very 
efficient  packaging  Lor  the  trusses.  Because  of  the  angle  shape  of  all  its 
constitutive  members,  a  truss  can  be  nested  into  another  one  by  adding  to 
the  package  only  slightly  more  than  the  0.2  cm  of  its  section's  wall  thickness. 
The  packaging  ratio  is  approximately  10  times  the  one  corresponding  to  a 
similar  packaging  of  trusses  using  tube  members.  The  stability  study  of  a 
typical  truss  member  showed  that  it  would  require  a  relatively  high  stress 
(or  equivalently,  axial  force)  to  reach  the  calculated  buckling  load.  Under 
normal  conditions  the  existing  forces  are  anticipated  to  be  well  below  this 
value  and  the  design  was  deemed  satisfactory.  The  graphite-epoxy  truss 
assembly  lias  a  mass  of  approximately  1  kg/nr. 

The  30  meter  high  tower  supporting  the  receiver  package  was  designed 
to  impart  to  the  structure  a  10  Hz  natural  frequency  in  lateral  bending. 

As  it  avoids  obscuration  for  the  reflector,  the  cantilevered  tower,  11.7  m 
in  width,  extends  beyond  the  reflector's  edge  in  order  to  develop  the  neces¬ 
sary  bending  stiffness.  Its  triangular  shape  has  a  maximum  of  8  meters  at 
its  base  where  it  meets  with  the  antenna  surface.  The  individual  members, 
longer  than  those  of  the  back-up  trusses,  arc  graphite-epoxy  thin  wall  Lubes, 

10  cm  in  diameter  by  0.15  cm  in  thickness.  Figure  2-6  illustrates  the  struc¬ 
tural  arrangement  of  the  tower  and  its  attachment  to  the  back-up  trusses' 
assembly . 

Although  it  is  anticipated  that  a  substantial  amount  of  equipment  is 
needed  for  the  satellite  operations  such  as  data  processing,  communication, 
power  generation,  it  was  too  early  to  assign  extra  masses  to  the  structure. 

The  300  kg  assigned  to  the  tip  of  the  tower  is  larger  than  the  mass  required 
for  the  feed-receiver  equipment  and  is  viewed  as  .including  some  of  the  extra 
equipment.  The  location  of  this  mass  is  not  favorable  in  terms  of  dynamic 
responses  and  should  lead  to  conservative  results. 

2 . 2 . A  Masses  and  Inertia 

The  summary  of  the  masses  estimated  for  this  design  is  as  follows: 


Reflector  Panels  at  4  kg/m^ 

=  6300  kg 

2 

Back-up  Trusses  at  1  kg/m“ 

=  1575  kg 

Feed  Support  Structure 

625  kg 

Equipment  in  Feed 

=  300  kg 

TOTAL 

=  8800  kg 
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igure  2-6.  Structural  arrangement  of  antenna  components. 
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The  mass  moments  of  inertia  have  been  computed  with  respect  to  the 
reference  x,y,z  system  of  axes  located  at  the  center  surface  point  of  the 
reflector.  Since  the  overall  structure  is  not  completely  symmetrical, 
principal  axes  of  inertia  have  been  determined  (at  the  center  of  mass)  and 
the  associated  principal  mass  moments  of  inertia  have  been  evaluated. 
Figure  2-7  shows  the  orientation  of  x,y,z  system  and  the  principal  axes 
of  inertia  (1,2,3)  at  the  center  of  mass.  The  values  of  the  various  mass 
moments  of  inertia  appear  in  the  same  figure. 
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TOTAL  MASS  =  8800  KG 


LOCATION  OF  C,M. 


x  =  0.0 
y  =  -1.12m 
z  =  4.10m 


•  ABOUT  CENTER  OF  MASS 

I XX  =1.154  106  KG-M2 
IYY  =  2.274  106  KG-M2 
IZZ  =  2.715  106  KG-M2 
IXY  =  0 

IYZ  =  7.718  104  KG-M2 
ixz  =  0 

•  PRINCIPAL  MOMENTS  OF  INERTIA 
W.R.T.  CENTER  OF  MASS 

lX  =  1.154  106  KG-M2 
I2  =  2.262  106  KG-M2 
t  =  2.728  106  KG-M2 

0  =  9.6° 


Figure  2-7.  Mass  moments  of  inertia. 
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SECTION  3 


DISTURBANCE  IDENTIFICATION 


3.1  Introduction  and  Approach 

Disturbances  affecting  th°  structural  performance  of  the  MMW  antenna 
system  derive  from  various  external  and  internal  sources.  The  external 
sources  of  disturbance  consist  of  solar  and  terrestrial  related  effects. 

The  resultant  torques  and  forces  applied  to  the  structure  are  functions  of 
the  satellite's  position  in  its  orbit,  and  its  size,  shape,  and  mass  distri¬ 
bution.  While  the  identification  and  quantification  of  these  external  load¬ 
ings  can  be  analyzed  through  well  understood  derivations  and  available  pro¬ 
grams,  the  internal  loadings  generated  by  satellite  components  are  more 
difficult,  at  this  point,  to  identify  and  evaluate.  They  might  include 
torques  and  forces  from  reaction  wheel  mass  imbalances  and  bearings,  drive 
mechanisms,  solar  array  drives  and  compensating  jet  thrusters.  Since  the 
satellite  design  is  yet  to  be  finalized,  many  of  its  mechanical  part  charac¬ 
teristics  have  not  been  identified.  As  an  example,  the  effects  of  the  movable 
solar  panels  and  communication  antennas  can  only  be  estimated. 

From  a  structural  analyst's  point  of  view,  it  is  desirable  to  categorize 
the  disturbances  by  their  frequency  content.  Once  the  dynamic  characteristics 
of  the  structure  have  been  determined,  the  various  loadings  can  usually  be 
classified  into  quasi-static  and  dynamic  (vibrational)  disturbances.  In 
the  first  case,  a  static  analysis  is  sufficient  to  predict  the  structural 
distortions,  whereas  a  more  involved  dynamic  analysis  is  required  for  evaluat¬ 
ing  responses  to  dynamic  disturbances.  Table  3-1  attempts  to  identify  and 
classify  under  quasi-static  and  dynamic  headings  most  of  the  loadings  expected 
to  ultimately  affect  the  structural  performance. 

Since  this  study  is  mainly  concerned  with  the  structural  behavior  of 
the  satellite,  it  appears  natural  at  first  to  identify  only  a  few  of  these 
disturbances  that  are  most  likely  to  be  critical  and  will  furnish  meaning¬ 
ful  response  information.  In  other  words,  this  analysis  does  not  pretend 
to  offer  exact  responses  to  all  exactly  identified  disturbances,  but  it  will 
evaluate  the  effects  of  some  selected  representative  loadings  on  the  struc¬ 
tural  performance  in  a  behavioral  sense.  While  the  disturbances  are  usually 
approximated,  the  structure  will  be  faithfully  modelled  so  as  to  make  it 
possible  to  assess  the  viability  of  the  structural  concept. 

With  this  approach  in  mind,  some  dominant  disturbances  arising  during 
the  orbit  period  of  approximately  1.56  hours  have  been  established.  In 
the  following  sections,  the  procedures  used  in  evaluating  the  effects  of 
solar  and  earth  thermal  fluxes,  atmospheric  drag,  solar-radiation-pressure 
and  gravity  gradient  are  discussed.  The  approaches  usually  employ  simplified 
models  and  the  assumptions  used,  as  well  as  the  effects  neglected,  are  care¬ 
fully  defined,  so  that  the  results  can  be  meaningfully  interpreted  and  dis¬ 
cussed. 
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Table  3-1.  Disturbances  identification 

SOURCES 

•  External  or  environmental 

•  Configuration  dependent 

•  Attitude/position  system 

•  Cryo-system  (if  required) 

GROUPED  BY  FREQUENCY  CONTENT 
QUASI-STATIC 

•  Thermal  loadings  (steady  state)  from  various  sun  angles,  eclipsing,  earth  radiation 

•  Solar  &  earth  radiation  pressures 

•  Gravity  gradient 

•  Solar  panel  motion  reaction  torques 

•  Geomagnetic  field 

•  Solar  wind 

•  Atmospheric  drag 

DYNAMIC 

•  Thrusters  transients  for  attitude  control,  station  keeping 

•  Thermal  shock  (transient)  from  eclipsing 

•  Solar  array  drive 

•  Communication  antenna  motion 

•  Refrigeration  machinery 

•  Coolant  mass  shift 

•  Turbulence  in  pipes 
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3.2  Orbital  Thermal  Disturbances 


3.2.1  Introduction 

A  satellite  in  earth  orbit  is  exposed  to  internal  heat  sources,  as  well 
as  to  several  thermal  radiation  fluxes  emanating  from  different  external 
sources  While  the  solar  flux  is  the  predominant  one,  its  reflection  from  the 
earth,  known  as  albedo,  and  the  IR  flux  irradiated  from  a  warm  earth  may  also 
be  significant.  Incident  flux  levels  on  the  satellite  vary  as  a  function  of 
orbit  geometry,  satellite  orientation,  and  time.  Such  variations  give  rise  to 
time-varying  temperature  distributions,  which  in  turn,  cause  time-varying  de¬ 
formations  of  the  structure.  It  is  the  task  of  the  analyst  to  make  suitable 
assumptions  in  building  a  mathematical  model  capable  of  adequately  predicting 
the  thermal  behavior  of  the  spacecraft.  From  this  model,  temperature  distri¬ 
butions  which  are  a  function  of  both  spatial  coordinates  and  time  are  predicted 
throughout  the  structure. 

Since  the  orbital  thermal  environment  and  the  satellite  temperature  re¬ 
sponses  thereto  are  time-varying,  thermal  storage  effects  are  generally  very 
significant.  Therefore,  suitable  mathematical  modeling  requires  the  inclusion 
of  thermal  storage  capacitances  for  each  subdivision  of  the  satellite  structure. 
The  degree  of  fineness  in  subdividing  the  structure,  necessary  for  adequate 
modeling  accuracy,  is  an  important  trade-off  and  is  an  item  requiring  experi¬ 
ence  and  intuitive  judgement  of  the  analyst. 

As  the  satellite  experiences  the  various  external  radiation  fluxes,  three 
modes  of  thermal  exchange  can  take  place  between  its  subdivisions  and  to  space: 
conduction,  convection  and  radiation.  In  space  applications,  gas  convection 
and  conduction  between  structural  components  is  negligible  because  of  the  rare¬ 
fied  atmosphere  and,  therefore,  is  disregarded. 

Solid  conduction  is  easily  included  in  thermal  models,  but  its  contribu¬ 
tion  to  temperature  changes  depends  largely  on  the  particular  problem.  It  is 
found  to  be  very  small,  for  example,  along  typical  truss  members  because  of  the 
thinness  of  their  wall  sections.  In  the  design  of  high  precision  mirrors,  con¬ 
duction  becomes  a  significant  factor  because  it  establishes  thermal  gradients 
across  the  plate  thickness  and,  thus,  has  a  predominant  contribution  to  the 
surface  distortions.  In  some  cases  conductance  couplings  parallel  to  the 
mirror  surface  are  negligible,  while  in  other  cases  they  are  significant. 

By  far  the  most  important  mode  of  thermal  exchange  for  the  large  space 
structure  is  radiation.  A  satellite  structure  can  be  viewed  as  an  aggregate 
of  members  and  surfaces  that  are  exposed  to  external  incident  radiant  energy 
but  at  the  same  time  emit  radiant  energy.  All  these  surfaces  interact  with 
each  other  to  some  degree  which  for  typical  diffuse  surfaces  is  accounted  for 
bv  the  use  of  gray-body  radiation  exchange  view-factors.  Special  methods  are 
used  for  specular  surfaces.  As  the  satellite  progresses  along  its  orbital  path, 
the  angles  of  incidence  of  the  incident  fluxes  are  changing  constantly  (in¬ 
cluding  shadowing  effects)  while  individual  surfaces  usually  remain  in  constant 
orientation  with  respect  to  one  another.  There  exists  a  number  of  thermal  an¬ 
alyzer  computer  programs  capable  of  generating  thermal  models  and  calculating 
the  temperatures  for  any  predefined  orbit.  CSDL  uses  a  combination  of  pro¬ 
grams  shown  in  Figure  3-1. 
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ORBIT  1  CSDL  MODEL  DEFINITION 


Since  the  characteristics  of  each  orbit  cycle  change  relatively  slow 
with  time,  it  is  expected  that  the  thermal  disturbance  will  become  cyclic. 

After  several  revolutions  the  initial  transient  effects  will  disappear  and  the 
temperature  histories  of  satellite  surfaces  will  become  almost  identical  for 
each  subsequent  orbit.  This  quasi-steady-state  cyclic  condition  usually  pro¬ 
vides  a  good  measure  of  the  expected  thermal  fluctuations,  and  will  give  rise 
to  slowly  changing  mechanical  responses  which  are  calculated  using  static 
structural  analysis.  Under  special  conditions,  however,  the  effect  of  rapidly 
changing  thermal  conditions  must  be  calculated  by  a  dynamic  mechanical  analysis 
method.  This  could  be  created,  for  example,  by  the  sudden  passage  from  total 
eclipse  to  sunlight  incident  to  surfaces  of  low  modal  frequencies. 

For  this  specific  study  of  the  MMW  satellite,  many  simplifying  assumptions 
were  made  in  order  to  obtain  feasibility  assessment  of  the  design,  while  mini¬ 
mizing  costly  detailed  thermal  modeling. 

Because  the  structure  is  composed  of  thin  sheet  materials  of  large  surface 
areas,  or  of  thin-walled  truss  tubing  of  large  lengths,  thermal  solid  conduc¬ 
tances  along  the  surfaces  or  truss  element  axes  are  negligible  in  comparison 
to  radiative  couplings,  and  were  ignored.  The  repetitive  design  pattern  of 
the  antenna  structure  allowed  simplified  one-dimensional  analyses  which  ignore 
edge  effects  (these  effects  can  be  minimized  by  suitable  edge  shielding). 

Also,  truss  members  were  treated  as  being  isothermal  through  their  cross- 
sections  at  any  particular  axial  location. 

The  orbital  thermal  inputs  were  for  a  50U-km  altitude  equatorial  orbit, 
with  the  satellite  oriented  with  the  antenna  beam  always  directed  toward  the 
center  of  the  earth.  Orbital  geometry  was  simplified  by  the  worst-case  assump¬ 
tion  that  the  solar  vector  rotated  in  the  plane  of  a  principal  axis  of  the 
satellite.  The  orbit  period  is  1.58  hours,  of  which  38%  of  the  time  the  sun 
is  eclipsed  by  the  earth.  For  this  study  the  combined  earth  IR  and  albedo  flux 
was  taken  as  a  constant  value,  intended  only  as  an  aid  in  achieving  approximate 
absolute  temperature  levels.  The  disturbances  due  to  earth  flux  variations 
are  minor  in  comparison  to  those  due  to  direct  solar  loading.  This  constant 

2 

flux  incident  to  the  satellite  was  taken  to  be  344  watts/m  ,  while  the  solar 

2 

incident  flux  was  taken  as  1350  watts/m  . 

3.2.2  Antenna  Panels  and  Trusses 
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3. 2. 2.1  Model  Definition 


Thermal  modeling  of  the  antenna  panels  and  trusswork  was  based  on  the 
general  arrangement  of  Figure  3-2.  The  radiative  effects  of  cross  members  be¬ 
tween  nodes  4  and  5  of  the  trusswork  were  ignored,  and  the  reflector  panels 
were  treated  as  being  continuous  over  the  whole  antenna.  By  ignoring  curva¬ 
ture  and  edge  effects  for  the  whole  antenna,  the  treatment  was  vastly  simpli¬ 
fied  to  the  one-dimensional  7-node  problem  shown.  Curvature  effects  are  ap¬ 
proximately  accounted  for  separately,  and  edge  effects  can  be  minimized  by 
suitable  shielding. 

The  thermal  blanket  is  composed  of  top  (7)  and  bottom  (6)  layers  of 
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Figure  3-2.  Antenna  panel  and  truss  thermal  model. 


0.1  nan  thick  teflon  or  other  fabric  covers,  separated  by  alternating  layers  of 
dacron  "bridal  veil"  separators  and  0.01  mm  thick  aluminized  mylar  shields 
(possibly  5  layers) .  The  exterior  surface  of  the  top  cover  is  metalized  to 
achieve  reasonably  low  solar  absorptivity  and  emissivity  (this  could  be  a 
radiatively  selective  surface  for  temperature  level  adjustment).  In  some 
cases  analyzed,  the  blanket  shielding  layers  were  "short-circuited",  essen¬ 
tially  making  nodes  6  and  7  a  single  node,  to  represent  a  simpler  and  less- 
costly  single  sheet  of  material. 

A  number  of  different  design  alternatives  were  studied  for  the  reflector 
panels.  The  alternatives  are  described  in  their  variations  from  a  baseline 
design  used  for  final  analysis.  This  baseline  design  consists  of  a  sandwich 
panel,  all  made  of  graphite-epoxy.  The  reflector  sheet  (2)  is  1.5  mm  thick, 
while  the  back  sheet  (3)  is  1  mm  thick.  The  core  is  comprised  of  webs  of  1  mm 
thickness  spaced  in  a  square  matrix  on  88  mm  centers.  Web  graphite  fiber 
orientation  is  assumed  normal  to  panel  sheets,  for  maximum  thermal  conductance 
between  sheets.  (The  design  variations  are  shown  in  Figure  3-3). 

3. 2. 2. 2  Design  Trade-offs 

Preliminary  design  trade-off  studies  were  made,  results  of  which  are  sum¬ 
marized  in  Figure  3-3.  For  each  design  condition,  results  are  shown  as  three 
different  temperature  difference  values,  pertaining  to  the  panels  or  the  truss- 
work.  The  values  "MAX  ATG"  and  "MIN  ATG"  are  the  maximum  and  minimum  cyclic 
gradient  temperature  differences  experienced  across  the  panels  or  trusswork. 

The  node  notation  shown  defines  the  calculation.  For  example,  nodes  "2-3" 
indicates  the  temperature  at  node  3  is  substracted  from  the  temperature  at  node 
2.  The  values  "AVG  ATU"  are  the  amplitudes  of  the  cyclic  temperature  excur¬ 
sions  experienced  by  the  panels  or  trusswork.  The  values  are  for  the  excur¬ 
sion  of  the  average  temperatures  of  the  two  associated  nodes;  for  example, 
the  average  of  nodes  2  and  3. 

For  the  comparison  of  various  reflector  panel  designs  it  was  desired  to 
subject  them  to  a  large  disturbance,  and  therefore  several  simulations  were 
run  without  solar  eclipse  and  without  an  insulation  blanket.  Runs  2,  3,  6  and 
7  show  these  comparisons.  It  will  be  shown  later  that  for  small-diameter 
panel  cells  the  ATG  values  are  the  significant  distortion  influences,  and  so 
for  a  given  material  the  smallest  ATG  values  result  in  the  least  distortion. 
Since  the  expansion  coefficient  for  aluminum  is  much  greater  than  for  graphite- 
epoxy.  Run  6  is  not  a  good  candidate  on  the  basis  of  relative  distortions. 
Therefore,  the  Run  3  design  was  chosen  as  the  baseline  design. 

There  is  obviously  a  trade-off  between  panel  weight  and  thermal  distor¬ 
tion,  as  indicated  by  Runs  2  and  3.  Another  trade-off  is  apparent  for  the 
real  orbital  mission  (with  solar  eclipse)  in  comparing  Runs  4  and  5.  Thermal 
distortions  of  both  panels  and  trusswork  may  be  traded  for  cost  and  weight  of 
the  insulation  blanket.  It  was  decided  to  use  Run  4  as  the  case  to  analyze 
in  further  detail,  to  see  if  distortions  would  be  acceptably  low  without  the 
insulation  blanket. 

3. 2.2. 3  Baseline  Run  4  Results 

The  timelines  of  orbital  temperatures  for  Run  4  are  shown  by  Figure  3-4. 
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TEMPERATURES  (°K) 


RUN 

ECLIPSE 

BLANKET 

PANEL  TYPE 

1 

NO 

YES 

Half  the  web 
mass  and 
conductance 

NODES 


MAX  ATG 
MIN  ATG 
AVG  'ATU 


2 

NO 

NO 

Half  the  web 
mass  and 
conductance 

NO 

NO 

Baseline 

4 

YES 

NO 

Baseline 

5 

YES 

YES 

Baselinp 

6 

NO 

NO 

Material  all 
aluminum 
(equal  mass) 

■ 

NO 

NO 

Air-mat 

Design 
(equal  mass) 

2-3 


+5.1 

-2.4 

22.3 


4-5 


+0.2 

-0.1 

15.0 


MAX 

ATG 

+8.4 

+0.9 

MIN 

ATG 

-4.3 

-0.6 

AVG 

ATU 

10.9 

12.7 

MAX 

ATG 

+0.9 

MIN 

ATG 

-0.6 

AVG 

ATU 

10.8 

11.6 

MAX 

ATG 

+1.8 

+0.6 

MiN 

ATG 

-2.7 

AVG 

ATU 

15.6 

17.4  | 

MAX 

ATG 

+0.7 

+0.1 

MIN 

ATG 

-0.4 

-0.1 

AVG 

ATU 

2.6 

0.7 

MAX 

ATG 

+1.8 

mmm 

MIN 

ATG 

-0.9 

AVG 

ATU 

10.8 

11.6  | 

MAX 

ATG 

+9.2 

+0.9 

MIN 

ATG 

-7.0 

-0.6 

AVG 

ATU 

10.8 

11.6 

Figure  3-3.  Preliminary  design  trade-offs. 
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Figure  3-4.  Temperatures  at  center  of  antenna  (run  4) 


The  time  scale  is  arbitrary,  representing  the  fourth  orbit  cycle  from  an  iso¬ 
thermal  initial  condition.  As  may  be  seen  by  comparison  of  temperatures  at 
-90°  and  +270°  of  sun  angle,  the  orbit  is  virtually  in  quasi-equilibrium. 

It  is  of  interest  to  observe  the  perturbations  that  occur  at  the  beginning  and 
end  of  eclipse,  where  the  sun  shines  on  the  reflector  panels  for  a  short  time. 
This  may  be  compared  to  Run  3  results,  shown  by  Figure  3-5,  which  shows  the 
effects  of  much  greater  reflector  panel  solar  exposure  when  there  is  no  eclipse. 

While  these  results  are  for  the  center  region  of  the  antenna,  based  on 
the  assumption  of  a  flat  plane,  the  real  antenna  curves  through  an  arc  of 
±30  degrees  from  the  center,  in  the  x-plane  (60  m  dimension  of  the  antenna)* 

An  approximation  was  made  of  the  effects  of  this  curvature  by  adjusting  the 
Run  4  temperature  results  by  introducing  phase  angle  off-sets  in  the  results. 
This  was  done  for  four  zones,  the  center  of  which  were  at  angles  of  +22.5, 

+7.5,  -7.5  and  -22.5  degrees  from  the  center.  In  addition,  corrections  were 
made  in  the  reflector  panel  temperatures  to  account  for  the  geometry  and 
shadowing  effects  of  solar  loadings. 

These  curvature  adjustments  were  made  with  the  Z-axis  at  three  sun  angles 
of  -92,  +68  and  +118  degrees,  estimated  as  being  the  worst  thermal  distortion 
excursion  possibilities,  for  which  distortion  calculations  would  be  made. 

These  results  are  tabulated  in  Figure  3-6.  New  nodal  definitions  are  indica¬ 
ted  in  the  sketch  to  identify  the  four  zones.  Distortion  results  based  on 
this  data  are  reported  in  Section  4.5. 

3.2.3  Tower  Structure 

3.2. 3.1  Model  Definition 

The  accurate  thermal  modeling  of  complex  truss  structures  that  are  di¬ 
rectly  exposed  to  solar  loadings  is,  in  general,  tedious  and  costly.  For  a 
first-order  estimate  of  thermal  distortions,  a  simplified  model  was  defined 
for  thermal  analysis.  Since  the  side  tower  truss  sections  (in  the  Y-Z  plane) 
are  parallel  it  was  decided  to  assume  a  solar  orbit  in  this  plane  so  that  one 
side  section  could  be  modeled,  with  all  shadowing  by  truss  members  in  that 
plane  and  with  the  other  side  responding  identically  by  symmetry.  The  modeling 
in  other  solar  orientations  would  be  more  difficult,  but  would  generally  be 
expected  to  produce  lower  total  distortion  excursions  than  for  this  case  where 
the  solar  shadowing  effects  are  severe  in  causing  tower  bending. 

The  thermal  model  is  shown  by  Figure  3-7.  Thermal  storage  capacitances 
are  included  in  this  model,  but  truss  member  axial  solid  conductances  are 
assumed  zero  because  they  are  very  raall  relative  to  radiative  couplings  to 
space  (members  are  assumed  isothermal  in  cross-section).  Also,  due  to  the 
high  surface  reflectivities  and  generally  small  view  factors  between  truss 
members,  there  are  negligible  radiative  couplings  between  members.  Each  member 
was  therefore  treated  independently,  having  only  direct  solar  loading  timelines 
and  coupling  to  space.  Solar  loadings,  including  shadowing  effects,  were 
calculated  in  detail  for  each  node.  However,  the  truss  members  shown  by 
dotted  lines  were  accounted  for  in  their  shadowing  effects,,  but  not  included 
in  the  thermal  model.  It  was  assuioed  cnat  the  major  impact  on  tower  bending 
was  due  to  temperature  variations  of  the  perimeter  members  modeled.  A  more 


-92 


+68 


+118 
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287.3 
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22 
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32 
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23 
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33 
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14 
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24 
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34 
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44 
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283.5 

285.1 

15 

300.6 

286.7 

282.7 

25 

300.5 

285.3 

282.5 

35 

299.1 

284.0 

283.3 

45 

297.2 

282.8 

285.1 

significant  but  still  minor  error  was  introduced  by  omission  of  radiative 
interaction  between  the  truss  members  and  the  antenna  (which  also  has  high 
reflectivity) . 

3 . 2 . 3 . 2  Tower  Temperature  Results 

The  temperature  histories  for  the  tower  truss  members  are  shown  by  Figure 
3-8.  Although  the  simulation  was  run  through  four  orbit  cycles,  it  is  evident 
that  temperature  levels  are  still  decreasing  (compare  -90  and  +270  degree 
values).  However,  these  results  give  a  good  approximation  to  the  cyclic 
temperature  excursions  in  quasi-equilibrium. 

It  is  evident  from  Figure  3-8  that  the  thermal  response  is  quite  complex, 
this  mainly  being  a  result  of  the  shadowed  solar  loading  timelines.  Obviously, 
responses  would  vary  noticeably  with  orbital  orientation,  and  these  results 
must  be  considered  a  first-order  estimate.  Data  from  Figure  3-8  was  used  to 
calculate  tower  deformations  at  a  number  of  different  times,  these  being  pre¬ 
sented  in  Section  5.2. 

3.2.4  Discussion 

The  thermal  modeling  work  presented  here  was  intended  for  or der-of -mag¬ 
nitude  estimation  of  orbital  thermal  distortions,  to  assess  design  feasibility 
and  problems.  It  is  advisable  that  the  thermal  modeling  work  be  more  compre¬ 
hensive,  in  any  future  detailed  design  effort. 

The  results  for  the  assumed  baseline  design  antenna  and  associated  truss- 
work,  and  for  the  tower  trusswork,  assume  rather  simple  and  passive  thermal 
hardware.  With  negligible  cost  and  weight  penalties  for  thermal  design,  the 
resulting  distortions  presented  in  Section  5.2  are  acceptable. 

However,  it  would  be  quite  easy  to  improve  the  thermal  design  further  by 
passive  means,  of  low  cost  and  weight  penalty.  For  example,  it  may  be  quite 
feasible  to  achieve  lower  surface  emissivities  on  surfaces  exposed  to  solar 
loadings.  Also,  it  is  obvious  (from  comparision  of  Runs  4  and  5)  that  great 
reduction  of  antenna  distortions  would  result  from  use  of  an  insulation  blanket 
behind  the  antenna  trusswork.  A  similarly  large  reduction  of  distortions 
would  occur  at  the  tower  trusswork  if  it  were  enclosed  in  an  insulation  blanket, 
thereby  transforming  the  large  point-source  solar  thermal  loadings  into  at¬ 
tenuated  diffuse  loadings. 

3.3  Atmospheric  Drag 

3.3.1  Introduction 

The  MMW  satellite  in  orbit  at  a  height  of  500  km  encounters  significant 
aerodynamic  forces  as  it  flies  through  the  rarefied  upper  atmosphere.  These 
forces  create  orbital  perturbations  that  must  be  carefully  evaluated.  As  a 
result  of  the  braking  effect  of  the  atmosphere,  the  satellite  experiences  an 
orbital  decay  and  usually  a  deterioration  of  its  pointing  accuracy.  Although 
these  perturbing  forces  are  so  small  that  the  spacecraft  structure  is  not 
appreciably  affected  in  its  configuration,  their  overall  effect  must  be 
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counteracted  by  an  attitude  control  system.  It  is  the  disturbance  associated 
with  that  control  law  that  must  be  carefully  designed  and  tailored  so  as  to 
minimize  the  dynamic  structural  responses  of  the  delicate  spacecraft. 

This  section  briefly  reviews  the  important  factors  entering  the  calcula¬ 
tion  t  f  the  upper  atmosphere  density  used  in  the  evaluation  of  drag  forces  and 
torques.  A  brief  study  of  the  orbital  decay  and  the  MMW  antenna,  performed 
by  a  special  purpose  computer  program  is  also  presented  for  two  levels  of 
solar  activity.  It  emphasizes  the  problems  associated  with  station-keeping 
and  indicates  the  propulsion  requirements  to  restore  the  orbit.  Finally  a 
simplified  calculation  of  force  and  torque  resultants  with  respect  to  the 
center-of-mass  of  the  satellite  is  presented  for  some  extreme  cases  of  two 
atmospheric  densities  corresponding  to  an  average-  and  a  high-level  of  solar 
activities. 

3.3.2  Upper  Atmosphere  Density  Model 

The  determination  of  the  upper-atmospheric  density  is  usually  based  on  the 
study  of  the  orbital  decay  of  a  large  number  of  earth  satellites.  The  aim  is 
to  define  a  simple  analytic  relationship  for  determining  the  upper-atmospheric 
density  allowing  for  its  main  variations,  such  as  diurnal  and  semi-annual 
variations  as  well  as  its  correlation  with  the  intensity  of  the  solar 
activity  at  the  wavelength  of  10.7  cm.  The  density  p  can  be  computed  as  a 
product  of  several  factors  each  of  them  corresponding  to  one  or  other  forms 
of  the  known  density  variations: 

p  =  PH  K1  K2  K3  K4* 

where  PH  is  a  night-time  vertical  profile  of  the  atmospheric  density;  is 

a  factor  describing  the  density  change  with  the  change  of  solar  activity  F 
at  10.7  cm  relative  to  some  average  level  of  radiation  F^;  is  a  coefficient 

considering  the  diurnal  effect  in  the  density  distribution;  is  a  correction 

for  the  semi-annual  effect;  is  a  factor  taking  into  account  the  correlation 

between  the  atmospheric  density  variations  and  geomagnetic  perturbations. 

The  coefficients  change  with  a  change  of  the  mean  level  of  solar  activity 
Fq  over  the  eleven-year  cycle  of  solar  activity.  Other  density  variations 
are  not  taken  into  account  since  they  have  a  smaller  influence.  From  the  ob¬ 
servation  of  the  motion  of  existing  earth  satellites,  the  available  experi¬ 
mental  material  is  usually  divided  in  groups  corresponding  to  mean  levels  of 

solar  activity  F  of  75,  100,  125,  150,  200,  250  units  (10-22  Wnf2  Hz-1). 

All  the  K  factors  above  are  described  by  involved  analytical  relationships. 
With  the  atmospheric  density  information  and  the  definition  of  the  orbit  param¬ 
eters,  the  orbital  dynamics  and  the  station  keeping  requirements  of  the  MMW 
antenna  are  evaluated  with  the  help  of  a  specialized  computer  program.  The 
following  study  employed  the  current  semi-analytical  satellite  theory  imple¬ 
mentation  in  the  GTDS  R/D  program  (Goddard  Trajectory  Determination  System 
(GTDS),  Research  and  Development  version). 
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Orbital  initial  conditions  and  spacecraft  (S/C)  parameters  were  assumed 
as  follows: 


Initial  Conditions: 

a  =  6878.0  km 
e  =  0.001 
i  =  100° 

=  0.0° 

=  0.0° 

M  =  0.0° 

epoch  =  January  31,  1979 

S/C  Parameters: 

A  =  642 . 0 
M  =  8800.0  kg 


The  orbital  conditions  provide  a  500  km,  circular  orbit.  The  orbital  incli¬ 
nation  of  100°  is  designed  to  provide  global  coverage  and  to  be  consistent 
with  flight  operations  from  WTR.  The  value  2.0  of  the  coefficient  of  drag 
(Cj-j)  corresponds  to  the  default  option  in  the  GTDS  R/D  program. 

Two  flight  conditions  have  been  studied.  The  first  case  takes  the  above 
data  and  assumes  a  Harris/Priester  atmospheric  density  model  corresponding 
to  a  moderate  level  of  solar  activity  (?1()  ?  =  150).  This  represents  a  typi- 

cei  flight  condition  over  the  eleven  year  solar  cycle.  The  second  case  assumes 

an  atmospheric  density  model  with  a  high  level  of  solar  activity  (F  =  250). 

10  *  7 

This  high  level  of  solar  activity  better  represents  the  conditions  expected 
in  1979  and  1980. 

3.3.3  Orbital  Decay 


semi-major  axis 
eccentricity 
inclination 
ascending  node 
argument  f  perigee 
mean  anon  .  ly 


For  the  average  solar  activity  (F^q  ^  =  150),  the  key  elements  of  the 
orbital  history  are  summarized  in  Table  3—2.  The  semi— major  axis,  a. 

Table  3-2  ?10  ?  =  150 


T(days) 

a  (km) 

PH(perifocal) 
height, km 

0 

6878. 

492.982 

10 

6872.792 

493.233 

20 

6867.475 

482.476 

30 

6861.969 

471.230 

is  decreasing  in  a  nearly  linear  manner  due  to  the  atmospheric  drag.  The 
perigee  height  is  also  decreasing  due  to  the  decay  in  a.  However,  there  is 
also  a  long  periodic  oscillation  due  to  gravitational  perturbations  of  the 
orbital  eccentricity  (this  can  be  seen  more  clearly  in  the  GTDS  print-out). 

Over  the  thirty  day  period, 

Aa  =  -16.031  km 
APH  =  -21.75?  km 

The  average  daily  Aa  =  -0.5344  km. 

For  the  high  solar  activity  (F^  ^  =  250) ,  the  key  elements  of  the  orbital 
history  are  summarized  in  Table  3-3.  The  primary  difference 

Table  3-3  =  250 

XU  •  / 


T(days) 


a  (km) 


6878. 

6859.584 

6837.117 

6807.003 


PH(perifocal 
height, km) 


492.982 

478.167 

446.441 

409.072 


with  the  previous  case  is  that  the  serai-major  axis  is  decaying  more  rapidly 
due  to  the  increased  atmospheric  density  ass  ~ Lated  with  the  higher  level  of 
solar  activity.  The  character  of  the  decay  is  still  nearly  linear  but  the 
decay  rate  increases  significantly  over  the  thirty-day  span.  Over  the  thirty 
day  period, 

Aa  =  -70.997  km 
APH  =  -83.910  km 

The  average  daily  Aa  =  -2.367  km. 

3.3.4  Stationkeeping  Requirements 

Active  altitude  control  is  necessary  to  maintain  the  initial  orbital  con¬ 
figuration  of  the  MMW  antenna.  Nickerson  et.  al.  provide  linearized  varia¬ 
tional  equations  for  the  maneuvers  about  a  reference  circular  orbit.  These 
equations  assume  that  the  propulsive  impulses  are  along  an  action  line  normal 
to  the  radius  vector  and  in  the  plane  of  the  orbit.  These  equations  are 


Aa  =  2a0  V  T  avt 


Aex=2V~  cos  9  AVt  . 


Ae  =  2~\  / —  sin  0  AV 
y  \  V  T 
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where 


Ae^  =  change  in  e  cos  id 
Ae^  =  change  in  e  sin  <D 
y  =  gravitational  constant 
ao  =  reference  circular  orbit  semi-major  axis 
Avt  =  tangential  velocity  change  impulse 
0  =  argument  of  latitude 

Inspection  of  these  results  shows  that  the  application  of  two  equal  impulses  in 
the  same  direction,  with  the  impulses  being  separated  by  180°  in  0,  results 
in  a  net  change  in  the  semi-major  axis  but  no  change  in  the  magnitude  of  the 
eccentricity. 

For  example,  consider  the  propulsive  requirements  to  restore  the  semi¬ 
major  axis  of  the  orbit  for  the  average  solar  activity  case  (F^q  j  =  150) 
after  thirty  days.  The  data  for  the  maneuvers  equations  is 

Aa  =  8.0155  km 
a  =  6861.969  km 
y  =  398600.8  kmJ/sZ 

Substitution  of  these  quantities  into  the  first  maneuvers  equation  gives 
AV^  =  A. 45  m/s^ 

Thus  the  total  AV  requirement  to  restor  the  orbit  after  thirty  days  is 

8.9  m/s.  Over  a  five-year  lifetime,  the  total  AV  would  be  approximately 
542  m/s. 

The  AV  requirements  are  much  larger,  for  the  high  solar  activity  case. 

The  data  for  the  maneuver  equations  is 

Aa  =  35.50  km 

a  =  6807.003  km 

°  3  2 

y  =  398600.8  km  /s 

These  give 

AVt  =  19.95  m/s 

Thus  the  total  AV  requirement  to  restore  the  orbit  after  thirty  days  is 

39.9  m/s.  Over  a  five-year  lifetime,  the  total  AV  to  be  applied  is  approx¬ 
imately  2428.9  m/s. 

3.3.5  Density  at  the  Orbital  Altitude 

The  maximum  density  at  the  nominal  orbital  altitude  is  of  interest  in 
the  calculation  cf  aerodynamic  torques.  The  values  are: 
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Maximum  density  at  500  km 


F1q  7  =  150  0.2042  E-2  kg/km3 

F10  7  =  250  0.5645  E-2  kg/km3 

3.3.6  Forces  and  Torques  Resultants 

For  the  complex  geometric  shapes  of  a  large  satellite  such  as  the  MMW 
antenna,  the  "exact"  evaluation  of  drag  forces  becomes  extremely  difficult  and 
involved.  First,  the  drag  coefficient  for  each  exposed  component  needs  to  be 
determined.  The  shape  and  size  of  the  element  and  even  the  quality  of  its 
surface  affect  the  coefficient.  In  addition,  it  is  necessary  to  take  into 
account  the  shadowing  effect  and  possibly  the  multiple  collision  of  reflected 
particles.  While  there  exists  specialized  programs  capable  of  accounting  for 
some  of  these  factors,  this  approach  was  not  followed  here.  Because  the  in¬ 
tent  was  only  to  establish  an  order-of-magnitude  for  this  highly  fluctuating 
disturbance,  a  simplified  model  was  deemed  sufficient.  The  spacecraft  y-axis 
was  assumed  to  coincide  with  the  velocity  vector  and  the  satellite  was  repre¬ 
sented  by  only  two  exposed  surfaces.  The  projection  of  the  curved  reflector 

2 

into  the  x-y  plane  determined  the  first  surface  of  599  m  .  The  second  surface 
represented  the  tower.  Since  the  structure  is  not  a  continuous  surface,  an 

2 

equivalent  solid  area  of  43  m  was  calculated  on  the  basis  of  its  solidity 
ratio.  The  relative  positions  of  the  centers  of  pressure  of  the  two  surfaces 
with  respect  to  the  center  of  mass  of  the  spacecraft  were  used  in  determining 
the  torques.  In  the  familiar  aerodynamic  relations  used  for  the  drag  calcu¬ 
lation,  the  drag  coefficient  was  assumed  to  be  2.0.  Sketches  of  the  simpli¬ 
fied  model,  the  mathematical  expressions  used  and  the  calculation  results  are 
presented  in  Figure  3-9.  These  numbers  are  viewed  as  an  estimate  of  the  typi¬ 
cal  disturbance  that  can  result  from  the  atmospheric  drag  encountered  at  a 
500  km  orbit.  The  design  of  the  thruster's  impulses  of  the  attitude  control 
system  was  based  partially  on  these  findings.  (See  Section  3.6.) 

3.4  Solar  Radiation  Pressure 

3.4.1  Introduction 


Solar  radiation  pressure  is  a  major  disturbance  to  satellite  attitudes 
since  the  forces  are  not  balanced  about  the  center  of  gravity.  The  attitude 
control  system  jets  can  provide  the  necessary  corrective  torque  to  prevent  the 
LOS  error  from  becoming  larger  than  the  stipulated  tolerance.  In  order  to  be 
able  to  determine  this  disturbance,  a  complete  time-history  of  the  resultant 
solar  radiation  pressure  torques  and  forces  was  produced  for  the  entire  orbit. 
This  section  first  presents  the  approach  to  solar  radiation  pressure  analysis 
where  equations  involving  specular  and  diffuse  reflection  are  discussed. 

Next,  a  model  for  the  MMW  satellite  is  proposed  and  the  simplifications  neces¬ 
sary  to  bring  the  study  down  to  a  tractable  form  are  presented.  As  a  result, 
plots  of  resultant  forces  and  torques  are  obtained  for  the  entire  orbit  time. 


3.4.2  Solar  Radiation  Pressure  Analysis 


.3 

The  momentum  per-unit-volume  (units  Nsm  )  which  is  carried  by  a  colli- 

2  _2 
mated  beam  of  light  is  given  by  H/c  where  H  is  the  power  per-unit-area  (W  m  ) 

and  c  is  the  velocity  of  light  (ms  ^).  When  the  light  is  totally  absorbed  on 
a  material  body  the  force  exerted  on  the  body  is  the  time-rate  of  deposition 
of  momentum.  Hence,  the  force  per-unit-area  (area  is  normal  to  the  beam) 
exerted  in  the  direction  of  travel  of  the  light  is 

F  =  (momentum/volume)  x  (velocity  of  light) 

F  =  —■  *  c  =  H/c  (units  N  m~ ^) 
c 

For  a  reflecting  surface,  the  rate-of-change  of  momentum  of  the  photons  can 
be  higher  because  the  photons  leaving  the  surface  have  a  component  of  momentum 
in  a  direction  opposite  to  a  component  of  the  arriving  photons.  In  computing 
the  radiation  pressure  on  a  spacecraft  it  is  particularly  significant  to  de¬ 
termine  the  net  resultant  force  per-unit-area  when  the  incident  radiation  is 
diffusely  reflected.  The  term  "diffuse  reflector"  implies  that  the  intensity 
of  radiation  leaving  the  surface  is  independent  of  direction  and  this  is  known 
as  Lambert's  cosine  law.  Because  of  the  azimuth  symmetry  of  the  diffuse  radi¬ 
ation,  there  exists  only  a  net  pressure  force  in  a  direction  opposite  to  that 
of  the  surface  normal  (P^).  This  component  is,  for  normal  incident  light, 

v  H  (1-“> 

o 

where  H  is  the  incident  radiation  flux  intensity  (W/nf)  and  a  is  the  absorp¬ 
tivity  of  the  surface. 

Figure  3-10  summarizes  the  three  types  of  pressure  forces  that  can  be 
modeled  as  part  of  the  orbital  heat  flux  program.  Case  I  represents  a  black 
(totally  absorbing)  surface  and,  as  such,  is  just  a  special  case  of  either 
Case  11  (specular  surface)  or  Case  III  (diffuse  surface).  Figure  3-10  provides 
the  resultant  force  as  resolved  into  its  components  perpendicular  and  parallel 
to  the  surface  (P.^),  per  unit  surface  area. 

Figure  3-11  summarizes  the  radiation  pressure  on  a  flat  black  plate  viewing 
the  earth  as  a  function  of  altitude  for  three  sources  of  radiation,  i.e., 
solar,  reflected  solar  and  earth  IK,  The  location  of  the  flat  plate  over  the 
subsolar  point  of  the  earth  was  chosen  to  maximize  the  pressure  force  due  to 
the  sunlight  reflected  by  the  earth. 

It  is  apparent  from  Figure  3-11,  that  the  solar  pressure  forces  on  a  pay- 
load  can  be  over  an  order-of-magnitude  higher  than  the  pressure  forces  due  to 
earth  IR  and  reflected  sunlight,  particularly  at  the  higher  altitudes. 

3.4.3  MMW  Antenna  Model  and  Results 

The  antenna  structure  was  idealized  as  a  simple  model  of  a  flat  rectangu¬ 
lar  plate  representing  the  slightly  curved  reflector.  The  tower  structure  was 
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3714  n® 


Altitude  (1000  nm) 

Solar:  1353  W/m2 

Albedo:  0.3 

Earth  IK:  237  W/m2 

Figure  3—11 -  Solar,  maximum  albedo  and  earth  IR 
radiation  pressures  on  a  flat  black 
plate  as  a  function  of  altitude. 


not  included  in  the  model,  since  its  intercepting  area  is  negligible  compared 
with  the  reflector  area. 

The  ocher  assumptions  used  in  the  calculation  of  the  sun  radiation  pres¬ 
sure  resultants  are  summarised  here: 

(1)  Sun  radiation  vector  is  in  the  circular  orbit  plane 

(2)  Direct  sun  radiation  only;  albedo  and  1R  effects  are  relatively 
small  and  are  not  included  (See  Figure  3-11) 

(3)  Ail  diffuse  surfaces  are  assumed  to  have  xero  absorptivity  (ct  **  o) 

(A)  No  shading  from  tower  or  due  to  antenna  curvature  is  accounted  for 

The  program  computes  the  total  incident  flux  on  the  surface  as  a  function 
of  orbit  position  and  spacecraft  orientation.  With  the  angle  of  the  incoming 
ray  known  and  the  surface  optical  properties  available,  the  magnitude  and  di¬ 
rection  of  the  resulting  force  are  obtained.  Torques  with  respect  to  the 
satellite  center  of  mass  are  easily  calculated  from  geometrical  relationships. 
Since  the  sun  is  assumed  to  he  in  the  plane  of  the  orbit,  and  owing  to  the 
model  symmetry,  only  two  force  resultants  Fx  and  Fs,  and  one  torque  resultant, 
l'y,  are  present.  Figure  3-12  defines  the  various  parameters  describing  the 
a.itenna-orhit-sun  relationships.  The  two  forces  and  the  one  torque  resultants 
are  plotted  as  a  function  of  time  over  the  entire  95  minute  orbit.  Due  to 
assumption  (2)  above  the  force  resultants  are  identically  xoro  as  the  satellite 
moves  in  the  earth  shadow. 

'  •  5  Gravity-Gradient 

3.5.1  Introduction, 

Thu  effects  of  ‘-.rav it y -grad lent  have  often  been  used  as  a  simple  means  of 
stabilising  small  satellites.  For  large  and  asymmetric  spacecraft  such  as  the 
MMW  antenna,  however,  gravity-gradient  torques  become  significant  and  are  a 
major  source  of  environmental  disturbance.  For  the  MMW  satellite  the  calcu¬ 
lated  torques  are  even  larger  than  those  associated  with  the  effects  of  atmos- 
pho\  It:  drag  at  the  500  km  altitude. 

This  section  discusses  the  treatment  of  the  effects  on  an  orbiting  body  of 
the  earth’s  gravitational  field.  With  proper  selection  of  body  axes,  the  equa¬ 
tions  for  the  torques  can  bo  reduced  to  simple  expressions.  The  torques  for 
the  MMW  spacecraft  are  evaluated  simply  in  terms  of  orbit  geometry  and  mass 
moments-of-inertia.  It  is  interesting  to  note  that  the  gravity-gradient 
torques  for  an  earth  pointing  satellite  are  almost  constant  (the  MMW  antenna 
Is  earth  pointing),  whereas  they  would  vary  constantly  for  an  earth  orbiting 
spacecraft  that  is  Inortlally  fixed. 

3 .5.2  G  rav 1 tv-G  rad  ten  t  Ann 1 y  sis 

While  a  rigorous  treatment  of  the  effects  of  gravity-gradient  can  ho  quite 
elaborate  if  the  gravitational  field  is  expressed  by  an  infinite  series  of 
harmonics,  this  approach  is  less  complicated  as  it  takes  advantage  of  a  number 
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earth  shadow 


Figure  3-12.  Solar  radiation  resultants 


of  legitimate  simplifications.  The  anomalies  due  to  irregularities  of  the 
earth  mass  distribution  and  its  surface  features  can  be  disregarded  in  a  first- 
order  estimation.  The  evaluation  of  torques  accounting  at  least  for  the 
oblateness  of  the  earth  is  usually  recommended,  but  this  approach  considers 
the  simple  case  of  a  spacecraft  in  circular  orbit  about  a  spherical  earth. 

The  derivation  of  the  basic  equations  for  this  "ideal"  case  is  presented  below. 


For  a  body  in  a  circular  orbit  of  radius,  R  ,  about  the  earth,  (see  Fig- 

-3. 

ure  3-13)  a  reference  coordinate  system  (i^,  k^.)  is  defined  such  that  kr 

is  always  directed  radially  toward  the  earth,  and^  its  origin  coincides  with 
the  mass'  center  of  the  body.  The  vector  radius  R  has  components  in  the  refe- 
ence  coordinate  frame  as  follows 


R  =  Rokr 


=  R  i,  +  R  j,  +  R„k 
x  b  yJb  2  b 


(3-1) 


The  components  of  the  vector  R  in  the  body  axes  are  obtained  from  the  rotation 
matrix  relating  the  two  systems  of  axes  by  use  of  the  Euler  angles  ,0  and  (j>. 
They  are: 


R  =  F  sin  0 
x  o 

R  =-R  sin  <f>  cos  0 

y  o 

R  =-R  cos  d>  cos  0 
zo  T 


(3-2) 


If  denotes  the  radius  vector  from  the  body's  center  of  mass  to  a  generic 
mass  element,  dm,  it  has  the  form 

p  =  xib  +  yjfa  +  2k^  (3-3) 

and  the  radius  vector  from  the  earth's  center  to  this  mass  element  becomes 

(3-3) 


^  ^  i 

r  =  R  +  p 


The  force  due  to  gravity  on  a  mass  element  is 

dF-  -if? 

r 

where 


(3-4) 


\l  =  GM 

O 

=  g  R2 
°e  e 

=  398,601  km3/sec2 
G  =  universal  constant  of  gravity 
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1—13 .  Reference  and  body  axes  for  a  body  in  a  Cx.cular  orbit 
about  the  earth. 


This  causes  a  torque  dL  about  the  mass  center  of  the  vehicle 

(Ifxt) 


dL  ='(?  x  dF  =  -  ^ 


Great  simplifications  of  the  equations  are  obtained  because  the  origin  of 
the  body  axes  coincide  with  the  center  of  mass.  Further  simplifications  can 
be  gained  if  these  same  axes  are  the  principal  axes  of  inertia.  With  these 
assumptions  some  algebraic  manipulations  lead  to  the  final  result  for  the 
gravitational  torques  in  scalar  form: 


(3-6) 


L 

3i. 

(I 

-  I  ) 

X 

2R3 

z 

y 

o 

L  = 

y 

3p 

2R3 

(IZ 

-  v 

o 

L  = 

3U 

(I 

- 1 ) 

z 

2R3 

X 

y 

(3-7) 


(3-8) 


Here 


x 


/ 


(y  +  z  )  dm,  etc. 
and  the  product-of-inertia  terms 


xy 


■/ 


xydm 


are  zero,  since  the  principal  axes  of  inertia  are  chosen  to  coincide  with  the 
body  axes.  Note  that  in  the  expression  of  the  torques,  only  tw£  Euler  angles 
(6,  ({))  appear,  because  the  third  angle  </j,  a  rotation  about  the  k  vector  has 
no  effect  on  the  gravitational  force  on  a  mass  element  dm.  r 


3.5.3  MMW  Antenna:  Gravity-Gradient  Torques 


The  MMW  antenna  is  an  earth  pointing  spacecraft  on  an  assumed  circular 
orbit  at  an  altitude  of  500  km.  Because  of  its  asymmetric  configuration  and 
due  to  its  orientation,  the  satellite  principal  axes  of  inertia  do  not  coincide 
with  the  reference  coordinate  system  of  the  orbit.  By  inspection  of  the  equa¬ 
tions  (3-6),  (3-7),  and  (3-8)  derived  above,  it  is  clear  that  gravity-gradient 
torques  will  therefore  be  present.  Figure  3-14  indicates  the  relative  orien¬ 
tation  of  the  two  systems  of  axes  and  the  orbit  parameters  needed  for  the  eval- 
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Figure  3-14a.  Definition  of  satellite  axes  cor  gravity-gradient 
torque  evaluation. 
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•  ORBIT  REFERENCE  AXES 

•  SATELLITE  BODY  AXES 
»  RELATIVE  ORIENTATION 


tr,tr,*kr  UNIT  VECTOR  TRIAD 
VV*b  UN,T  VECTOR  TRIAD 

<t>  =  -4.7° 

0=0 
tp  =  0 


•  I,.  l2,  l3  SATELLITE  PRINCIPAL  MASS  MOMENTS  OF  INERTIA 


t 

1  2R3 

(l3-  l2) 

sin  20  cos3  0 

R  =  Re  +  h 

=  6378  +  500  =  6878  km 

T  -  3K 
t2 

<‘3-‘l> 

sin  20  cos  0 

K  =  398,601  km3/s2(=/i  IN  TEXT) 

2R'* 

ORBITAL  PERIOD  T  =  94.6  min 

T  -  3K 

>3 - T 

2R3 

(1,-12) 

sin  20  sin  0 

STEADY  STATE  TORQUE  T1  =  0.14  n-m 

INCREMENTAL  TORQUES  FOR  PEAK  ATTITUDE  TOLERANCE  OF  2  ror 
AT,  =  0.0034  n-m 

AT2  *  0.0015  n-m 


Figure  3-14i>.  Definition  of  satellite  axes  for  gravity-gradient 
torque  evaluation. 


uation  of  the  torques.  Because  of  the  orientation  of  reference  and  body  axes, 
the  Euler  angle  0  is  equal  to  zero.  As  a  result,  the  only  gravity-gradient 
torque  that  remains  is  about  the  body  x.^  -  axis  (same  as  reference  xf  - 

axis).  This  disturbing  torque  is  called  steady-state  as  opposed  to  incre¬ 
mental  torques  AT^  and  AT2  that  could  be  present  in  case  another  disturbance 

misaligns  the  body  with  respect  to  the  reference  axes.  Incremental  torques 
AT^  and  AT2  have  been  estimated  for  a  peak  attitude  tolerance  of  2  mr.  The 

torque  T^  is  always  zero. 

3.6  Dynamic  On-Board  Disturbance 

3.6.1  Disturbance  Causes 

The  disturbances  on  board  the  MMW  satellite  include  vibrating  equipment 
such  as  a  cryo-cocler,  reaction  forces  from  irregularities  in  the  unbalanced 
pointing  motions  of  communications  antennas  and  reaction  forces  from  attitude 
and  orbit  control  systems.  In  this  case  both  the  attitude  control  and  the 
orbit  control  are  assumed  to  be  done  with  thrusters.  In  order  to  determine 
if  there  can  be  a  serious  problem  due  to  disturbances  it  is  assumed  that 
thrusters  are  not  throttled,  but  are  turned  on  and  off  abruptly.  In  this 
situation  the  thrusters  are  expected  to  produce  much  more  disturbance  than  the 
other  items  mentioned  above,  so  in  the  following  only  the  thrusters  are  con¬ 
sidered. 

3.6.2  Summary  of  Thruster  System  Requirements 

The  thrusters  are  required  to  counteract  the  relatively  steady  drag  force 
and  attitude-disturbing  torques  caused  by  atmospheric  drag,  gravity  gradient, 
and  solar  pressure.  Since  the  only  interest  at  this  time  is  to  roughly  deter¬ 
mine  the  disturbance  which  the  thrusters  cause  in  the  course  of  maintaining 
the  desired  attitude  and  orbital  altitude,  no  specific  control  law  will  be 
derived,  and  the  requirements  on  the  system  will  not  be  specified  in  any  de¬ 
tail.  Table  3-4  summarizes  the  major  external  disturbance  forces  and  torques 
acting  on  the  satellite.  These  have  been  estimated,  as  treated  in  other 
sections.  In  the  presence  of  these  disturbances  the  spacecraft  attitude  must 
be  maintained  to  a  tolerance  of  +  2  mrad. 


Table  3-4.  External  disturbance  forces  and  torques 


Cause 

Drag  Force 

Torque 

Atmospheric  Drag 

(Axis  1  normal  to  orbit  plane) 

Average  Solar  Activity 

(Axis  2  along  velocity  vector) 

Several-Day  Average 

0.0285  N 

Tx  =  0.0081  N-m 

Peak  within  an  Orbit 

0.0760  N 

T1  =  0.0223  N-m 
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Table  3-4  (cont.)  External  disturbance  forces  and  torques 


Cause 

Drag  Force 

Torque 

High  Solar  Activity 

Several-Day  Average 

Peak  within  an  Orbit 

i 

0.105  N 

0.210  N 

T  =  0.030  N-m 

Tx  =  0.0617  N-m 

Solar  Radiation  Pressure 

Average  over  an  Orbit 

Peak  within  an  Orbit 

«  0.007  N 

0.007  N 

Max  T1  =  0.017  N-m 

Gravity-Gradient 

(limits  shown  apply  for 
]e^  or  e2|  <  2mrad) 

T.,  =  0.140  +  1.7e1  = 

0.140  +  0.0034  N-m 

T2  =  5.8e2  =  +  0.0115  N-m 

3.6.3  Control  Rationales 

The  determination  of  th'-nster  force  magnitudes  and  locations  depends  on 
several  satellite  design  decisions.  As  is  shorn  in  Table  3-4  the  gravity- 
gradient  torque  has  a  comparatively  large  steady  component  tending  to  rotate 
the  satellite  in  the  orbit  plane,  due  to  the  fact  that,  with  the  mass  distri¬ 
bution  assumed,  the  satellite’s  inertia  ellipsoid  is  tilted  in  the  orbit  plane 
away  from  the  local  vertical.  It  might  be  practical  to  redistribute  the 
masses  and  essentially  eliminate  the  steady  component  of  gravity-gradient 
torque,  and  in  principle  the  fir.t-order  gradient  torques  could  be  made  zero, 
or  changed  from  destabilizing  to  stabilizing  torques. 

Similarly  it  might  be  possible  to  align  the  center  of  atmospheric  drag 
with  the  center  of  gravity  to  reduce  the  torque  due  to  the  drag  force.  Neither 
of  these  possibilities  has  been  explored  at  this  point,  so  the  results  can  be 
considered  as  conservative  assessments.  In  any  case,  there  is  actually  no 
direct  propellant  penalty  for  not  balancing  the  torques,  since  by  using 
single  unbalanced  thrusters  rather  than  opposing-pair  attitude  thrusters  all 
attitude  control  thrust  for  in-orbit-plane  torques  contributes  to  orbit  al¬ 
titude  maintenance,  which  is  necessary.  (Some  penalty  can  result  from 
requiring  a  lower  thrust  level.) 

With  regard  to  maintenance  of  the  orbit  altitude,  two  options  are  possible, 
and  these  will  be  treated  in  the  following  two  subsections.  One  is  to  quasi- 
continuously  maintain  the  orbit  throughout  both  the  operating  and  non-operating 
parts  of  the  orbit,  using  an  on-off  control  of  selected  thrusters  arranged 
in  an  array.  The  average  thrust  would  be  controlled  to  maintain  the  orbital 
altitude  and  the  low  desired  eccentricity.  The  average  point  of  application 
within  the  perimeter  of  the  thruster  array  would  be  controlled  to  maintain  the 
spacecraft  within  2  mr  of  the  desired  attitude. 
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The  other  option  is  to  do  intermittent  orbit  corrections,  which  would  be 
needed  every  three  to  ten  days  according  to  the  estimated  orbit  decay  rate 
and  altitude  tolerance.  Presumably  these  corrections  could  be  done  in  the 
parts  of  the  o  bit  when  the  system  is  not  operating  so  the  disturbance  to 
the  antenna  structure  would  not  degrade  system  performance.  The  attitude  con¬ 
trol  could  be  done  with  smaller  thrusters,  thus  causing  less  disturbance  than 
in  the  quasi-continuous  option. 

Both  of  these  options  are  treated  in  a  very  simplified  way  in  the  two 
following  subsections,  the  purpose  being  to  roughly  characterize  the  structure 
disturbances  rather  than  design  the  reaction  control  system. 

3.6.4  Quasi-Continuous  Orbit  Maintenance 

As  seen  in  Table  3-5  the  average  thrust  required  for  orbit  maintenance 
varies  from  0.015  N  for  HSA  (high  solar  activity)  to  0.0285  N  for  ASA  (average 
solar  activity)  and  to  much  lower  values  for  low  solar  activity  conditions. 

To  properly  choose  the  thruster  magnitudes  and  locations  such  that,  with  the 
above  average  thrust,  the  average  line  of  application  is  located  so  as  to 
produce  the  torque  needed  to  control  attitude,  requires  appropriate  information 
of  the  variabilities  of  the  various  forces  and  torques  on  the  spacecraft. 

Since  tne  gravity-gradient  torque  (assumed  not  to  be  balanced,  as  noted  pre¬ 
viously)  is  the  largest  torque  component,  variation  of  drag  requires  that 
both  the  magnitude  and  the  point  of  application  of  the  average  thrust  be 
changed.  This  close  coupling  between  attitude  control  and  orbit  control  does 
appear  to  be  quite  undesirable,  but  nonetheless  is  considered  in  assessment 
of  possible  disturbances  to  the  spacecraft. 

The  thrust  calculations  for  the  two  cases  noted  above  are  tabulated  as 
follows,  assuming  all  torques  are  in  the  same  direction. 

Table  3-5.  Thrust  displacement  for  attitude  control 


High  Solar 

Average  Solar 

Activty 

Activity 

Average  thrust 

over  several  orbits 

0.105  N 

0.0285  N 

Required  Torque 

Atmospheric  drag 

0.0617 

0.0223 

Components 

Solar  Radiation 

0.017 

0.017 

Gravity  Gradient 

0.1434 

0.1434 

Total  Torque 

0.222  Nm 

0.183  Nm 

Required  displacement  from  c.g. 

0.222  .  . 
0.105  2'lm 

0.183 

0.0285  b 

This  gives  some  idea  of  vhe  thruster  magnitudes  and  the  distance  over 
which  the  average  thrust  line  must  be  able  to  be  moved  by  varying  the  duty 
cycle  of  the  several  thrusters.  Possibly  a  pair  of  0.2  N  thrusters  located 
about  2m  from  the  c.g.  plus  lower  thrust  units  at  larger  distances  above  and 
below  the  c.g.  would  be  typical  design  results.  Pairs  would  be  spaced  later¬ 
ally  for  yaw  control.  For  the  present  purposes  it  is  enough  to  conclude  that 


the  turning  on  or  off  of  a  thrust  of  0.4  N  is  a  reasonable  estimate  of  the 
disturbance. 


The  next  question  is  how  often  such  a  disturbance  may  occur.  Considering 
a  time  interval  in  which  the  thrusters  are  all  off,  if  we  assume  a  to  and  from 
excursion  AG  =  Imr  is  executed  with  a  torque  T  =  0.2  Nm  (see  Table  3-5)  the 
total  time  would  be 


t 


001)  1.154  x  106 
0.2 


=  215  sec 


If  the  torque  exerted  by  the  thrusters  is  double  the  externally  applied  torque, 
the  thrust  duration  would  be  equal  to  the  off  time  given  above.  Actually  there 
will  be  additional  transients  due  to  the  requirements  to  control  the  other 
spacecraft  axes  (which  have  much  smaller  external  torques),  but  more  impor¬ 
tantly,  it  may  not  be  practical  to  allow  A0  to  be  so  large,  allowing  for 
sensor  errors  and  suitable  performance  margins.  Also,  if  the  thrusters  torque 
is  rather  larger  than  the  externally  applied  torque,  the  on  time  will  be  pro¬ 
portionately  less.  In  conclusion,  it  appears  that  a  disturbance  model  should 
be  the  turn-on  or  turn-off  of  0.4  N  of  thrust  at  about  one  minute  intervals. 


3.6.5  Intermittent  Orbit  Maintenance 


In  this  case  the  orbit  maintenance  burns  are  done  intermittently  at  three 
to  ten  aay  intervals,  during  times  when  the  satellite  is  not  required  to  be 
operating  as  a  sensor.  Therefore  the  disturbance  peak  is  of  no  concern,  and 
only  has  to  decay  sufficiently  before  the  system  must  again  be  operational. 

The  attitude  control  thrusters  can  then  e  placed  at  a  maximum  distance  from 
the  c.g.  so  the  force  levels  can  be  reduced.  Using  the  same  external  torque 
values  as  before,  and  a  thruster  torque  of  double  that,  if  the  thruster  is 
placed  at  the  antenna  feed  (level  arm  of  about  25  meters)  the  force  required  is 


0.4  Nm 
F  =  25^T~ 


0.010  N 


This  is  an  order-of-magnitude  less,  and  the  time  between  transients  could  be 
the  same,  about  one  minute. 


To  get  an  idea  of  the  magnitude  of  thruster  required  for  orbit  maintenance, 
if  the  drag  of  0. 105  N  over  a  three  day  interval  is  compensated  by  two  burns , 
each  occupying  10  degrees  of  orbital  motion  and  spaced  by  180  degrees  around 
the  orbit,  the  thrust  required  is  given  by 

T-  (0.105  N)  (3odays„x  24_x_3600.  =  86  „ 

'360°'  *  ^  m^nutes  x  60 

where  the  orbital  -od  is  95  minutes.  Clearly  there  is  considerable  freedom 
in  the  choice  of  such  parameters. 

Incidentally,  the  propellant  required  for  orbit  maintenance  at  500  km 
altitude  is  appreciable.  The  average  drag  of  0.0285  N  implies  a  loss  of  about 
510  m/sec  velocity  for  a  spacecraft  of  8800  kg  over  a  year  life.  For  a 
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specific  impulse  of  260*  sec  this  means  a  propellant  mass  given  by 


m  =  8800  \e 


510 
260  x  9 


••-J- 


=  1950  kg 


It  would  therefore  be  of  interest  to  consider  raising  the  orbital  altitudes 
somewhat . 

3.6.6  Summary 


The  specific  disturbance  values  given  in  the  preceding  sections  are  pro¬ 
vided  for  use  in  an  exploratory  assessment,  and  are  not  to  be  considered  as 
firm  values.  As  has  been  noted,  the  balancing  of  gravity-gradient  would 
greatly  reduce  the  steady  torque,  and  the  torque  due  to  atmospheric  drag 
could  also  be  largely  balanced  out.  In  addition,  an  increase  in  the  orbital 
altitude  would  save  propellant  mass  and  reduce  the  drag  torque.  If  the  result 
were  to  be  that  thruster  disturbances  were  still  inconveniently  large,  reaction 
wheel  control  could  be  considered  with  momentum  dumping  either  by  small 
thrusters  or  by  active  balancing  of  the  external  torques. 


*Typical  of  shuttle  vernier  RCS  engine,  at  26  lb  thrust. 
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SECTION  4 


STRUCTURAL  MODEL  AND  ANALYSIS 


4.1  Introduction  and  Approach 

The  approach  to  the  mathematical  model  of  the  MMW  structure  was  influenced 
by  both  the  need  for  a  proper  structural  representation  and  the  requirements 
imposed  by  the  evaluation  method  of  antenna  performance.  It  was  the  intent 
here  to  incorporate  into  the  model  all  those  structural  characteristics  that 
could  significantly  affect  the  antenna  behavior.  Another  driving  requirement 
for  this  model  was  the  need  to  represent  the  reflector  surface  by  a  large 
number  of  information  points.  In  terms  of  the  mathematical  model,  such  a 
requirement  rapidly  increases  the  size  of  the  system  matric*  s,  to  a  point 
where  the  solution  becomes  both  susceptible  to  errors  and  abnormally  costly. 
This  conflict  was  easily  circumvented  by  a  substructuring  approach  used  in 
both  the  static  and  the  dynamic  analyses.  Since  the  antenna  structure  con¬ 
sists  of  a  main  assembly  supporting  124  panels,  identical  in  the  geometry  and 
structural  characteristics,  a  substructuring  technique  was  even  more  desirable. 
Cnee  the  accuracy  problems  inherent  to  the  uncoupling  procedures  have  been 
assessed,  the  antenna  analysis  can  be  reduced  to  the  analysis  of  two  basic 
structures:  (a)  a  typical  panel  and  (b)  the  back-up  trusses/receiver  tower 

assembly.  Although  the  combination  of  the  distortions  of  the  individual  com¬ 
ponents  require  additional  computation,  these  calculations  are  relatively 
easy  since  they  are  simple  geometrical  manipulations. 

The  response  calculation  to  the  quasi-static  loading  is  deemed  "exact". 

The  panels  being  supported  on  three  kinematic  mounts,  it  is  legitimate  to 
uncouple  them  structurally  from  the  underlying  back-up  trusses.  In  the  dynam¬ 
ic  analyses,  however,  the  same  approach  leads  to  some  error  in  the  evaluation 
of  the  elastic  distortions.  This  error,  however,  can  be  satisfactorily  re¬ 
duced  if  the  dynamic  characteristics  of  the  panels  are  such  that  they  can 
minimize  the  dynamic  coupling.  To  this  end,  the  design  of  the  panel  was  con¬ 
cerned  with  imparting  adequate  bending  stiffness  to  the  structure  in  order  to 
achieve  a  specified  lowest  frequency  of  approximately  30  Hz.  Since  the  back¬ 
up  structure  will  exhibit  most  of  its  significant  natural  frequencies  below 
the  panel's  fundamental  frequency,  the  strong  interaction  between  these  struc¬ 
tures  is  limited  to  the  higher  frequencies  of  the  trusses.  The  corresponding 
mode  shapes,  although  excited,  show  significantly  smaller  amplitudes  and  are 
expected  to  be  of  negligible  consequence. 

This  section,  therefore,  describes  two  structural  models.  The  back-up 
trusses  and  the  tower  represent  the  first  model  for  which  a  modal  analysis 
uncovers  their  dynamic  characteristics.  The  second  model  is  a  representation 
of  the  typical  hexagonal  reflector  panel.  The  natural  frequencies  and  mode 
shapes  are  presented  for  three  different  panel  constructions,  thus  permitting 
trade-offs  for  the  most  desirable  selection.  The  evaluation  of  the  structural 
responses  is  presented  in  Section  5. 


55 


4 . 2  Modcj.  Descript i on  of  Trusses/Tower  Assembly 
4.2.1  Mathematical  Model 

The  main  structure  of  the  MMW  antenna  includes  a  grid  of  trusses  forming 
a  spherical  surface  and  a  tapered  tower  frame  supporting  the  receiver  array. 

The  mathematical  model  used  is  a  faithful  representation  of  the  actual  struc¬ 
ture,  as  evidenced  by  the  large  size  of  the  stiffness  matrix.  However,  to 
keep  the  number  of  degrees-of-freedom  (DOF)  “o  a  reasonable  level,  the  large 
platform  was  modelled  with  axial  members  whereby  only  three  DOF’s  are  present 
at  each  joint.  Although  the  truss  member  connections  usually  exhibit  good 
bending  capability,  this  simplification  is  legitimate  since  the  axial  load 
path  is  overwhelming  in  the  truss  configurations.  For  the  tower,  however, 
frame  members  were  used  to  insure  stability  of  some  joints.  The  distributed 
mass  from  the  reflector  panels  were  lumped  at  only  34  surface  joints  so  that 
the  dynamic  problem  remained  at  a  tolerable  size.  While  there  are  1131  DOF’s 
in  the  stiffness  matrix  K,  the  total  number  of  entries  along  the  diagonal  mass 
matrix  M  is  limited  to  141.  It  is  useful  to  recall  that  an  accurate  represen¬ 
tation  of  the  elastic  properties  of  a  structure  requires  a  more  refined  mathe¬ 
matical  model  than  the  representation  of  the  inertia  properties  does.  This  is 
due  to  the  fact  that  the  inertia  properties  (or  kinetic  energy)  depend  directly 
on  the  displacements  of  the  structure,  whereas  the  stiffness  properties  (or 
strain  energy)  are  functions  of  the  derivatives  of  the  displacements,  and  it 
is  well  known  that  the  derivatives  have  larger  components  of  error  than  the 
displacements  themselves. 

The  geometrical  description  of  the  mathematical  model  is  best  understood 
by  considering  the  plots  shown  on  Figure  4-1.  The  presence  of  outriggers  at 
the  periphery  and  the  minimum  number  of  diagonal  members  result  in  reducing 
the  system  structural  mass  without  impairing  overall  stability.  The  pattern 
selected  for  the  orthogonal  trusses  (in  projection)  is  such  that  the  inter¬ 
sections  offer  hard-points  to  accept  the  three  kinematic  mounts  for  each  re¬ 
flector  panel.  The  1.5m  deep  trusses  have  a  mass  of  1  kg/m^  including  an 

2 

allowance  of  0.1  kg/m  for  the  thermal  blanket,  while  the  applied  mass  from 

2 

the  hexagonal  graphite-epoxy  panels  amount  to  4  kg/m  . 

Because  of  the  preliminary  nature  of  this  study,  only  two  different  mem¬ 
ber  cross-sections  have  been  used  in  this  model:  (a)  an  open  thin-wall  section 
for  all  the  truss  members  and  (b)  a  thin-wall  tube  for  the  tower  frame  elements. 
The  geometry  of  these  cross-sections,  the  material  used  and  the  various  dynamic 
and  stability  characteristics  of  the  individual  members  are  summarized  on 
Table  4-1. 

The  tower  member  lengths  range  from  5  m  to  15  m  with  most  of  the  members 
at  about  10  m.  In  this  design  the  thin-wall  tube  with  10  cm  diameter  produced 
a  relatively  low  fundamental  frequency  of  11.4  Hz  for  a  1C  m  long  member.  The 
effects  of  the  vibrations  of  single  members  in  the  tower  structure  are  dif¬ 
ficult  to  anticipate  at  this  point,  however,  it  is  always  possible  to  use  a 
larger  cross-section,  or  provide  intermediate  bracing  points  if  this  low  level 
of  natural  frequency  becomes  a  problem.  Table  4-1  compares  the  merits  of  two 
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1  model . 


TRUSS  MEMBERS 

•  MATERIAL  GRAPHITE-EPOXY 


•  FUNDAMENTAL  FREQUENCY 

•  OVERALL  FLEXURAL  BUCKLING 

•  FLEXURAL-TORSIONAL  BUCKLING 


E  =  14  X  1010  N/m2 
w  =  0.32  kg/m' 

A  =  1 .86  cm2 
lx  =  l.j  =  3.74  cm4 
ly  =  12  =  1.38  cm4 

rmin  -  r2  =  0862  cm 
MAXIMUM  LENGTH  L  =  5.16  m 

f  =  10.4  Hz 
Pcr  =  2800  N 
Pcr  =  4400  N 


TOWER  FRAME  MEMBERS 
•  MATERIAL 


•  CROSS  SECTION 
GEOMETRY 


•  MAXIMUM 
MEMBER  LENGTH 

=  15.4  m 

•  CROSS-SECTIONAL 
PROPERTIES 


•  FUNDAMENTAL 
FREQUENCY 


GRAPHITE-EPOXY 


A  =  4.71  cm2 
I  =  58.9  u.i4 

J  =  117.8  cm4 
r  =  3.53  cm 

f  =  4.8  Hz 
PCT=  1360  N 


E  =  14  X  1010N/m2 


A  =  12.57  cm2 
I  =  628  cm4 

J  =  1256  cm4 
r  =  7.07  cm 

f  =  9.6  Hz 
Pcr  =  14500  N 


•  FLEXURAL 
BUCKLING 


tube  cross-sections  for  the  maximum  member  length  of  15.4  (tower  diagonal 

/  2  2W 

diagonal  length  (10  +  11.7  }  =  15.4  m.)  Note  that  the  tower  support  mem¬ 

bers  extending  below  the  reflector  are  even  longer  in  this  model,  but  in  real¬ 
ity,  they  will  be  braced. 


Graphite-epoxy  composite  was  used  for  all  the  structural  members.  The 
excellent  E/p  ratio  (E  =  modulus  of  elasticity,  p  =  density)  coupled  with  a 
very  low  coefficient  of  thermal  expansion  have  made  it  a  prime  candidate  for 
space  applications.  These  desirable  characteristics  are  evidenced  when  com¬ 
pared  to  those  of  other  possible  candidates  are  shown  in  Table  4-2. 

4.2.2  Modal  Analysis 

The  mathematical  model  of  the  overall  MMW  antenna  structure  containing  141 
dynamic  DOF'  .  (47  times  3  translational  masses)  was  analyzed  as  a  free-free 
structure.  he  results  of  the  modal  analysis  show  frequencies  ranging  from 
0.89  Hz  to  x36  Hz  with  60  elastic  modes  having  frequencies  below  25  Hz.  Since 
it  is  realized  that  the  reliability  of  the  results  in  the  higher  modes  is 
questionable  due  to  the  limitation  of  the  modeling,  the  first  60  modes  (non- 
rigid  body  modes  below  25  Hz)  were  used  in  evaluating  both  transient  and  steady- 
state  responses.  This  deletion  is  further  substantiated  by  the  natural  ex¬ 
pectation  that  the  higher  modes  are  less  significant  in  response  calculation. 
Table  4-3  describes  the  nature  of  the  mode  shapes  associated  with  the  first 
10  elastic  modes.  The  plots  of  some  mode  shapes  shown  in  Figure1'  4-2  (a)  to 
(i)  can  be  quite  useful  to  anticipate  the  seriousness  of  performance  degradation 
associated  with  a  particular  mode.  Finally,  it  should  be  noted  that  the 
natural  frequencies  of  the  system  are  uniformly  distributed  and  that,  for  this 
reason,  the  tuning  of  a  mechanical  disturbance  to  avoid  resonance  appears 
nearly  impossible. 


Table  4-2.  Characteristics  of  other  possible  material  candidates 


— 

Material 

E  x  1010 
(n/m2) 

P 

(xlO3) 

•a 

kg/ni 

E/p 

6 

(*10  ) 

m 

V 

cx 

(*10~6) 

m/m/°K 

Aluminum 

7.2 

2.77 

2.7 

0.33 

23.0 

Beryllium 

30.3 

1.83 

16.9 

0.025 

12.0 

Graphite-Epoxy 

13.8 

1.72 

8.2 

0.3 

0.2 

Magnesium 

4.5 

1.83 

2.5 

0.35 

26.0 

Molybdenum 

32.4 

1.02 

3.2 

0.307 

6.6 

Titanium 

11.0 

4.6 

2.4 

0.34 

7.5 

60 


MODE  11 
FREQ  =  2.702  HZ 


Figure  4-2e.  Millimeter  wave  antenna 


The  response  calculation  to  dynamic  disturbances  were  obtained  through 
modal  superposition  in  which  a  uniform  modal  damping  ratio  of  0.005  was  as¬ 
signed  to  each  mode.  These  results  are  presented  in  Section  5.3. 

The  calculation  of  the  structural  responses  to  quasi-static  disturbance 
such  as  the  thermal  loading  was  obtained  by  simple  inversion  of  the  stiffness 
matrix  K.  Before  inverting  the  matrix  however,  it  was  necessary  to  support 
the  structure  in  order  for  K  to  be  positive  definite.  While  the  distortions 
of  the  structural  joints  are  a  function  of  the  choice  of  the  support,  the  ar¬ 
bitrary  selection  of  supports  is  acceptable  since  only  the  relative  displace¬ 
ments  are  of  interest.  Note  that  the  supports  must  keep  the  structure  stati¬ 
cally  determinate  (externally) .  Results  of  quasi-static  loadings  are  presented 
in  Section  5.2. 


Table  4-3.  Natural  frequencies  and  mode  shapes  (MMW) 


Mode 

Frequency 

HZ 

Descrip  tion 

1-6 

0 

Rigid  body  modes 

r— — 

7 

0.9 

~ 

Asymmetric  lateral  surface  beading 

Lateral  tower  bending  1st  surface  torsion  mode 

8 

1.4 

1st  surface  bending  mode 

9 

1.7 

1st  feed  support  bending  mode 

10 

1.8 

Tower  bending  -  surface  twisting 

11 

2.7 

Surface  lateral  bending 

12 

3.4 

Asymmetric  surface  bending 

13 

4.2 

Surface  bending  mode 

14 

4.9 

Surface  twist-tower  bending  and  twisting 

15 

4.9 

Asymmetric  surface  bending 

16 

5.7 

Tower  torsion  mode 

4.3  ?anel  Design  axs.  Analysis 
4.3.1  Introduction 

In  the  attempt  to  design  an  effL  lent  reflector  surface  structure  for  the 
2 

large  area  (about  1600  m  )  of  the  MMW  antenna,  the  designer  is  faced  with  many 


interacting  requirements.  The  selection  of  a  concept  is  affected  by  the  Col 
lowing  driving  parameters. 


(1)  structural  concept:  rigid  segment  or  inflatable 

(2)  reflector  weight  and  segment  size 

(3)  required  surface  accuracy 

(4)  expected  thermal  and  dynamic  loadings 

(5)  active  or  passive  control  of  surface  accuracy 

(6)  exposure  to  adverse  environment 

(7)  packaging  and  deployment  techniques 

Trade-offs  among  these  parameters  should  converge  towards  an  adequate 
compromise. 

This  section  de. is  with  a  segmented  reflector  concept,  and  eo  dots  trade¬ 
off  studies  among  parameters  such  as  weight,  materials,  thickness  uni  funda¬ 
mental  frequency  for  a  hexagonal! y  shaped  panel  of  pre-determined  size.  Three 
panel  designs  are  presented:  (a)  a  graphite-epoxy  egg-crate  structure,  (b)  a 
honeycomb  aluminum  sandwich  construction  and  (c)  an  inflatable  such  as  an  air- 
mat  sandwich  structure  using  Dacron  weave  and  drop  threads.  Since  the  graphite- 
epoxy  panel  appears  to  fulfill  all  the  requirements,  it  was  chosen  to  be  the 
candidate  for  the  MMW  antenna.  The  other  'wo  concepts  have  shortcomings  that 
additional  study  and  design  manipulations  could  possibly  circumvent. 


4.3.2  Hexagonal 

The  choice  of  a  segmented  reflector  concept  appears  appropriate  1  use 
of  the  large  size  of  the  reflector  (approx i'mu-ely  30  m  by  60  m)  and  the  r.-ther 
stringent  requirements  on  Lite  reflector  surface  accuracy  ^figure  error  standard 
deviation  of  0.05  mm)  to  he  maintained  by  a  passive  structure.  The  size  of  the 
segment  is  strongly  dependent  on  the  p;  kaglng  technique  whereby  all  segments 
must  fit  together  in  the  most  efficient  manner  in  the  space  allocated  by  the 
shuttle  cargo  bay  (cylindrical  shape  18.28  m  long  by  4.57  m  in  diameter). 

The  form  or  tlu  shape  of  the  segments  must  be  such  that,  when  assembled,  they 
form  a  uniform  reflector  surface.  Finally,  from  a  manufacturing  point-of-view, 
it  is  desirable  to  fabricate  a  .large  number  of  similar  panels,  thus  Improving 
cost-effectiveness.  All  these  requirements  combined  to  point  toward  the  design 
of  a  hexagonal  "rigid"  panel  with  a  4.5  m  maximum  diametrical  dimension.  i..o 
thickness  had  to  he  Limited  to  6  cm  so  as  to  allow  the  back-to-back  packaging 
of  124  panels  within  a  7.8  m  length  of  the  shuttle  cargo.  Geometrically,  the 
hexagonal  panels  when  deployed  and  assembled  can  form  the  desired  surface  of  an 
el! ipt ‘eai-torus,  and  their  spatial  recurrent  ■  facilitates  the  choice  of  a 
pattern  for  the  design  of  the  underlying  trusses.  Because  of  these  desirable 
characteristics,  the  hexagonal  shape  of  the  panel  remained  an  invariant  in  the 
following  trade-off  studies.  The  ultimate  choice  was  based  on  the  inerLts  of  a 
particular  material  and  structural  concept  responsible  for  the  pi.nei  thermal 
and  dynamic  characteristics. 


4.3.3  Mathematical  Model  and  Analysis  of  the  Hexagonal  Panel 

The  hexagonal  panel  was  analyzed  as  a  plate  represented  by  44  finite  ele¬ 
ments  (30  rectangular,  14  triangular  elements).  The  selection  of  the  struc¬ 
tural  joints  was  such  that  they  coincide  with  optical  grid  points  used  in  the 
antenna  performance  evaluation,  thus  circumventing  costly  interpolations.  The 
model  geometry  and  its  characteristics  are  given  in  Figure  4-3.  Although  the 
underlying  truss  configuration  offers  support  points  for  the  "corners"  of  the 
hexagons,  this  choice  of  support  is  not  optimum  for  the  panel.  In  this  study, 
the  three  kinematic  mounts  are  located  on  radial  lines,  0.9  m  inward  from  the 
corners.  This  location  was  chosen  because  it  is  ve’ y  close  to  the  optimum  sup¬ 
port  condition  in  terms  of  the  high  fundamental  frequency.  Three  translations,, 
masses  at  each  unrestrained  DOF  account  for  the  147  dynamic  DOF's. 

The  mechanical  properties  of  an  equivalent  plate  have  been  calculated  for 
the  three  design  candidates.  The  same  finite-element  model  can  be  used  for 
the  three  panel  construction  once  the  bending  and  shear  stiffness  are  known 
and  the  proper  material  properties  ar^  available  Table  4-4  summarizes  the 
features  of  the  3  panel  candidates  anc.  their  respective  plate  stiffness.  While 
the  graphite-epoxy  egg-crate  and  the  aluminum  honeycomb  construction  concepts 
are  well  understood,  the  inflatable  airmat  is  somewhat  unconventional.  A 
brief  description  of  this  structural  concept  and  its  dominant  characteristics 
is  therefore,  given  below. 

Packageability  and  automatic  erectability  make  inflatable  construction 
especially  attractive  for  satellites  and  space  stations.  An  example  of  in¬ 
flated  construction  is  an  inflated  fabric  plate  called  Airmat.  Airmat,  which 
was  developed  by  Goodyear  Aircraft  Corporation,  consists  of  two  covers  of  coated, 
woven  material  connected  by  closely  spaced  flexible  fibers  called  drop  cords 
(or  drop  threads).  The  airmat  takes  on  its  predetermined  conf igurat.' on  as  a 
result  of  the  internal  pressure.  A  theory  for  the  structural  analysis  of  in¬ 
flatable  fabric  plates  has  been  developed  to  provide  information  on  their 
load-carrying  capabilities  and  structural  characteristics.  Static  tests  as 
well  as  vibration  tests  have  been  performed  on  inflatable  fabric  plates  in 
order  to  assess  the  validity  of  the  theory.  The  design  of  the  hexagonal  panel 
using  the  airmat  construction  is  based  on  typical  values  for  the  choice  of 
airmat  geometry,  woven  material,  drop  cords  and  pressure  levels.  As  evidenced 
by  comparing  the  plat  :  stiffness  in  Table  4-4,  the  12  cm  deep  airmat  exhibits 
an  adequate  bending  stiffness,  but  has  a  very  low  carrying  capability  in  shear. 
Since  the  shear  stiffness  is  proportional  to  the  internal  pressure  p,  an  in¬ 
crease  in  pressure  seems  reasonable.  However,  such  a  move  is  limited  by  rhe 
amplitude  of  the  surface  ripples  that  could  be  caused  by  the  deflection  of  the 
face  sheet  between  the  point  supports  provided  by  the  drop-threads.  A  reason¬ 
able  compromise  is  therefore  needed.  Table  4-4  lists  the  various  airmat  par¬ 
ameters  used  in  the  proposed  panel  design. 

The  results  of  the  finite-element  model  analysis  are  summarized  in  Table 
4-5.  A  modal  analysis  run  for  each  candidate  revealed  their  dynamic  charac¬ 
teristics  with  a  list  of  the  first  eight  natural  frequencies.  An  estimate  of 
their  thermal  behaviour  was  obtained  on  the  basis  of  the  thermal  model  developed 
in  bection  3.2  for  the  graphite-epoxy  panel.  Through  extrapolations,  the 
maximum  thermal  gradient  AT^  in  the  aluminum  honeycomb  and  the  airmat  have 
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Figure.  4-3.  Mathematical  model  of  the  reflector  surface  panel. 
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been  established  for  an  identical  thermal  flux.  The  corresponding  deflections 
at  the  center  of  the  plate  are  listed  in  the  same  Table. 

As  mentioned  above,  the  egg-crate  panel  made  of  graphite-epoxy  appears  to 
satisfy  both  dynamic  and  thermal  requirements  and  consequently  was  chosen  as 
the  candidate.  The  aluminum  honeycomb  design  shows  a  lower  fundamental  fre¬ 
quency,  principally  because  of  its  lower  modulus  of  elasticity.  It  is  inter¬ 
esting  to  note  that  the  thermal  gradient  is  almost  3  times  smaller  due  to  the 
aluminum  high  conductivity  that  tends  to  even  out  the  face  sheet  temperatures. 
This,  however,  is  not  enough  to  cancel  the  effects  of  the  aluminum  large  co¬ 
efficient  of  thermal  expansion  a.  A.n  increase  in  the  panel  thickness  could  im¬ 
prove  both  its  dynamic  and  thermal  characteristics,  but  it  would  affect  ad¬ 
versely  the  packaging  ratio.  The  greatest  advantage  of  the  airmat  concept  lies 
in  its  high  packageability  and  relatively  low  mass.  Based  on  perhaps  old  in¬ 
formation  (Goodyear  Corp.  studies  date  back  to  the  mid-sixties)  the  design  of 
the  airmat  does  not  show  great  promise.  Its  ability  to  minimize  the  thermal 
gradient  AT^  through  its  thickness  is  greatly  impaired  by  the  presence  of  only 
drop  threads,  thus  creating  relatively  large  thermal  deformation.  These  con¬ 
clusions,  based  on  somewhat  old  data  should  not  be  considered  definitive.  New 
'studies  of  the  possibilities  offered  by  inflatables,  in  general,  and  airmats, 
in  particular,  should  be  undertaken  for  this  kind  of  structural  applications. 

If  the  shortcomings  in  their  thermal  behavior  can  be  overcome,  the  inflatables 
may  well  prove  superior  to  other  concepts,  especially  because  of  their  initial 
"flexibility"  from  a  packaging  point  of  view. 

The  hexagonal  panel  model  approximately  represents  all  124  reflector 
panels,  and  will  be  used  in  estimating  the  absolute  thermal  distortions  of  the 
surface.  In  order  to  show  the  performance  degradation  as  a  function  of  a 
particular  design,  the  three  models  will  be  used.  For  the  dynamic  response, 
the  panels  are  considered  to  vibrate  outside  the  range  of  interest  and  there¬ 
fore,  are  assumed  to  be  "rigid".  Their  deformations  of  rigid-body  nature  are 
imparted  by  the  relative  motions  of  their  3  kinematic  mounts  (from  underlying 
truss  displacements).  The  results  of  the  response  calculations  are  presented 
in  Section  5. 


76 


SECTION  5 


RESPONSE  EVALUATION 


5.1  Introduction  and  Approach 

Section  4  described  the  mathematical  modeling  of  both  the  trusses/tower 
assembly  and  the  hexagonal  reflector  panels.  The  structural  responses  have 
been  evaluated  as  the  geometric  combination  of  the  elastic  distortions  of  the 
component  substructures.  For  the  static  case,  the  uncoupling  of  the  reflector 
panels  from  the  underlying  trusses  was  possible  because  of  the  kinematic 
connection.  Each  panel  was  first  analysed  independently.  The  rigid-body  dis¬ 
placement  components  associated  with  deformations  of  the  supporting  back-up 
trusses  were  added  to  the  elastic  distortions.  With  124  panels,  each  having 
51  surface  points,  the  deformed  reflector  surface  was  finally  expressed  by 
three  displacement  components  of  6324  points. 

Ir.  the  response  calculation  to  dynamic  disturbance,  the  124  reflector 
panels  were  assumed  to  be  "rigid".  (The  rationale  for  this  assumption  was 
discussed  in  Section  4.)  As  the  back-up  trusses  vibrate,  they  move  the  panels 
thus  imposing  rigid-body  displacements  and  rotations  on  them.  The  elastic 
distortions  of  the  panels  are  disregarded  because  it  is  assumed  that  the  ampli¬ 
tudes  of  panel  vibrations  are  negligible  compared  to  the  amplitudes  of  truss 
vibrations.  The  absolute  displacements  of  all  6324  surface  points  were  ob¬ 
tained  by  evaluating  the  rigid-body  motions  of  each  one  of  the  124  panels. 
Figure  5-1  illustrates  the  steps  involved  in  the  computation  of  the  surface 
displacements.  The  antenna  performance  evaluation  based  on  the  reflector  de¬ 
formed  surface  is  presented  in  Section  6. 

Section  3  established  from  a  list  of  disturbances,  those  that  are  expected 
to  be  predominant  in  distorting  the  reflector  surface  and  in  affecting  the 
antenna  pointing  accuracy.  The  responses  calculation  therefore  was  limited 
to  four  loadings  capable  of  illustrating  the  antenna  behavior  under  normal 
environmental  conditions.  The  only  quasi-static  loading  is  the  thermal 
disturbance  that  corresponds  to  a  68°  sun  flux  incidence  with  respect  to  the 
z-axis  in  the  X-Z  plane  of  the  antenna  model.  Responses  to  this  loading  have 
been  evaluated  for  the  three  different  panel  constructions  presented  and  dis¬ 
cussed  in  Section  4.3.  Three  dynamic  disturbances  have  been  considered:  the 
first  two  are  step-impulse  functions  associated  with  the  firing  of  thrusters 
for  two  possible  attitude  control  schemes  and  the  third  one  is  a  harmonic  dis¬ 
turbance  simulating  a  cryo-cooler  reciprocating  machine.  A  detailed  descrip¬ 
tion  of  these  loading  and  the  associated  structural  responses  are  presented 
below. 

5.2  Response  to  Sun  Thermal  Flux  (68°) 


FOR  PANEL  k.  k  -  1,  2...124: 


(1)  ELASTIC  DISTORTIONS  ON  3  KINE¬ 
MATIC  MOUNT  SUPPORTS 

6f'  (i  -  1.  50) 

NOTE:  Sf1  s  0  IN  DYNAMIC  ANALYSIS 

(2)  PANEL  RIGID  BODY  DISPLACEMENTS  (A) 
AND  ROTATIONS  (</>)  FROM  BACK-UP 
TRUSSES  DISPLACEMENTS 

THEN  OBTAIN 

6,RB  (i  -  1,  50) 


Figure  5-1.  Reflector  surface  distortion  calculation. 


critical,  because  it  is  the  only  orbital  position  where  the  sun  impinges  di¬ 
rectly  onto  the  reflector  surface.  This  condition  creates  the  greatest  thermal 
gradient  across  the  reflector  panels.  Since  the  orbit  inclination  is  not  de¬ 
finitively  known,  the  sun  vector  orientation  in  the  X-Z  plane  was  chosen  to 
represent  a  worst  case  condition.  In  addition,  due  to  the  curvature  of  the 
reflector  along  its  longer  dimension  (about  60  m)  the  sun  incidence  with 
respect  to  the  various  panels  shows  the  widest  variation  of  about  30°  from  one 
edge  to  the  other.  In  order  to  account  for  this  variation,  the  reflector 
surface  was  arbitrarily  divided  into  four  equal  width  zones  of  approximately 
15  m.  The  thermal  conditions  of  the  structures  were  assumed  to  be  identical 
throughout  each  zone.  Although  the  panels  of  a  zone  theoretically  are  not 
under  exactly  the  same  thermal  condition,  this  restriction  was  found  acceptable 
and  greatly  decreased  the  amount  of  computations.  In  a  similar  manner,  the 
underlying  truss  assembly  was  assigned  four  temperature  distributions  for 
the  same  four  15  m  wide  zones.  Figure  5-2  illustrates  schematically  the  geom- 
try  of  the  zones  considered.  Since  the  tower  members  were  not  included  in 
this  thermal  model,  information  on  tower  thermal  distortions  are  only  approx¬ 
imations  in  this  analysis.  A  uniform  temperature  of  220°K  was  arbitrarily 
used  for  the  tower  members  of  the  model,  but  a  separate  analysis  for  the  pre¬ 
diction  of  tower  distortions  was  carried  out  based  on  a  special  thermal  model 
of  the  tower  and  will  be  presented  later  in  this  section. 

The  thermal  response  calculations  included  a  static  analysis  of  the  hex¬ 
agonal  panel  for  the  four  different  thermal  distributions  associated  with  the 
four  zones.  In  a  parallel  analysis,  the  truss-tower  assembly  was  subjected 
to  four  temperature  zones  as  obtained  by  the  thermal  model  described  earlier. 
The  temperatures  used  in  the  analysis  appear  in  Table  5-1.  The  distortions 
have  been  computed  for  an  arbitrary  reference  temperature  of  280°K.  A  special 
purpose  computer  program  was  developed  in  order  to  integrate  into  the  overall 
structure  the  124  panel  distortions  in  their  proper  locations.  Each  panel 
was  assigned  three  truss  joints  for  its  three  support  points.  Once  the  truss 
joints  displacements  are  known,  the  panel  support  motions  are  readily  avail¬ 
able.  On  the  basis  of  these  displacements,  the  rigid-body  motions  of  the 
panel  are  determined.  As  a  consequence,  each  panel  surface  point  (there  are 
51  points  per  panel)  undergoes  additional  displacements  as  a  function  of  the 
panel  three  rigid-body  translations,  three  rigid-body  rotations  and  location 
within  that  panel.  The  panel  distortions  finally  are  computed  as  the  sum  of 
elastic  and  rigid-body  displacements  and  then  expressed  in  the  global  coor¬ 
dinate  system  of  the  entire  reflector.  This  information  was  stored  on  tapes 
for  subsequent  use  in  the  performance  evaluation.  Note  that  the  computation 
described  above  was  repeated  two  more  times,  to  evaluate  the  surface  distor¬ 
tions  for  the  aluminum  honeycomb  panel  and  the  airmat  construction.  Figure 
5-3  depicts  the  out-of-plane  deformations  of  three  longitudinal  lines  on  the 
reflector  surface.  A  maximum  deflection  of  0.8  mm  is  recorded  at  the  reflector 
edge. 


The  response  of  the  receiver  tower  to  the  sun  thermal  flux  was  the  sub¬ 
ject  of  a  special  analysis  as  explained  in  Section  3.2.  In  this  case  the 
sun  vector  was  allowed  to  sweep  the  Y-Z  plane  as  shown  in  Figure  5-4.  The 
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TEMPERATURE  DISTRIBUTION 

o  Back-Up  Trusses:  4  Zones 
o  Panels:  4  Zones 

a)  Graphite-epoxy  a  =  0,5  x  10_6/°K 

o)  Aluminum  Honeycomb  a  =  23  x  10~6/°k 

o  Tower:  Uniform  Temperature  Assumed  at  220°K 


Figure  3-2.  Thermal  loading  for  sun  at  68°  with  z-axis 


Table  5-1.  Trusses  and  Panel  Temperatures  (°K) 
(Sun  Incidence  Angle  of  68°  with 
Respect  to  Z  Axis) 


Zone 

Trusses 

Panels 

Top 

Bottom 

Graphite 

i-Epoxy 

Aluminum 

Top 

Bottom 

Top 

Bottom 

1 

287.3 

286.8 

288.1 

285 

286.06 

285 

2 

285.8 

285.2 

287.0 

285 

285.68 

285 

3 

284.5 

283.9 

286.3 

285 

285.44 

285 

4 

283.4 

282.7 

285.8 

285 

285.28 

285 

PANEL  TOP 


PANEL 

BOTTOM 


1:7*1  I 


LINE  A 


WMMMMMmmwmrMmrm 

mmmammmmmmmmmm 


LINE  B 


Hi 

vmmmmmlMmmmmm 


J 


shadowing  effects  of  the  opaque  large  reflector  were  accounted  for  as  well 
as  the  shadow  that  one  tower  number  might  have  on  its  neighbor.  Since  the 
thermal  model  was  reduced  to  a  planar  problem,  only  the  displacement  Ay  and 
Az  were  calculated.  They  are  expressed  in  terms  of  the  angle  4>  for  the  en¬ 
tire  360°  sweep.  One  should  note  that  with  respect  to  a  mean  value  the  tower 
top  displacement  Ay  would  vary  by  +0.4  mm  while  Az  would  remain  within  ±0.2 
mm.  Whereas  the  Az  displacement  is  acceptable  when  compared  to  the  1  mm  tol¬ 
erance,  the  dec.enter  Ay  exceeds  by  far  the  stipulated  limit  of  0.15  mm.  In 
order  to  alleviate  the  deleterious  effects  of  the  sun  on  the  tower  decenter, 
several  possibilities  can  be  suggested.  As  an  example  of  passive  measures, 
an  improvement  of  the  coefficient  of  thermal  expansion  by  a  factor  of  3,  or 
possibly  a  judicious  combination  of  a  and  e  for  the  various  surface  quality 
could  reduce  Ay  sufficiently.  As  a  last  resort,  it  is  conceivable  to  design 
a  semi-active  control  system  that  would  displace  the  receiver  array  according 
to  a  predetermined  law  along  the  orbital  path.  While  none  of  these  sugges¬ 
tions  have  been  studied  in  detail,  it  can  be  concluded  that  the  implementation 
of  one  method  or  a  combination  of  them  should  keep  the  receiver  within  the 
tolerable  travel.  Figure  5-4  illustrates  the  tower  tip  displacement  histories 
as  a  function  of  sun  incidence  angle  d>. 

5.3  Approximate  Distortion  Calculations 

The  outcome  of  the  thermal  response  analysis  of  the  mathematical  model 
can  be  predicted  somewhat  based  on  the  distortion  evaluation  of  the  following 
ideal  case.  The  problem  of  the  distortion  of  a  homogeneous  spherical  reflec¬ 
tor  exposed  to  a  thermal  flux  can  serve  to  illustrate  the  relative  effects  of 
uniform  temperature  rise  and  thermal  gradient.  Consider  a  spherical  reflector 
of  curvature  1/R  and  having  a  circular  shape  of  radius  a.  The  plate  (thick¬ 
ness  h)  is  homogeneous  and  the  material  exhibits  a  uniform  coefficient  of  ther¬ 
mal  expansion  a.  When  exposed  to  the  sun  thermal  flux,  the  reflector  under¬ 
goes  deformation  that  can  be  essentially  attributed  to  a  uniform  temperature 
rise  AT^  and  a  thermal  gradient  AT^  across  its  plate  thickness.  Figure  5-5 

shows  the  reflector  geometry  and  the  expressions  for  the  AT  and  AT  as  a 

U  b 

function  of  plate  geometrical  and  mechanical  properties. 

Assuming  a  free-free  structure  or  equivalently,  a  plate  simply  supported 
along  the  periphery,  the  maximum  sag  6  occurring  at  the  center  is  expressed 
by  two  terms  as 


2  /AT  AT  \ 

a  a  | _ u  _ G  j 

2  \  R  h  / 


As  expected,  a  uniform  rise  in  temperature  causes  a  deflection  only  if 

there  is  a  curvature  1/R.  In  the  case  of  a  flat  plate,  R  =  and  no  distortion 

takes  place.  The  gradient  AT  is  usually  a  more  significant  source  of  distor- 

b 

tions.  For  a  given  AT^,  the  sag  5  is  inversely  proportional  to  the  plate 
thickness  h. 
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Figure  5-4.  Receiver  support  tower  thermal  distortions  (sun 
in  z-y  plane) . 


ABSOLUTE  TEMPERATURE 


Applying  the  expression  for  «$  in  the  case  of  an  approximate  circular 
plate  model  of  the  typical  hexagonal  panel,  the  following  data,  based  on 
thermal  model,  are  used: 


a 

= 

4.5  m  (0.5)  =  2.25  m 

R 

= 

60  m 

h 

= 

0.06  ir. 

a 

= 

0.5  x  10  6  m/m°K  (conservative) 

AT 

u 

= 

O 

O 

atg 

= 

3.1° 

6 

= 

h(2.2S)2  *  0.5  x  10  6[;“  +  o'oej 

= 

0.8  pm  +  65.4  pm 

- 

66.2  pm 

— 

0.066  mm 

The  numbers  show  that  the  contribution  of  the  gradient  AT  is  overwhelm- 

o 

ing  in  the  computation  of  the  sag.  A  large  uniform  rise,  (or  fall),  such  as 
the  one  associated  with  the  orbit  injection,  that  is  200°  would  only  cause 
a  sag  of  4.2  pm. 

5.4  Dynamic  Responses 

5.4.1  Introduction 


The  dynamic  responses  of  the  antenna  surface  were  obtained  through  the 
vibrational  distortions  of  the  underlying  trusses  only,  since  the  surface 
panels  have  been  assumed  to  undergo  rigid-body  motions  only.  This  assumption 
was  made  possible  by  comparing  the  dynamic  character-' sties  of  the  truss/tower 
assembly  and  the  typical  surface  panels.  The  natural  frequencies  of  the 
lower  modes  of  the  truss  structure,  that  are  predominant  in  the  distortion 
calculations,  fall  well  below  the  fundamental  frequency  of  the  reflector  panel. 
It  is  therefore  anticipated,  that  while  the  trusses  vibrate  at  these  low 
frequencies,  the  panel  would  indeed  behave  as  a  rigid  body.  At  higher  fre¬ 
quencies,  the  panel  is  expected  to  undergo  small  vibrational  distortions,  but 
its  rigid-body  motions  become  negligibly  small.  One  should  emphasize  that 
these  qualitative  statements  would  require  substantiation  if  the  calculated 
distortions  are  of  the  order  of  the  allowable  ones.  In  this  study,  it  appears 
that  the  dynamic  disturbances  cause  the  reflector  to  distort  well  within 
allowable  limits. 
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It  was  also  assumed  that  the  panels  could  be  replaced  by  only  transla¬ 
tional  masses.  While  the  trusses  and  the  panels  can  be  uncoupled  elastically 
(as  was  properly  done  in  the  static  analysis)  they  remain  coupled  inertially. 
Although  it  can  be  expected  that  the  panels  rotational  masses  would  not  affect 
the  responses  significantly,  additional  study  would  be  necessary  to  substan¬ 
tiate  this  fact  or  establish  its  importance. 

The  response  calculations  have  been  performed  for  three  dynamic  disturb¬ 
ances.  As  discussed  in  Section  3.6,  the  attitude  control  of  the  spacecraft 
can  be  achieved  by  two  different  thruster  control  laws.  One  possibility  con¬ 
sists  of  frequent,  but  small  bursts,  while  another  requires  two  substantial 
thrusts  approximately  180°  apart  executed  once  every  several  days.  Figure 
5-6  shows  the  location  of  the  thrusters  and  the  amplitudes  of  the  step-functions 
used  in  the  evaluation  of  the  responses.  The  third  dynamic  disturbance  con¬ 
siders  a  harmonic  forcing  function  that  could  represent  the  reciprocating 
machine  of  a  cryo-cooler.  One  of  the  possible  designs  for  the  receiver  array 
requires  the  temperature  to  be  kept  at  40°K,  and  this  can  only  be  achieved 
through  the  help  of  an  active  cooling  system.  The  location  for  the  disturb¬ 
ance  was  deliberately  chosen  to  produce  significant  and  unsymmetrical  distor¬ 
tions  (see  Figure  5-b).  The  characteristics  of  a  cryo-cooler  with  an  ampli¬ 
tude  of  15  N  at  a  frequency  of  10  Hz  corresponds  to  a  typical,  state-of-the- 
art  machine. 

The  methods  used  in  evaluating  the  responses  are  briefly  explained  and • 
results  presented  for  selected  locations  show  that  the  dynamic  distortions 
are  relatively  small  and  less  serious  than  the  thermal  ones. 

5.4.2  Dynamic  Response  to  Thrusters  Impulses 

The  firing  of  the  attitude  control  system  thrusters  was  believed  to  be 
one  of  the  critical  dynamic  disturbances  affecting  the.  MMW  antenna  satellite. 

The  characteristics  of  the  forcing  function  can  be  represented  mathematically 
by  a  step  impulse  with  amplitudes  and  duration  based  on  the  attitude  control 
philosopny  presented  in  Section  3.6.  Since  the  principle  of  superposition  is 
valid  for  this  time- invariant  linear  mechanical  system,  a  response  to  a  par¬ 
ticular  step  impulse  can  be  scaled  up  or  down  for  a  similar  loading  of  differ¬ 
ent  amplitude.  In  other  words,  the  surface  distortion  computed  for  the  86  N 
impulses  can  be  scaled  down  proportionally  for  the  other  impulse  loading  of 
0.4  N. 


Owing  to  the  decoupling  property  of  the  modal  analysis,  each  mode  was 
considered  as  a  single  DOF  system.  It  is  instructive  to  examine  the  response 
characteristics  of  such  a  system  to  a  step  impulse  in  terms  of  natural  fre¬ 
quency  and  damping  ratio.  Figure  5-7  depicts  the  response  function  and  the 
expression  for  the  maximum  overshoot  of  the  response,  its  peak  time  and 
settling  time  (5%  remaining  oscillation  about  the  static  deflection)  are  pre¬ 
sented.  Since  the  peak  time  of  the  response  is  inversely  proportional  to  the 
natural  frequency  u>,  the  peak  response  of  the  lowest  mode  (with  the  smallest 
u))  occurs  after  all  othar  modes  have  peaked.  For  this  reason,  the  response 
including  the  first  53  modes,  was  calculated  for  1.2  seconds,  a  time  corres¬ 
ponding  to  approximately  twice  the  time  to  the  peak  response  of  the  fundamental 
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•  ATTITUDE  CONTROL  THRUSTERS  SIMULATED 
BY  STEP  IMPULSE  FUNCTIONS  F(t)  *  FQ  U$(t) 


A)  QUASI-CONTINUOUS  ORBIT  COMPENSATION 

F0  “  MO  N 

B)  INTERMITTENT  ORBIT  COMPENSATION 

W /  2  BURNS.  180°  APART 
Fq  =  86  N 

•  POSSIBLE  CRYO-COOLER  FOR  RECEIVER  REFRIGERA¬ 
TION,  SIMULATED  BY  A  HARMONIC  DISTURBANCE 
F(t)  *  FQ  SIN  0  t 

WITH  Fq  *  15  N^O  =  10  Hz 


Figure.  5-6.  Description  of  dynamic  loadings. 


mode.  A  numerical  integration  using  the  Runge-Kutta  technique  with  a 
time-step  size  At  =  0.044  s  calculated  each  modal  response  in  which  a  critical 
damping  ratio  of  0.005  (%  percent)  was  assumed.  The  superposition  of  the 
various  modal  responses  was  routinely  performed  to  provide  the  physical 
truss  displacements  as  a  function  of  time. 


Figure  5-7.  Response  characteristics  to  a  unit  step  impulse. 


In  the  evaluation  of  the  structural  response  to  the  86  N  thruster  step 
impulse,  the  integration  time  was  selected  such  that  the  responses  of  all  53 
modes  under  consideration  could  reach  their  peak.  Because  the  fundamental 
frequency  is  0.89  Hz,  the  time  of  maximum  overshoot  for  this  mode  is 

t  =  -  =  0.56  s 

p  u/l-t;* 

and  all  other  modes  reach  their  peak  modal  response  in  the  interval  0  <  t  < 
0.56  s.  Since  it  is  possible  that  a  combination  of  modal  responses  can  result 
for  a  particular  DOF  in  a  larger  displacement  some  time  after  0.56  s,  the 
integration  was  carried  out  to  approximately  twice  as  long,  i.e.,  a  time 
t  =  1.2  s.  It  is  reasonable  to  anticipate  that  in  this  interval,  the  ampli¬ 
tudes  of  the  displacements  include  their  maximum  from  which  they  will  de¬ 
crease  steadily  due  to  the  energy  absorption  by  damping.  It  is  interesting 
to  follow  the  z-displacement  of  structural  joints  along  a  longitudinal  cross- 
section  of  the  antenna.  As  expected,  several  joints  reach  their  maximum  z- 
displacement  in  the  neighborhood  of  0.6  to  0.8  s  which  indicates  that  the 
lower  modes  are  predominant  for  the  out-of-plane  displacement  of  the  re¬ 
flector.  Figure  5-8  indicates  the  location  of  the  joints  considered  and 
Figure  5-9  shows  the  corresponding  response  histories.  At  the  reflector’s 
edge  the  maximum  z-distortion  reaches  0.1  mm,  while  distortions  inside  the 

-4 

reflector  area  remain  well  below  10  m. 

The  joints  located  at  the  "corners"  of  the  back-up  trusses  assembly  are 
expected  to  correspond  to  the  location  of  maximum  dynamic  response.  Figure 
5-10  shows  the  response  histories  of  these  joints  where  the  maximum  Az  is 
approximately  0.16  mm  for  3  joints  (261,  50,  and  270)  and  0.22  mm  for  joint 
41.  These  displacements  are  used  later  in  the  performance  evaluation  pre¬ 
sented  in  Section  6. 

The  dynamic  response  of  the  tower  tip  to  the  86  N  thrusters  impulse  needs 
special  consideration.  Figure  5-11  illustrates  the  responses  for  the  Ax,  Ay, 
and  Az  displacement  that  can  be  interpreted  as  2  decenters  and  a  defocus.  At 
the  start  of  the  response,  the  maximum  decenters  of  0.22  mm  (x-direction)  and 
0.32  mm  (y-direction)  are  slightly  larger  than  the  tolerance  (0.15  mm),  while 
the  defocus  of  0.06  mm  falls  below  the  tolerable  1  mm.  Since  these  displace¬ 
ments  decrease  with  time,  they  are  not  considered  critical  to  the  antenna 
performance.  Thus,  the  receiver  array  motions,  as  a  result  of  an  86  N 
thruster's  impulse  remain  acceptable. 

It  is  concluded  that  even  the  largest  thrust  level  does  not  cause  unac¬ 
ceptable  distortion,  so  the  highest  level  orbit  maintenance  thrust  need  not 
be  restricted  to  nonoperational  parts  of  the  orbit,  and  the  system  can  operate 
during  the  orbit  maintenance  thrust  periods.  The  lower  thrust  levels  for 
quasi-continuous  orbit  maintenance,  and  the  even  lower  levels  for  only 
attitude  control  likewise  do  not  cause  significant  antenr.a  deformation. 
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Figure  5-8.  Truss/tower  model.  Joints  identification  numbers. 


z-displacement  (m)  histories  for  the  first  1.2 
of  response  to  86  N  thruster's  impulse. 


z-displacement  (ra)  histories  for  the  first 
response  to  86  N  thruster's  impulse. 


Figure  5-10b.  z-displacement  (m)  histories  for  the 
response  to  86  N  thruster's  impulse. 


. . . 


5.4.3  Steady-State  Response  to  Harmonic  Disturbance 

The  case  of  the  steady-state  response  to  a  harmonic  forcing-function 
is  of  particular  interest.  When  the  exciting  frequency  coincides  with  a 
natural-mode  frequency  of  the  system,  the  amplitude  of  the  excited  mode  tends 
to  grow  infinitely  large.  However,  damping  places  an  upper  limit  that  de¬ 
fines  the  resonance  psak.  While  it  is  reasonable  to  anticipate  that  the 
largest  deflections  still  occur  at  the  time  of  the  largest  amplitude  response 
of  the  resonant  mode  when  many  modes  are  present,  it  is  of  interest  to  compute 
the  displacements  histories  along  a  full  cycle  of  the  steady-state  response. 
This  approach  was  followed  in  the  response  evaluation  of  the  truss/tower 
assembly  to  the  harmonic  disturbance  of  15  N  at  approximately  10  Hz.  To  be 
exact,  the  forcing  frequency  was  selected  to  coincide  with  the  system's 
natural  frequency  of  9.945  Hz  (flexible  mode  number  22),  the  closest  one  to 
10  Hz. 

In  order  to  establish  the  steady-state  response  of  the  system's  physical 
coordinates,  the  responses  of  the  individual  modes  were  first  evaluated  as 
single  DOF  systems.  It  is  therefore  desirable  to  examine  the  response  charac¬ 
teristics  of  a  simple  oscillator  to  a  sinusoidal  forcing  function.  Figure  5-12 
describes  the  system  variables.'  The  expression  for  the  total  response  includ¬ 
ing  transient  and  steady-state  components  is  discussed  below. 

The  basic  equation  of  motion  of  the  single  DOF  system 


mu  +  cu  +  ku  =  p  sin  Jit 
o 


(5-1) 


can  be  rewritten  in  the  more  convenient  form 


where 


•  2  *o 

u  +  2ruu  +  w  u  =  —  sin  Pit 
m 


(5-2) 


=  =  natural  frequency  of  the  undamped  system 

=  exciting  frequency 


=  2mw  =  critical  damping  value 

Q 

-  — —  =  damping  ratio 

r 


-f 


fbt  ■  p0*flt 


Figure  5-12.  Simple  oscillator  parameter  definition. 


The  general  solution  including  the  complementary  and  particular  solution  is 


u(t)  =  S.  (A  sin  (D^t  +  B  cos  u^t)- 


sin  fit  -  It,  —  cos  fit 


(5-3) 


where  o)p  =  m/l-C  is  the  damped  natural  frequency,  and  A  and  B  could  be  eval¬ 
uated  for  any  given  initial  conditions.  The  first  term  of  the  equation  repre¬ 
sents  the  transient  response  to  the  loading.  Since  it  damps  out  quickly 

according  to  the  exponential  £  it  is  of  little  interest.  The  steady-state 
part  of  the  response,  at  the  frequency  fi  of  the  applied  loading  but  out  of 
phase  can  be  rewritten  in  the  following  simple  form 


where 

is 


/  \  o 

u(t)  =  — j 

|  mco 

steady 

state 

the  phase  angle  0,  by  which  the  response  lags  behind  the  applied  load. 


r  ia\2  2  iu\2  r1/2 

|  -  (!)  j 


sin(fit  -  0) 


(5-4) 


0  - 


(5-5) 


for  the  special  case  when  fi  =  w,  the  response  reaches  a  maximum  amplitude 
called  the  resonance  peak,  u 


upeak 


(5-6) 


Each  one  of  the  53  flexible  modes  considered  exhibits  the  response  shown 
by  equation  (5-4).  While  all  modes  vibrate  with  the  same  frequency  fi,  they 
are  out  of  phase  with  the  exciting  function  by  the  phase  angle  0  (equation 
5-5).  It  is  a  simple  task  to  take  an  arbitrary  time  span  of  2ir/fi  correspond¬ 
ing  to  a  full  cycle  and  sum  up  the  contribution  of  all  the  modes,  including 
the  effects  of  phase  lags.  Thus,  complete  information  was  obtained  for  the 
antenna  displacements  during  the  steady-state  response.  The  numerical  values 
used  were: 


9.7 


15  N 


9.945  Hz 

62.49  rad/s 

0.005  for  all  modes 

natural  frequencies  (r  =  1,2,..., 53) 

and  the  location  of  the  disturbance  is  shown  in  Figure  5-6.  In  using  equation 
(5-4)  for  the  modal  response,  the  amplitude  of  the  forcing  function  pQ  must 

be  multiplied  by  the  corresponding  modal  coordinate  and  the  mass  m  is  the 

2 

generalized  mass  of  that  mode.  The  term  mw^  can  be  understood  to  be  the  spring 

constant  k  of  the  mode  r.  Although  the  results  included  respon  histories 
of  three  displacement  components  of  all  the  truss/tower  assembly  joints,  it 
was  deemed  unnecessary  to  present  them  all  Since  they  remained  negligibly 
small  in  terms  of  antenna  performance.  Since  a  tolerance  on  rms  surface  devia¬ 
tion  was  tentatively  set  at  150  pm,  the  maximum  surface  point  deviations  of 
about  +5  pm  are  not  anticipated  to  cause  any  serious  performance  degradation. 

As  an  illustration,  the  z-dlsplacement  of  the  reflector  centerline  (in 
the  long  60  m  direction)  was  plotted  for  three  time  snapshots  during  the 
steady-state  cycle  (see  Figure  5-13).  Also,  the  z-displacements  of  the  ex¬ 
treme  points  of  the  reflector  are  shown  against  time  for  the  full  0.1  s 
period  of  the  steady-state  responses.  As  expected,  the  responses  follow  a 
sinusoidal  curve  with  the  same  frequency  0  of  the  excitation,  (see  Figure 
5-14). 

The  distortions  of  the  tower  tip  need  special  attention  since  they 
correspond  to  the  location  of  the  receiver  array.  The  three  components  of 
displacements  Ax,  Ay,  and  Az  of  the  tower  tip  center  joint  801  are  plotted 
for  the  same  0.1  s  cycle.  The  maximum  distortion  for  any  such  component 
hardly  gjes  beyond  2  pm  and  thus  falls  well  within  the  stipulated  feed  toler¬ 
ance  (allowable  Az  =  1000  pm.  Ay  =  Ax  -  150  pm)  on  defocus  and  decenter  (see 
Figure  5-15). 

In  summary,  the  antenna  surface  distortions  and  feed  motions  resulting 
from  the  harmonic  disturbance  characterizing  a  cryo-cooler  remain  well  below 
critical  levels.  As  a  consequence,  the  antenna  performance  should  remain 
satisfactory  for  this  loading. 
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SECTION  6 


DEGRADATION  OF  ANTENNA  PERFORMANCE 
DUE  TO  SOLAR  FLUX 


6. 1  Introduction 

This  section  describes  the  degradation  of  point-spread  functions  (PSF) 
of  a  30  x  60  m  elliptic  torus  antenna  due  to  its  deformations  produced  by 
solar  flux.  The  antenna  is  assumed  to  be  operating  at  a  wavelength  of  3  mm 
with  a  cosine  feed  taper,  and  the  materials  considered  are  graphite-epoxy 
and  aluminum.  It  is  shown  that  the  graphite-epoxy  antenna  performs  much 
better  than  the  aluminum  antenna.  Strehl  ratios,  line-of-sight  error  (LOS), 
power  on  detector,  and  profiles  of  the  PSFs  are  given  for  three  regions  of 
the  antenna.  Airmat  antennas  were  also  considered,  but  owing  to  their  large 
deformations  (~X)  under  solar  flux,  they  were  not  investigated  in  any  detail. 

6.2  Antenna  Description 

Smoke,  dust,  fog,  and  clouds  frequently  prevent  satellite-based  observa¬ 
tion  of  earth  targets  with  visual  and  infrared  sensors.  The  ability  to  pene¬ 
trate  these  cover  conditions  and  to  perform  equally  well  by  night  or  day  makes 
microwave  passive  sensing  an  attractive  candidate  to  significantly  extend  the 
capabilities  of  satellite-based  observations  beyond  their  present  limitations. 

Passive  microwave  sensors  respond  to  millimeter  wave  thermal  radiation 
emitted  by  all  objects  in  proportion  to  their  temperature  and  surface  emissi- 
vity.  Useful  sensing  of  military  targets  results  predominantly  form  high 
contrasts  of  target  emissivity  with  respect  to  the  background  surface.  Oxygen 
and  water  vapor  are  the  only  atmospheric  gases  with  significant  absorptions. 
These  combine  to  establish  the  clear  air  transmission  windows  centered  at 
frequencies  35  GHz,  95  GHz,  140  GHz,  and  230  GHz. 

The  reflector  antenna  under  Investigation  is  an  elliptic  torus  having 
dimensions  of  30  x  60  m.  The  feed  horns  are  deployed  on  a  cylindrical  sur¬ 
face  passing  through  the  focal  circle  of  the  ellipses.  The  other  focal 
circle  lies  on  the  surface  of  the  earth.  The  design  of  the  feed  horns  pro¬ 
vides  compensation  for  the  aberration  caused  by  circular  cross  track  re¬ 
flector  profile.  Each  horn  accepts  incident  energy  from  approximately  half 
of  the  reflecting  surface  and  provides  for  sidelobe  control  through  proper 
tapering  of  the  energy  accepted. 

As  indicated  in  Figure  6-1,  the  antenna  is  made  up  of  124  hexagonal 
segme^t-R,  the  side  of  a  hexagon  being  2.25  m.  The  pupil  is  assumed  to  be 
circular  and  30  m  in  diameter.  Its  apodization  due  to  tapering  of  the  ac¬ 
cepted  energy  is  assumed  to  be  a  cosine  distribution  such  that  it  is  unity 
.at  the  center  and  zero  at  the  edge  of  the  pupil.  The  taper  distribution  out¬ 
side  the  pupil  was  not  available,  and  therefore  it  was  assumed  to  be  zero. 

The  antenna  is  assumed  to  be  operating  at  a  wavelength  of  3  mm  corresponding 
to  a  frequency  of  95  MHz.  The  system  parameters  are  tabulated  in  Table  6-1. 
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Table  6-1.  System  Parameters 


Antenna  Shape 
Antenna  Size 
Focal  Length 
Pupil  Diameter 
Taper 

Wavelength 


Elliptic  Torus 
30  x  60  m 
30  m 
30  m 
Cosine 
3  mm 


Figure  6-1.  Elliptic  torus  30  x  60  m  segmented 
antenna.  The  total  number  of 
segments  is  124. 


103 


6.3  Point-Spread  Function  Evaluation 


The  evaluation  of  the  PSF  on  an  imaging  system  in  the  presence  of  thermal 
and  dynamic  disturbances  is  described  in  "Computer  Simulation  of  a  Large 
Adaptive  Optical  System",  by  V.  Mahajan  and  J.  Govignon,  reprinted  from  the 
proceedings  of  the  Society  of  Photo-Optical  Instrumentation  Engineers,  January 
29-February  1,  1979,  Tuscon,  Arizona.  In  this  section,  the  evaluation  of  the 
PSF  and  associated  quantities  is  considered  for  a  circular  pupil  with  cosine 
apodization. 

Assuming  an  aberration-free  antenna  in  the  absence  of  thermal  loading, 
the  point-spread  function  (PSF),  i.e.,  the  irradiancc  distribution  of  the 
image  of  a  point  target,  in  the  presence  of  thermal  loading  is  given  by 

I(x,y)  =  -ppr  |//  P(u,v)  exp  [-  (xu  +  yv)]  dudv(2  (6-1) 

where  X  is  the  wavelength  of  the  radiation  and  f  is  the  focal  length  of  the 
antenna.  The  effective  pupil  function  P(u,v)  at  a  point  (u,v)  on  the  pupil 
is  given  by 


P(u,v)  =  Aq  cos  [ 


ti(u~+v~) 


r2ni 


]  exp  [~jp  6W(u,v) } ,  inside  the  pupil 


=  0,  outside  the  pupil 


(6-2) 


where  D  =  30  m  is  the  diameter  of  the  pupil,  Aq  is  the  irradiance  at  the  pupil 
from  the  target,  the  cosine  factor  represents  the  feed  taper,  and  6W  is  the 
wave  aberration  function  of  the  antenna.  The  irradiance  distribution  I(x,y) 
is  the  effective  PSF  of  the  antenna  system. 

To  evaluate  the  diffraction  integral  in  Eq.  (6-1),  by  a  discrete  Fourier 
transform,  we  let 


u 

=  jAu 

V 

=  kAv 

X 

=  £Ax 

y 

=  mAy 

(6-3) 

where  j,  k,  £,  and  m  are  integers,  and  (Au,  Av)  and  (Ax,  Ay)  are  the  spatial 
periods  in  the  pupil  and  image  planes,  respectively.  Therefore,  Eq.  (6—1) 
can  be  written 


I(£Ax,mAy)  = 


MuAvV 
l  Xf  j 


A2o  |  T. 

jk 


V  expl~ 


2TTi 

Xf 


(j£AxAu  +  kmAyAv)] 


(6-4) 


where 


p  =  ■—  P(jAu,kAv) 
jk  A0 


(6-5) 


In  our  fast  Fourier  transform  (FFT)  program,  the  PSF  is  given  by 


ltm  '  Z2  I1  VI’['W(ftf!)]| 

N  jk  u  v  1 


(6-6) 


where  N  *  N  is  the  number  of  points  in  the  FFT  grid  and  N  is  the  number  of 
u  v  55 

points  Inside  the  pupil.  Note  that  Eq.  (6-6)  is  normalized  in  such  a  way 

that  the  irradiance  at  the  center  of  the  image  loo  =  1  when  there  are  no 

aberrations  and  the  magnitude  of  the  pupil  function  is  constant  and  equal  to 

unity.  Comparing  Eqs.  (6-4)  and  (6-6)  we  find  that 


I(£Ax,mAy)  =  Ao  N2  I{ 


(6-7) 


Ax  =  Xf/N  Au 
u 


(6-8) 


Ay  =  Xf/N^Av 


(6-9) 


For  uniform  pupil  irradiance  A{j  from  the  target,  the  total  power  in  its 
image  is  given  by 

D/2 

P  =  2itAo  /  cos2(rr/D)r  dr 
0 


=  0.2335  D2a2 


(6-iO) 


The  fraction  of  power  incident  on  a  detector  placed  in  the  image  plane  is 
given  by 


P,  =  —  fS  I(x,y)  dx  dy 
t  det 


(6-11) 


where  the  integration  is  carried  over  the  detector  area.  If  the  detector  is 
centered  on  the  aberration-free  image  center  and  its  size  is  2jAx  x  2jAy, 
then  the  fraction  of  power  incident  on  it  may  be  written 


P  =  ^  E  E  w  I(£Ax,mAy) 
3  t  £=-j  m=-j 


(6-12) 


jjl  A. 


where  w  are  the  weights  of  the  integration  rule  used.  Substituting  Eqs. 
(6-7)  anS  (6-10)  into  Eq.  (6-12),  we  obtain 


P. 

3 


4.28 


N2AuAv 

N  N  B1 
u  v 


3 

E 

Z=-3 


j 

I 

m=-j 


w  1 
£m  £m 


(6-13) 


6.4  Antenna  Aberrations  Due  to  Solar  Flux 


As  indicated  in  Figure  6-2,  the  thermal  loading  from  the  sun  is  at  a 
maximum  when  the  antenna,  in  a  500  km  orbit  from  the  surface  of  the  earth, 
is  at  an  angle  of  ±68°  from  the  line  joining  the  sun  and  the  center  of  the 
earth.  The  antenna  is  heated  nonuniformly  and,  therefore,  deforms  causing 

degradation  of  its  performance.  Assume  that  a  certain  point  on  the  antenna 

surface  moves  to  r^  because  of  the  thermal  loading.  Let  the  displacement  of 

the  point  be  written 


Sr.  =  r!  -  r.  (6-14) 

ill 

For  small  6r  ,  the  corrocponding  wave  aberration  contributed  to  a  ray  incident 

y 

at  the  point  i\  is  given  by* 

*  *4- 

6W.  =  -2n.  •  6r.  cos  8.  (6-15) 

i  i  i  x 

where  n^  is  a  unic  vector  along  the  normal  to  the  antenna  surface  at  the  point 

r^,  and  8^  is  the  angle  of  incidence  of  the  ray  on  the  undeformed  surface. 

The  wave  aberration  is  positive  if  the  path  length  is  increased  and  negative  if 
ic  is  decreased.  For  simplicity,  we  assume  that  cos  0i  =  1  so  that 


6W±  «  -2n±  ♦  6r^  (6-16) 

We  also  assume  that  a  point  (u,v)  on  the  pupil  is  identical  to  a  corresponding 
point  on  the  antenna  surface. 

In  the  evaluation  of  the  antenna  performance,  each  hexagonal  segment  was 
represented  as  shown  in  Figure  6-3,  by  a  periodic,  array  of  30  points.  The 
spatial  periods  along  the  u  and  v  axes  were  Au  =  0.675  m  and  Av  =  0.650  m, 
respectively.  The  slight  difference  in  the  two  spatial  periods  has  to  be 
accommodated  in  order  to  achieve  continuity  of  the  array  from  segment  to 
segment.  The  corresponding  distribution  of  points  in  a  30  n  diameter  circular 
pupil  is  shown  in  Figure  6-4. 


V,  Mahajan  and  J.  Govignon,  "Computer  Simulation  of  a  Large  Optical  System", 
the  Proceedings  of  the  Society  of  Photo-Optical  Instrumentation  Engineers, 
January  29-February  1,  1979,  Tuscon,  Arizona,  p.  440,  Ec.  4. 


106 


SATELLITE  POSITION 
CONSIDERED 


Ui 

<2 

1| 

o  2 

£< 

Ui  > 

a.  lu 


* 

CD 

CD 

O 

6 

t— 

X 

X 

in 

ro 

X 

o' 

CM 

o 

ii 

It 

O 

CM 

a 

8 

<N 

<§* 

Q 

LU 

S 

3 

ui 

t- 

< 

CO 

CO 

< 

0 c 

ui 

o 

CO 

a: 

O 

a 

Ui 

> 

5 

O 

u 

> 

Ui 

3 

H 

< 

cc 

LU 

X 

o 

a. 

z 

o 

a. 

5 

ui 

ui 

X 

H 

CO 

i 

UI 

s 

2 

LU 

H 

D 

cc 

2 

O 

fsj 

I 

a. 

< 

H 

5 

D 

O 

LL 

2 

a: 

_l 

3 

. , 

a 

< 

CO 

-J 

a: 

Ui 

2 

< 

5 

ui 

5 

< 

o 

CL 

\— 

• 

• 

107 


Figure  6-2.  Thermal  loading  due  to  sun. 
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Figure  6-4.  Distribution  of  pupil  points.  The  total 
number  of  points  inside  a  circle  of  30  m 
diameter  is  1575.  The  spacings  between 
two  adjacent  points  along  the  u  and  v  axes 
are  Au  =  0.675  m  and  Av  =  0.650  m. 
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As  indicated  in  Figure  6-2,  the  antenna  was  divided  into  4  zones.  Only 
one  segment  within  a  given  zone  was  analyzed  for  its  faceplate  warping  due  to 
tl  =  solar  flux,  and  all  the  other  segments  within  that  zone  were  assumed  to 
ut.iave  similarly.  However,  the  rigid-body  movements  of  a  segment  due  to  the 
trusses  were  calculated  for  each  segment.  The  combined  effect  of  faceplate 
warping  and  truss  bending  gave  the  net  deformati,.i  of  a  segment. 

The  numbers  of  points  along  the  u  and  v  axes  are  chosen  to  be  odd  so  as 
to  define  a  center  of  the  pupil  and  they  are  equal  to  47  and  45,  respectively. 
The  total  number  of  points  in  the  pupil  are  1575,  i.e.,  N  =  1575.  The  deforma¬ 
tions  due  to  the  solar  flux  were  determined  at  antenna  points  lying  in  an  array 
of  47  x  97  points.  As  indicated  in  Figure  6-5,  the  center  of  the  pupil  was 
chosen  at  three  different  positions: 


(1)  (0,  19) 

(2)  (0,0) 

(3)  (0,-19) 

The  peak-to-peak  and  standard  deviation  of  aberration  of  the  whole  antenna 
and  portions  of  it  corresponding  to  pupils  centered  at  the  three  positions 
for  both  graphite-epoxy  and  aluminum  antennas  are  given  in  Table  6-2. 

Table  6-2.  Antenna  Aberrations  in  Millimeters  due 
to  Solar  Flux 


Standard  Deviation  of 
Aberration 


Peak-to-Peak  Aberration 


Graphite-Epoxy  0.94  0.34  0.32  0.70  0,15  0.07  0.05  0.16 

Aluminum  1.71  0.54  0.76  1.45  0.30  0.11  0.15  0.28 


Pupi 

1  Position 

1 

2 

3 

0.07 

0.05 

0.16 

0.11 

0.15 

0.28 

6.5  Results 


The  PSFs  have  been  determined  by  using  an  FTT  grid  for  which  N  =  N  =  256. 

u  v 

The  various  evaluation  parameters  are  tabulated  in  Table  6-3.  The  effect  of 
the  solar  flux  on  the  PSFs  is  shown  in  Figures  6-6  and  6-7  for  the  graphite- 
epoxy  and  aluminum  antennas,  respectively.  In  these  figures,  the  ar.tenna 
performance  is  indicated  by  two  numbers:  the  Strehl  ratio  S  and  the  LOS  error. 
The  Strehl  ratio  represents  the  central  irradiance  of  the  PSF  relative  to  its 
value  in  the  absence  of  the  solar  flux.  The  quantity  100(1  -  S)  is  generally 
referred  to  as  the  percent  loss  in  antenna  gain.  The  LOS  error  was  defined 
as  the  movement  of  the  location  of  the  peak  irradiance  of  the  PSF.  Note  that 
the  Strehl  ratio  is  degraded  significantly  in  the  case  of  aluminum  antenna 
compared  to  the  graphite-epoxy  antenna.  Moreover,  position  3,  which  corre¬ 
sponds  to  maximum  solar  flux,  gives  the  smallest  Strehl  ratio  in  both  cases. 
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Table  6-3.  PSF  Evaluation  Parameters 


No.  of  points 
in  pupil 

N  = 

1575 

Pupil  sampling 
interval 

>  t> 

<  c 

ii  n 

0.675  m 

0.650  m 

PSF  sampling 
interval 

Ax  = 
Ay  = 

0.52  mm,  Aa  = 
0.54  mm,  AS  = 

17.4 

18.0 

prad 

prad 

FFT  grid  size 

N 

u 

256,  N  =  256 

V 

From  Eqs.  (6-8)  and  (6-9)  we  note  that  the  spatial  periods  of  the  PSF  data 
points  are  given  by 


and 


Ax 


Xf 

N  Au 
u 


®  0.52  mm 


Ay 


Xf 

N  Av 

V 


=  0.54  mm 


The  corresponding  spatial  periods  in  angular  units  are  given  by 


(6-17) 


(6-18) 


and 


Act 


Ax 

f 


=  17.4  prad 


(6-19) 


Ay 


A z 

f 


=  18.0  prad 


(6-20) 


Since  the  FFT  grid  size  is  256  x  256  the  total  image  area  covered  by  the  PSF 
data  is  133  *  138  mm,  or  4.43  x  4.61  mrad. 


The  LOS  error,  which  was  determined  in  units  of  (Aa,  AS),  can  be  deter¬ 
mined  more  accurately  by  determining  the  center  of  the  reference  sphere  with 
respect  to  which  the  variance  of  the  aberrations  is  minimum.  The  centroid 
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of  the  PSF  can  also  be  defined  as  the  LOS.  However,  in  this  report  the  LOS 
was  defined  simply  as  the  location  of  the  peak  irradiance  in  the  PSF  data  as 
determined  by  the  FFT  program.  It  is  evident  from  Figures  6-6  and  6-7  that 
the  LOS  error  is  much  larger  than  the  tolerance  value  of  5  prad. 


Details  of  the  PSF  for  a  circular  pupil  with  cosine  taper  are  shown  in 
Figure  6-8,  where  a  central  profile  along  the  x-axis  is  plotted.  The  PSF  is 
normalized  so  that  its  value  at  the  center  is  unity.  The  linear  and  angular 
units  of  the  horizontal  scale  are  0.52  mm  and  17.4  yrad,  as  given  b-  Eqs. 
(6-16)  and  (6-18),  respectively.  Figure  6-8  shows  approximately  20  sidelobes. 
Because  of  the  solar  flux  on  the  antenna,  the  irradiance  of  the  sidelobes 
increases.  This  is  shown  in  Figures  6-9a  through  6-9f  for  the  three  pupil 
positions  on  the  graphite-epoxy  antenna.  Central  profiles  along  the  x  and  > 
axes  are  shown.  The  linear  and  angular  units  of  the  horizontal  scale  for 
the  y-profiles  are  0.54  mm  and  18.0  prad,  respectively.  The  principal  side¬ 
lobes  in  these  profiles,  two  in  the  x  and  one  in  the  y  profiles,  occur  because 
of  the  grating  effects  produced  by  the  segmented  nature  of  the  antenna  and  the 
manner  in  which  the  deformations  in  a  given  zone  of  the  antenna  were  compiled. 
The  irradiance  value  noted  in  Figures  6--9a  through  6-9f  is  the  maximum  value 
of  irradiance  in  the  profile  relative  to  a  central  value  of  unity  for  the 
aberration-free  case.  Similar  results  for  the  aluminum  antenna  are  shown  in 
Figures  6-10a  through  6-10f.  These  figures  show  how  the  sidelobe  structure 
changes.  The  poor  performance  of  the  aluminum  antenna  is  evident  from  Figure 
6-10f.  Numerical  analyses  of  these  results  were  not  carried  out.  For  example 
the  irradiance  of  the  far-off  sidelobes  shorn  in  all  of  the  figures  cannot  be 
too  accurate  because  of  aliasing  in  the  FFT  operation. 


The  fraction  of  power  P.  incident  on  a  detector  in  the  image  plane  is 

shown  as  a  function  of  detector  size  j  in  Figure  6-11.  Again,  note  the  signif¬ 
icant  degradation  of  detector  power  (a  measure  of  antenna  performance)  in  the 
case  of  the  aluminum  antenna  compared  to  the  graphite-epoxy  antenna. 


6.6  Discussion  and  Conclusions 


The  performance  of  graphite-epoxy  and  aluminum  antennas  operating  at  a 
wavelength  of  3  mm  has  been  evaluated  when  they  are  irradiated  by  solar 
radiation.  The  degradation  of  the  PSFs  due  to  antenna  warping  caused  by  solar 
flux  is  given  in  terms  of  Strehl  ratio,  LOS  error,  power  on  an  image-plane 
detector,  and  sidelobe  structure.  It  is  shown  that  the  LOS  error  is  much 
above  the  tolerance  value,  and  that  a  graphite-epoxy  antenna  performs  much 
better  then  an  aluminum  antenna.  Airmat  antennas  were  also  considered,  but 
owing  to  their  large  deformations  (~X)  under  solar  flux,  they  were  not  investi¬ 
gated  in  any  detail. 
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Figure  6-9c.  PSF  x-profile  for  pupil  position 
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Figure  6-10a.  Sidelobe  structure  for  aluminum  antenna  under  solar 
flux.  PSF  x-profil.e  for  pupil  position  1. 


SECTION  7 


PACKAGING  AND  DEPLOYMENT 


7.1  Introduction  and  Issues 

As  the. size  of,  the  reflector  increases  and  surface  accuracy  requirements 
remain  relatively  stringent,  the  launching  and  orbit  injection  of  an  entire 
satellite  poses  a  very  challenging  problem.  The  severe  size  and  volume  con¬ 
straints  imposed  by  the  cargo  bay  of  the  Space  Shuttle,  in  combination  with 
the  required  high  deployment  reliability,  increase  the  complexity  of  the 
folding  and  deployment  technique.  Efficient  folding,  packaging  and  stowing 
of  the  space  structure  into  a  relatively  small  volume  become  important  factors 
in  the  satellite  design  trade-offs,  since  they  impact  all  the  structural 
components  and  may  affect  their  respective  performance.  (See  surface  panel 
design  trade-offs  in  Section  4.3). 

The  relative  merits  of  folding  and  deployment  concepts  can  be  evaluated 
as  a  function  of  the  following  parameters 

(1)  overall  packaging  ratio 

(2)  design  size  and  concept 

-  inflatable 

-  variable  geometry 

-  segmented 

(3)  symmetry  and  components  repetitiveness 

(4)  folding  mechanical  complexity 

(5)  reliability 

The  segmentation  approach  must  be  considered  in  close  connection  with  the 
problems  involved  with  the  contemplated  mode  of  in-orbit  assembly,  including 
the  choice  of  parking  or  final  orbit.  Automated  deployment  and  erection,  re¬ 
motely  controlled,  computer  controlled,  or  with  man-in-the-loop,  can  be  al¬ 
ternatives  resulting  from  trade-offs  between  cost,  design  complexity  and 
reliability. 

The  launch  vehicle  dynamic  characteristics  also  strongly  affect  the  fold¬ 
ing  technique  which  needs  to  provide  the  stowed  package  with  protection  from 
the  effect  of  launch  acceleration  and  vibration.  Temporary  supports  and  ties 
can  be  provided  for  launch  in  order  to  ensure  the  integrity  of  satellite 
components. 

The  deployment  ’  locedure  relies  on  mechanical  devices  such  as  hinges, 
spring  flexure  joir.s,  locking  mechanisms  that  will  undoubtedly  increase 
the  design  compI<", Jty.  For  a  precision  structure,  trade-offs  between  perform¬ 
ance,  reliabili'y  and  structural  deployment  complexity  must  be  evaluated.  It 


appears  possible  to  enhance  simplicity  by  enforcing  certain  symmetry  charac¬ 
teristics  of  the  folding  scheme.  A  reflector,  if  divided  and  packaged  in  a 
symmetrical  manner,  provides  a  desirable  repeatability  for  the  algorithm 
of  surface  accuracy  control  of  its  identical  panels.  The  preservation  of 
symmetry,  even  if  it  creates  some  additional  mechanical  complexity  in  deploy¬ 
ment,  is  strongly  recommended.  Manufacturing  complexity  and  cost  are  also 
reduced  by  the  repetitive  fabrication  of  similar  reflector  panels. 

It  is  evident  that  the  pac?.aging  and  deployment  schemes  affect  the 
structural  design  philosophy.  1  desirable  features  for  an  efficient  de¬ 
sign  of  a  large  deployable  struct  *e  and  the  issues  associated  with  the 
various  phases  from  design  to  deployment  are  summarized  in  Table  7-1.  It  is 
apparent  that  they  are  strongly  interactive,  often  equivalent,  and  sometimes 
complimentary.  The  overriding  need  is  to  maximize  the  stiffness  of  the 
structure  by  distributing  the  allocated  mass  in  the  most  efficient  manner. 

But  the  design  concept  must  remain  amenable  to  an  acceptable  packaging  tech¬ 
nique.  This  exemplified  the  need  for  an  integrated  approach  that  ultimately 
satisfies  the  simultaneous  requirements  of  performance,  packaging,  and  de¬ 
ployment  . 

7.2  Packaging  Concept  for  the  MMW  Antenna  Structure 

One  of  the  main  objectives  of  this  study  was  to  demonstrate  that  a  large 
MMW  antenna  performing  adequately  without  the  help  of  active  control  can  be 
packaged  and  launched  in  a  single  shuttle  cargo.  It  should  be  realized  that 
the  antenna  performance  relying  on  passive  means  of  control  only  translates 
naturally  into  increased  stiffness  and  thus,  additional  mass.  In  view  of 
the  large  reflector  size  and  relatively  heavy  structure,  any  packaging  tech¬ 
nique  becomes  a  challenging  problem. 

While  a  complete  packaging  and  deployment  procedure  can  not  be  presented 
here,  an  attempt  is  made  to  show  how  it  is  indeed  possible  to  stow  the  large 
MMW  antenna  structures  in  a  single  shuttle  cargo  bay  (the  volume  available 
is  a  cylinder  18.28  m  long  and  4.57  m  in  diameter).  A  procedure  for  stowing 
the  124  reflector  panels  and  the  packaging  of  the  back-up  trusses  and  other 
satellite  components  in  the  shuttle  cargo  bay  is  suggested  below. 

The  reflector  panels  have  been  designed  as  thin  curved  plates  with  a 
hexagonal  shape  so  as  to  most  efficiently  use  the  circular  shape  of  the 
cargo  bay.  When  stowed  back-to-back  as  "stacked  wafers"  the  124  panels,  6  cm 
in  thickness,  fill  up  7.80  m  of  the  cylindrical  bay.  The  first  panel  in 
the  stack  takes  up  another  8  cm  due  to  its  curvature  and  a  spacing  of  0.2  cm 
is  allowed  between  each  plate. 

The  remaining  volume  is  available  for  the  back-up  trusses,  the  tower 
frame  and  the  various  satellite  components  and  hardware.  An  estimate  of  the 
volume  required  by  the  packaging  of  truss  sections  indicates  that  it  would  fit 
well  within  the  remaining  10.5  m  of  the  cargo.  Unlike  some  folding  techniques 
requiring  hinges  for  the  collapsible  trusses  at  every  member  end,  the  tech¬ 
nique  recommended  attempts  to  limit  to  a  minimum  the  number  of  movable  joints. 
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Table  7-1.  Structural,  packaging  and  deployment  issues 
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This  can  be  reasonably  achieved  by  using  thin-wall  open  cross-sections  for 
all  the  truss  members.  The  shape  selected  is  an  angle  with  small  lips  in 
order  to  ensure  a  reasonable  factor  of  safety  against  instability.  Whereas 
the  traditional  tube,  (the  most  efficient  axial  member)  is  difficult  to  pack 
efficiently,  the  open  section  lends  itself  to  a  high  packaging  ratio.  While 
the  first  truss  section  has  a  width  of  5  cm,  each  additional  truss  can  be 
nested  into  the  previous  one,  adding  only  the  thin-wall  thickness  of  0.2  cm 
and  some  0.2  cm  clearance  to  the  package.  In  order  to  give  an  idea  of  the 
volume  required  to  stow  the  entire  truss  work,  the  longitudinal  and  trans¬ 
verse  1.5  m  deep  trusses  are  assumed  to  be  cut  into  10.5  m  long  sections. 
Approximately  100  pieces  are  needed  so  that  they  can  fit  in  a  parallelepiped 
of  1.5  m  x  0.40  m  *  10.5  m.  The  remaining  top  and  bottom  diagonals  could  be 
stacked  as  individual  members,  thus  requiring  very  little  space.  To  be 
conservative,  however,  it  was  assumed  that  the  top  and  bottom  faces  were  made 
of  3.9  m  deep  horizontal  trusses.  Using  the  same  approach  to  the  packaging 
of  these  horizontal  trusses,  the  required  volume  has  dimensions  of  3.9  m  * 

0.40  m  x  10.5  m. 

The  members  forming  the  tower  frames  are  thin-wall  tubes,  10  cm  in 
diameter.  Because  of  their  span,  these  members  have  been  designed  as  tubes 
to  exhibit  good  bending  capabilities.  Although  many  simple  folding  techniques 
can  be  advocated,  the  tubes  are  simply  shown  stacked  up  in  some  of  the  re¬ 
maining  available  space  of  the  cargo. 

Figure  7-1  shows  the  relative  position  of  the  satellite  components  stowed 
in  the  shuttle  cargo  bay.  It  can  be  visualized  that  the  suggested  folding 
procedure  still  leaves  ample  room  for  the  additional  hardware  such  as  the 
receiver  array  package,  the  communication  antennas,  the  thermal  blankets,  etc. 
Although  it  is  possible  that  a  more  complete  study  of  this  packaging  approach 
could  require  slightly  more  volume  to  accommoda-e  various  connection  hinges 
and  locking  mechanisms,  the  additional  space  mat  would  be  required  is  avail¬ 
able.  It  is,  therefore,  possible  to  assert  that  the  MMW  antenna  structure  can 
be  packaged  into  a  single  shuttle  cargo  bay. 

The  use  of  thin-wall  open  sections  for  the  truss  members  presents  a  number 
cf  interesting  possibilities  for  efficient  packaging  of  trusses.  Among  those 
considered.  Figure  7-2  illustrates  the  possible  folding  of  an  incomplete  space 
truss  in  which  some  folded  diagonal  members  make  efficient  use  of  the  avail¬ 
able  space  between  members  of  the  fixed  trusses.  Although  this  procedure  is 
preliminary  and  has  shortcomings  in  terms  of  structural  efficiency  (e.g., 
eccentricities)  it  is  worth  pursuing  since  it  may  lead  to  a  more  automatic 
deployment  while  minimizing  the  number  of  hinges  ana  ensuring  a  good  packaging 
ratio. 
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SECTION  8 


CONCLUSIONS 


The  objective  of  this  study  was  to  design  an  MMW  antenna  with  a  struc¬ 
ture  that  would  perform  adequately  under  several  constraints  imposed  by  vari¬ 
ous  sources.  The  most  important  constraints  arc  associated  with: 

(1)  Orbital  environment  disturbances. 

(2)  Use  of  passive  structural  systems  only. 

(3)  Packageability  in  a  single  shuttle  cargo  bay. 

The  last  two  constraints  derive  from  desirable  requirements  involving  simpli¬ 
city,  efficiency,  and  cost. 

The  design  of  the  structural  components  was  suggested  on  the  basis  of 
past  experience  in  dealing  with  a  satellite  structure  of  a  given  size  (60  m  * 

30  m  reflector)  associated  with  a  stipulated  deformational  tolerance  (on  the 
order  of  0.5  mm).  A  fundamental  frequency  of  1  Hz  for  the  overall  structure 
was  assumed  to  impart  to  the  satellite  enough  stiffness  such  that  it  would 
distort  within  tolerable  limits  when  exposed  to  the  typical  environment  of  a 
500  km  orbit.  But,  it  is  only  after  going  through  the  entire  analysis  cycle 
that  the  overall  merits  of  the  system  can  be  judged.  The  structure  suggested 
for  this  analysis  appears  slightly  stiffer  than  required  when  exposed  to  the 
expected  dynamic  disturbances.  Modifications  in  design  sizes  and  cross- 
sectional  areas,  with  corresponding  savings  in  masses  can  be  explored,  but 
other  limitations,  such  as  those  required  by  structural  stability  may  possi¬ 
bly  become  critical  in  the  process.  The  study  showed  that  the  most  severe 
loadings  are  those  associated  with  the  thermal  variations  of  the  satellite 
structures  due  co  the  thermal  fluxes  from  the  sun,  the  earth  IR  and  its  al¬ 
bedo.  The  design,  therefore,  must  rel>  on  a  material  exhibiting  a  favorable 
coefficient  of  thermal  expansion  such  as  a  graphite-epoxy  composite  (a  = 

0.2  x  10  6  m/m°k),  and  the  presence  of  a  thermal  shield  capable  of  reducing 
the  temperature  excursions  and  gradients  to  an  acceptable  level.  It  was 
shown  that  it  is  possible  with  state-of-the-art  materials  and  thermal  blanket 
design  techniques  to  keep  the  large  reflector  of  the  MMW  antenna  from  distort¬ 
ing  beyond  specified  values  in  the  changing  thermal  environment. 

In  the  trade-off  studies  for  the  design  of  the  reflector  panels  it  was 
found  that  a  graphite-epoxy  (G-E)  egg-crate  structure  is  by  far  superior  to 
an  aluminum  honeycomb,  mainly  because  of  the  large  difference  between  the  CTE 

(23  x  10  6  m/m°k  for  aluminum  versus  0.2  x  10  0  m/m°k  for  G-E).  The 
possibilities  offered  by  inflacables  such  as  airmats  should  he  pursued  for 
the  design  of  a  reflector  panel.  In  this  study,  a  tentative  design  showed  un¬ 
acceptable  thermal  deformations,  but  there  is  room  for  optimization,  especially 
since  the  restrictions  on  panel  depth  that  are  critical  for  rigid  designs, 
are  almost  negligible  for  inflatables.  The  receiver  tower,  however,  showed 
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thermal  deformations  that  can  cause  an  unacceptable  LOS  error.  This  fact  was 
expected  because  no  attempt  was  made  to  use  a  thermal  shield  for  the  tower. 

It  can  be  anticipated  that  with  thermal  control,  these  deformations  could  be 
brought  to  a  tolerable  level.  Alternatively,  a  semi-active  control  approach 
could  be  used,  based  on  the  predictability  of  the  slowly  changing  thermal 
conditions,  to  compensate  for  unacceptable  feed  movements  as  the  satellite 
progresses  along  its  orbital  path. 

As  evidenced  by  the  requirements  affecting  the  design  of  the  reflector 
panel,  the  issues  of  packaging  and  deployment  have  strongly  affected  the  de¬ 
sign  of  structural  components.  Segmentation  of  the  large  antenna  surface 
was  necessary  and  the  hexagonal  panels  emerged  as  the  best  solution  for 
high  packaging  ratio  and  efficient  surface  assemblage.  Since  they  fill  up 
about  42  percent  of  the  shuttle  cargo  bay,  the  stowing  of  the  remaining  back¬ 
up  trusses  and  30  m  tower  frame  required  efficient  folding  techniques.  An 
attractive  approach  was  suggested  for  tne  truss  members,  where  thir.-wall 
open  cross  sections  replaced  the  traditional  circular  tube.  Once  acceptable 
factors  of  safety  on  lateral  stability  of  individual  members  have  been  made, 
nesting  of  consecutive  trusses  that  can  improve  the  packaging  ratio  of  a 
factor  of  approximately  10  can  be  considered.  While  a  complete  packaging 
method  has  not  been  investigated  for  the  back-up  trusses,  the  above  approach 
presents  a  number  of  interesting  possibilities  and  should  be  pursued.  On  the 
basis  of  rough,  but  conservative  estimates,  the  packaging  of  the  entire  an¬ 
tenna  (including  equipment)  into  a  single  shuttle  cargo  bay  is  indeed 
possible. 

The  study  of  the  representative  design  for  the  MMW  antenna  showed  that 
a  reasonable  performance  can  be  expected  from  a  structural  concept  using  tra¬ 
ditional,  state-of-the-art  material  and  design  techniques.  Of  special  inter¬ 
est  is  the  simplicity  of  the  system  that  can  operate  with  a  minimum  of  con¬ 
trol,  thus  increasing  reliability.  It  should  be  noted  that  the  thermal 
blanket  on  the  tack  face  of  the  supporting  trusses  is  a  purely  passive  means 
of  controlling  the  temperatures.  While  it  is  anticipated  that  kinematic 
mount  actuators  for  panels  are  indispensable  in  the  surface  acquisition  phase 
of  deployment,  they  should  not  be  necessary  during  operation  thereafter. 

The  "exact"  evaluation  of  the  performance  of  the  MMW  antenna  was  ob¬ 
tained  through  the  processing  of  the  information  on  the  deformational  state 
of  the  entire  structure.  In  this  study,  PSF  arid  LOS  motions  have  been  cal¬ 
culated  for  a  number  of  design  cases.  These  values  must  be  assessed  by  MMW 
antenna  specialists  so  that  performance  improvements  through  structural  modi¬ 
fications  or  additional  control  can  be  evaluated.  But,  as  an  indirect  measure 
of  performance,  deformational  tolerances  also  have  been  set  for  the  rms  of 
reflector  surface  deviations  and  the  motion  of  the  feed  support  (decenters 
and  defocus).  As  mentioned  above,  the  merits  of  the  structural  design  have 
been  assessed  mainly  by  comparing  these  deformations  to  the  tolerances  and 
the  conclusions  based  on  this  comparison. 

The  desirable  simplicity  of  the  system  hinges  on  a  relatively  stiff 
structure  and,  therefore,  translates  into  a  substantial  amount  of  mass. 

(The  entire  satellite  is  expected  to  have  a  mass  of  8800  kg.)  This  approach 
using  a  "heavy  and  stiff"  structure  can  be  contrasted  with  another  solution 
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involving  a  "light  and  very  flexible"  structure  that  derives  its  necessary 
stiffness  from  active  control  means.  Using  an  electrostatically  shaped  mem- 
brone  as  a  reflector  is  an  attractive  solution  with  regard  to  mass,  but  re¬ 
quires  a  sophisticated  control  system  to  maintain  the  surface  accuracy  of  the 
flexible  membrane.  The  trade-off  between  the  two  systems  must  involve  a 
large  number  of  characteristics  that  are  sometimes  difficult  to  assess. 

While  the  total  mass  and  stowed  volume  are  easy  to  compare,  the  complexity  of 
the  control  system  and  its  reliability  are  examples  of  features  less  amenable 
tc  resolution. 

The  choice  of  a  500  km  orbit  for  the  design  study  was  somewhat  arbitrary. 
Although  the  atmospheric  density  decreases  drastically  above  an  altitude  of 
300  km  and  usually  does  not  present  a  critical  problem  for  small  satellites, 
the  large  size  of  the  KMW  antenna  poses  an  unusual  problem  for  its  orbital 
dynamics.  The  atmospheric  drag  forces  and  associated  torques  cause  a  rela¬ 
tively  serious  orbital  decay  (in  the  order  of  10  km  in  a  30-day  period)  and 
pointing  disturbances  must  be  compensated  for.  Based  on  a  5-year  life  span, 
the  need  for  approximately  2000  kg  in  propellant  to  maintain  operations  can¬ 
not  be  overlooked.  A  higher  orbit  could  be  considered  to  alleviate  this 
problem.  The  present  design  for  use  &-  this  altitude  should  be  revised  or 
modified  to  include  structures  needed  for  the  support  of  this  additional  mass. 
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Aba bract 

Villa  paper  describes  a  method  for  evaluating  the  optical  performance  of  opt  leal  systems 
consisting  of  largo,  flexible,  actively  controlled  mirrors  in  the  presence  of  thermal  and 
dynamic  disturbances.  The  computer  simulation  described  here  Is  divided  into  four  major 
blockai  optica,  thermal,  a true Lures ,  and  dynamics.  The  control  aspects  are  included  :n  the 
various  models  uaod  i,t  each  block.  The  output  of  the  simulation  is  a  dynamic  point' spread 
function  from  which  real-time  or  average  quantities  auch  as  the  l iiu—  of-s ight  error,  Strehl 
rat  io,  and  power  in  a  certain  area  are  caleul  tad.  As  an  example,  a  two-mirror  system  that; 
expands  a  htgh-powor  collimated  beat,  and  focuses  it  on  a  distant  plane  is  considered.  The 
effects  of  heating  of  the  mirrors  are  simulated  and  results  on  the  point-sproau  function  and 
associated  quantities  are  presented. 


Lr  introduction 

A  computer  simulation  for  evaluating  the  optical  performance  of  largo  adaptive  optical 
systems  -a  desr«-ibed.  Given  the  optical  system  design  data  (e.g,,  mirror  figures,  spacingo, 
and  ante..),  optical  ray-' race  programs  are  used  to  obtain  the  angles  and  the  points  of  inci¬ 
dence  of  *  he  rays  on  the  mirrors  and  the  wave  aberrations  of  tho  system.  The  deformations 
of  the  optical  system  (mirror  warping  and  misalignments)  due  lc>  some  disturbance  (e.g., 
thermal  loading,  or  a  dynamic  disturbance)  are  calculated  by  structural  programs -such  as 
f'TAKDYNH  and  NASTKAh ,  Any  control  aspects  of  the  system  are  included  in  this  calculation, 
the  residual  deformations  u.e  given  at  discrete  points  on  the  optical  elements,  which  gen¬ 
erally  .-o  not  the  same  as  the  points  of  incidence  of  the  rays  used  in  optical  ray-trace 
progra  ...  in  the  simulation,  the  do  format  ions  obtained  at  structural  points  are  interpo¬ 
lated  to  obtain  tho  deformations  at  the  optical  poihts. 

Tho  wave  abe-rat ions  produced  by  these  dofoimat ions  are  calculated  using  (he  ray- trace 
dal,'  of  the  mule  formed  system,  and  are  added  to  the  dial  urbance- free  aberrations  of  the  sys¬ 
tem  to  obtain  tho  total  wave  aberrations  ot  the  deformed  system.  A  pupil  rune Lion  is  formed 
irom  those  total  aberrations,  and  using  fast  Fourier  transform  (ITT)  programs,  the  point- 
spread  {unction  (PSD  ir.  calculated.  The  amount  of  image  energy  m  a  certain  area,  the  cen¬ 
troid  of  the  imago,  and  line-of-night  (l.OS)  or,-or  arc  obtained  from  the  point-spread  func¬ 
tion.  m  the  case  of  a  dynamic  disturbance,  both  the  dynamic  and  time -averaged  l’RF  and 
related  quantities  can  be  obtained.  Vne  image  of  an  i neoheront  extended  obj-set  lying  in  an 
isopl.tn.iLic  region  ol  the  system  can  be  obtained  by  convolving  tho  object  distribution  witli 
the  corresponding  PSP.  However,  images  of  extended  objects  are  not  considered  in  Utis  paper. 

An  example  of  a  two 'mirror  system  that  expands  a  high-power  collimated  beam  and  focuses 

it  on  a  distant  plane  is  considered.  The  effects  of  heating  of  the  mirrors  by  ‘.lie  beam  are 

simulated  and  results  on  tin.  PSP,  Strehl  ratio,  and  ensqu.nvd  power  arc  presented 

i  ■  Hist  ur!  i  a  nee  -  Pr  e  q  _  System 

Consider  a  point  obi eel  at  an  angle  a  from  the  optical  axis  of  an  imaging  system  consist¬ 
ing  of  m  mirrors.  The  optical  performance  of  the  system  is  evaluated  by  tracing  a  certain 
number  of  raja.  These  rays,  which  represent  light  emanating  from  the  point,  object,  are  dis¬ 
tributed  uniformly  in  a  square  at  ray  in  the  entrance  pupil  of  the  system.  Depending  on  Die 
sjko  of  the  entrance  pupil  and  its  associated  obscurations,  only  some  of  the  rays  are  i rans- 
mi t.'.e-.i  by  the  system.  bet  N  be  tile  number  of  transmitted  rays. 

The  rays  are  traced  through  the  system  striking  the  mirrors  and  finally  reaching  the  exit 
pupil  of  the  system.  Wu  find  the  end  points  of  all  Lite  rays  so  that  the  total  optical  path 

length  travelled  by  each  of  them  is  the  same.  For  example,  the  pat.lt  length  of  ail  Lite  lays 

may  be  equal  to  that  of  llio  ray  which  starts  and  cuds  at  the  centers  of  the  entrance  and 
exit  pupils  Ueeauso  of  obseu rat  ions  in  mirror  systems,  this  ray  is  usually  not  physically 
transmitted  by  the  system,  but  its  optical  path  length  can  be  '\i  leu  l  a  led  as  if  it,  wore. 

Tho  surface  passing  through  the  end  points  of  the  rays  is  called  tho  system  wavefront  for 
the  point  object  under  e  •nsideration,  i.e.,  wavefront,  for  a  field  angle  a.  If  tho  end 
points  lie  on  a  spherical  surface,  then  we  have  a  spherical  wavefront  and  the  system  us  said 
to  ho  aberration- free  for  that  field  angle.  The  extent  or  tho  boundary  of  tho  wavefront  is 
defined  by  the  rays  touching  the  boundary  of  the  exit  pupil.  The  center  of  curvature  of  the 


0  Tlte  Charles  Stark  Draper  Laboratory ,  Inc.,  1 8 ') . 


si'll:  t'.i /,  l-V  histniinetitatton  in  AsUonamy  III  119791 1  459 


mmsmnmwtn 


MAHAJAN.  GOVIGNON 

apht  wavefront  gives  the  focus  of  the  system  for  field  angle  a,  i.e.,  the  center  of 

the  i*: ’fiance  distribution.  The  irradiance  distribution  is  called  the  PSF. 

If  the  end  points  of  the  rays  do  not  lie  on  a  spherical  surface,  we  determine  a  spherical 
surface  passing  through  the  center  of  the  exit  pupil  such  that  the  path-lenath  deviations 
along  the  rays  from  it  have  a  minimum  variance.  These  deviations  are  called  wave  aberra¬ 
tions  or  wavefront  errors,  and  their  sign  convention  is  such  that  the  aberration  correspond¬ 
ing  to  a  certain  ray  is  positive  if  it  has  to  travel  an  extra  amount,  to  reach  the  spherical 
surface.  The  minimum  variance  spherical  surface  is  called  the  reference  sphere  and  its  cen¬ 
ter  of  curvature  is  called  the  best  diffraction  focus.  Let  the  best  diffraction  focus  be 
located  at  rf (x&, y£, z£;a) .  The  surface  passing  through  the  best  focii  for  various  field 
ancjles  (point-object  positions)  is  called  the  focal  or  the  image  surface. 

Let^ (u,v,w)  be  the^coordinates  of  a  point  where  a  ray  intersects  the  reference  sphere  such 
that  u2  +  v2  +  w2  =  R2 ,  where  R  is  the  radius  of  curvature  of  the  reference  sphere.  If 
W4(u,v;o)  is  the  wave  aberration  and  A(u,v;ci)  is  the  amplitude**  at  a  point  (u,v)  in  the 
pupil,  then  the  pupil  function  of  the  system  is  defined  by 


P(u,v;a)  = 


"’mi  » 

A(u,v;a)  exp  -4—  W°(u,v;a )  ,  in  the  clear  region 


of  the  pupil 


=  0,  otherwise 


the  pupil  function  is  nonzero  only  in  the  region  of  transmitted  rays.  The  amplitude 
variations  in  the  pupil  can  be  due  to  a  number  of  sources,  e.g.,  the  input  wave  to  the  sys¬ 
tem  has  amplitude  variation  or  the  mirrors  have  nonuniform  reflectivities.  In  Eq .  (1),  X  is 
the  wavelength  of  the  object  radiation  and  i  =  /-T. 

The  irradiance  or  PSF  value  at  a  point  (x,y;0>  with  respect  to  rf,  where  the  (x.y) 

■'lane  is  parallel  to  the  (u,v)  plane,  is  given  by! 


I  (x,y;a) 


t// 


P(u,v;a)  exp  - 


(xu  +  yv)  du  dv! 


The  integral  in  Eq.  (2)  is  a  Fourier  transform  of  the  pupil  function.  For  small  angular 
deviations  from  angle  a,  the  pupil  function  and  therefore  the  PSF  do  not  change,  except  for 
an  overall  displacement.  The  angular  extent  of  the  object  over  which  the  PSF  does  not 
change  significantly  is  called  an  j.soplanatic  angle  of  the  system.  It  is  evident  that 
in  a  computer  simulation,  a  discrete  Fourier  transform  is  used,  particularly  because  the 
aberrations  are  known  at  a  set  of  discrete  points,  the  points  at  which  the  rays  end  in  the 
exit  pupil.  An  FFT  program  is  used  to  evaluate  the  integral  in  Eq.  (2). 

Let  lu  and  Lv  be  the  spatial  periods  o*  the  aberration  data  points  along  the  u  and  v 
axes,  respectively.  The  coordinates  (u,"'  of  a  point  in  the  pupil  can  be  written 

u  =  i..u  v  =  k_v  (3) 


where  j  and  k  are  integers.  The  discrete  form  of  the  pupil  function  in  the  clear  region 
of  the  pupil  can  be  written 


P  ( u ,  v ;  :> )  =  P ( J2u ,k  .vs  .) 


o  - .  . 

=  A ( g  u,k  v;  .!  exp  =7-=-  ( j 'u.k ’ v;  «) 


It  is  assumed  throughout  that  the  aberrations  arc  small  so  that  minimization  of  their 
variance  leads  to  a  better  quality  image,  and  the  center  of  the  reference  sphere  is  a 
point  of  maximum  irradiance.  In  the  case  of  large  aberrations,  geometrical-optics-based 
criteria,  such  as  the  smallest  spot  diameter,  may  be  used.  In  that  case  there  is  no  dis¬ 
tinguishable  principal  maximum  in  the  PST,  and  the  LOS  may  be  determined  by  the  centroid 
of  the  image  rather  than  the  diffraction  focus. 

## 

For  a  point  object  at  a  distance  z>  from  an  imaging  system  with  a  uniformly  illuminated 
pupil,  A(u,v;a)  is  constant  over  the  pupil  and  equal  to  (Pj/4")  Vz.i.M,  where  P.  is  the 
total  power  radiated  uniformly  by  the  object  m  all  directions  and  M  is  the  pupil  magni¬ 
fication,  e.g.,  M  is  the  ratio  of  the  diameters  of  the  exit  and' entrance  pupils. 
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Or 


P (u, v;a) 


Ad  <a)A^k  (a)  exp 


2^W»k(a) 


A,  (a)Pjk(a) 


(4b) 


where  A0  is  a  normalization  amplitude  factor  so  that  A3k(a)  and  therefore  Pjk(a)  are  dimen¬ 
sionless  quantities.  The  A3k(a)  may  be  chosen  such  that,  for  example,  the  maximum  value  at 
a  certain  point  such  as  the  center  of  the  pupil  is  unity. 

If  Pt(a)  is  the  total  power  entering  the  pupil,  it  is  given  a y 


Pt(a)  =  JJ  |P(u,v;a)  | 2  du  dv  (5) 

pupil 


where  the  integration  is  carried  over  the  clear  region  of  the  pupil.  Substituting  Eq.  (4) 
into  Eq.  (S) ,  a  discrete  form  of  Eq.  (5)  can  be  written 


Pfc(a)  =  Ao  (a) AuAv  wj)J  P  jk  la>  I  * 

j,k 

=  A§  (a)  AuAv  ^  wjkAjk*a'  (6) 

j.k 


where  w.,k  are  the  weights  of  the  integration  rule  (e.g.,  trapezoidal,  Simpson's,  etc.). 

For  efficient  use  of  the  FFT  programs,  data  points  with  zero  pupil  function  are  added 
in  the  pupil  plane,  such  that  the  number  of  points  along  the  u  and  v  axes  are  each  equal 
to  certain  powers  of  2.  Thus,  for  example,  a  grid  of  NuxNv  points  is  generated,  where  nu 
and  Nv  are  the  numbers  of  rows  and  columns,  respectively,  of  the  data-point  matrix,  both 
Nu  and  Nv  being  equal  to  powers  of  2.  Of  the  NuxNv  data  points,  only  N  have  a  nonzero  value 
of  the  pupil  function. 

Let  Ax  and  Ay  be  spatial  periods  of  data  points  in  the  image  plane  along  the  x  and  y 
axes,  respectively.  Thus,  we  may  write 


x  =  £Ax 
y  =  mAy 


(7) 


where  i  and  m  are  integers.  Substituting  Eq.  (3),  (4),  and  (7)  into  Eq.  (2),  a  discrete 
form  of  Eq.  (2)  can  be  written 


I (£Ax,mAy;a) 


A§  (a) 


T#  K 


exp 


-  (jUAxAu  +  kmAyAv) 


(8) 


Buth  (£,m)  and  (j,k)  pairs  form  a  i#u«Nv  matrix  of  data  points.  In  our  FFT  program,  the  PSF 
value  at  an  (t,m)  data  point  is  given  by 


l£m(a)  =  d-;Epjk(a)  e*pi 

•j,k 


-2iu| 


111  +  huA 

k  \) 


(9) 


Comparing  Eq.  (8)  and  (9) ,  we  find  that 


IUAx,mAy;a)  =  A5(a)NJIJm(a) 


(10) 
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Moreover,  the  spatial  periods  in  the  pupil  and  image  planes  are  related  according  to 


Ax  * 

XR/N^AU 

Ay  * 

XR/Nv*v 

(11) 

For  an  exit  pupil  with  extents  Du  and  Dv  along  the  u  and  v  axes,  respectively,  Eq.  (11) 
can  be  written 


A*  =  XFunu/Nu 

iy  -  X£VVNv  (12> 

where,  for  example 

Fu  =  R/Du  (13) 


is  the  focal  ratio  of  the  image  along  the  u  axis,  and 


nu  »  Du/iu  U-4 

is  the  number  of  rays  in  a  length  of  0U.  Note  that  to  reduce  the  spatial  periods  in  the 
image  plane,  the  matrix  dimensions  Nu  and  Nv  must  be  increased.  Thus,  by  merely  adding 
data  points  with  zero  values  in  the  pupil  plane,  a  closely  spaced  distribution  of  PSF 
values  can  be  obtained. 

Frcra  Eq.  (9)  and  (10),  the  central  PSF  value  can  be  written 

ko, o;c)  =  (^jrf a§(a)|£  V-r  u5a) 

!  3  rk 

Substituting  for  A3 (a)  from  Eq.  (6),  it  can  also  be  written 


\T  P,kU>!- 

P  (u)iuiv  If-  3k  ! 

I  (0 , 0 ;  a)  =  — -  -  (15b) 

<XR>:  E  wjkAjk(a) 

i.k 


The  fraction  of  power  m  a  certain  area  of  the  image  plane  can  be  obta-r.ed  by  integrating 
the  -*SF  over  that  area.  For  example,  if  a  detector  is  placed  in  the  imay  plane,  the  frac¬ 
tion  of  total  image  power  collected  by  the  detector  is  given  by 


paet(*>  -  njr//I(s^>  (l6) 

r  det 


Since  the  aberration  function  W  is  given  at  a  uniformly  spaced  set  of  points,  the  PSF 
data  is  also  obtained  at  a  uniformly  spaced  set  of  points.  Therefore,  it  is  convenient  to 
calculate  the  onsquared  power,  i.e.,  the  power  on  a  square  (or  rectangular)  detector.  For 
a  circular  detector,  the  PSF  data  has  to  be  interpolated  to  obtain  encircled  power.  We  have 
shown  that  for  diffraction-limited  imaging  systems  with  centrally  obscured  circular  pupils, 
the  difference  between  corresponding  ensquared  and  encircled  powers  (square  having  the  same 
half-width  as  the  radius  of  the  circle)  is  less  than  10%  of  the  total  power,  regardless  of 
the  amount  of  obscuration  or  the  size  of  the  detector. : 
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Using  Eq.  (7),  a  discrete  form  of  Eq.  (16)  can  be  written 


Pdet(a) 


det 

E  UAx.mAyja) 

£,m 


(17) 


where  w,m  are  the  weights  of  the  integration  rule  used,  and  only  those  (i,m)  data  points 
are  used  in  the  summation  that  lie  on  the  detector.  Substituting  Eq.  (6),  (10),  and  (11) 
into  Eq.  (17) ,  we  obtain 


Pdet('J) 


u} 


det 

E  w?.mI?,m(a) 


f,  ,m 


N  N 
u  v 


E 

j.k 


wjkAjk(a> 


(18) 


If  the  light  entering  the  exit  pupil  has  uniform  phase  and  amplitude,  l.c.,  for  a  spheri¬ 
cal  wavefront  with  uniform  amplitude  across  it,  we  may  write  using  script  letters  for  the 


pupil  function,  PSF,  and  related  quantities 

i’(u,v;<) 

= 

A  Q 

(19a) 

»V‘> 

= 

V‘) 

= 

1  (at  N  (gk)  points  only) 

(19b) 

/  m<->  -  -iE 

exp 

■*)!!' 

(20) 

1  3<k 

1  (■  ‘x,m  y;  ,)  =  | 

NAuAv  ,  ,  , 

XR  J  A»(,,/.m(,) 

(21) 

A  V (  .)  *  l-  (  ,)/N'.u  v 


'’det<!) 


det 

fnr  £  w.m?.m<‘) 


(22) 


(23) 


/  <  0 , 0 ;  .) 


N/’t  (  <)  '  u  ‘  v 
(  R) 


(24a) 


(24b) 


Assuming  that  t ho  total  power  entering  the  pupil  when  the  system  wavefront  is  spherical  and 
the  amplitude  across  if  is  uniform  is  the  same  os  when  there  are  aberrations  and  amplitude- 
variations,  we  obtain  the  StrehL  ratio*  of  the  degraded  image  from  iiq.  (15b)  and  (24b) 

* 

The  SLrehl  ratio  is  generally  J.-ined  as  Lac  ratio  of  central  values  of  a  PSP  in  the  pres¬ 
ence  and  absence  of  aberration-  in  a  s.stom.  Here,,  wo  generalise  t.ns  definition  to  in¬ 
clude  Lae  degradation  of  central  irradiance  din-  to  aniplilade  vjrul,  ms  as  well. 
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which  is  less  than  1. 


!E  Va>  2 

_  1  (O.O)g)  _  |j,k _ 

‘  N  E  W«> 

l.k 


Disturbed  (Deformed)  System 


Now  suppose  that  the  imaging  system  is  disturbed  and  the  mirrors  are  deformed.  Assume 
that  the  disturbance  is  time  dependent/  although  it  may  contain  a  time-independent  component, 
e.g.,  a  dynamic  disturbance  in  the  presence  of  a  slowly  varying  thermal'  disturbance.  Con¬ 
sider  the  state  of  the  optical  system  at  a  certain  time  t,  assuming  that  the  control  system 
has  restored  the  system  to  the  best  of  its  ability. 

Suppose  that  a  certain  point  rim(a)  on  the  surface  of  the  mth  mirror  touched  by  the  ith 
ray  is  moved  by  the  disturbance  and  restored  by  control  to  a  point  rjn(o,t),  so  that  the 
residual  vector  displacement  of  the  point  at  a  time  t  is  given  by 

Sr._(a;t)  =  r!  (a;t)  -  r.  (a)  (26) 

im  im  im 

For  small  5r,  the  corresponding  wave  aberration  contribution  to  the  ith  ray  is  given  by 
(see  the  Appendix) 


5Wim(a;t)  =  2lniB(x) 


5r  (a;  t)  ]  cos  9  - {a ) 
in\  im 


where  n,m  is  a  unit  vector  along  the  normal  to  the  mirror  surface  at  the  point  r1;n,  and  Si!n 
is  the  angle  of  incidence  of  the  ray  on  the  undeformed  mirror.  The  total  aberration  con¬ 
tributed  by  all  the  mirrors  (because  of  their  deformations)  to  the  ith  ray  is  obtained  by 
summing  their  contributions,  i.e. 


;w(ui(v. jo;t)  = 


E  6Wim(‘”t> 


The  total  aberration  Wi  at  a  point  (u1(Vi;g)  in  the  exit  pupil  is  obtained  by  adding  GW,  to 
the  disturbance-free  aberration  wj ,  i.e. 


Wfu^  ,v^;ci;t)  =  W5  (ui,vi;g;t)  +  IWtu^v^joijt) 


Given  the  total  aberrations  Wt,  once  again  we  find  a  reference  sphere  such  that  when  the 
aberrations  are  measured  with  respect  to  it,  their  variance  is  minimized.  If  the  center  of 
curvature  of  this  reference  sphere  is  located  at  tp(.<:t! ,  the  amount  of  defocus  caused  by 
the  disturbance  is  given  by 


lrF(i;t)  =  rF(a;t)  -  r°(a) 


=  llXphit)  ,.".yF(c«;t)  , iZp! •« ; t)  ! 


where  [ixr (a ;  t) ,  lyF(o;t)l  is  the  transverse  defocus  and  izr(a;t)  is  the  longitudinal  de¬ 
focus.  The  transverse  defocus  represents  the  LOS  error. 

Let  Wf (uj^v^ajt)  be  the  wave  aberration  at  a  point  (ul,v1;i)  in  the  exit  pupil  with  re¬ 
spect  to  the  reference  sphere  centered  at  rr(."<;t).  Using  these  aberrations,  once  again  we 
go  through  the  equations  of  Section  2  to  determine  the  FSF  and  ensquarea  power.  By  deter¬ 
mining  P (a; t)  at  different  times,  -we  can  calculate  the  average  power  or  total  energy  in  a 
certain  area  during  a  certain  tine  interval. 
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^to  that  the  use  of  the  wave  aberrations  implies  that  the  irradiance  distribution  is 
cent' red  about  the  point  rp(a;t).  This  implies  that  there  is  an  LOS  error  and  longitudinal 
defocus  which  are  being  corrected  in  real  time  by  the  control  system.  If  they  are  not  cor¬ 
rected,  the  aberrations  W,  of  Eg.  (29)  must  be  used,  i.e.,  aberrations  with  respect  to  a 
reference  sphere  centered  at  (xp,yF,zF)  rather  than  at  (Xp,yp,zF;t)  must  be  used.  If 
tile  detector  surface  for  the  field  of  view  of  the  imaging  system  is  planar  rather  than 
curved  (as  determined  by  the  locus  of  the  best  diffraction  focii) ,  as  is  often  the  case 
(e.g.,  films)  then  the  (u,v)  plane  is  chosen  to  be  parallel  to  the  detector  plane  and  the 
reference  spheres  for  various  field  angles  u  must  be  centered  on  the  detector  plane. 


4.  Simulation  Architecture 


In  terms  of  evaluating  the  effects  of  thermal  and  dynamic  disturbances  on  the  optical 
performance  of  a  large  optical  system,  the  simulation  as  indicated  in  Figure  1  is  divided 
into  four  blocks:  thermal,  structures,  dynamics,  and  optics.  Given  a  certain  optical  sys¬ 
tem  with  its  associated  structure,  the  information  on  any  heat  sources  is  used  to  carry  out 
a  thermal  analysis.  The  output  of  this  analysis  is  the  temperature  distribution  within  the 
system.  The  temperature  distribution  is  used  to  determine  the  deformations  of  the  optical 
elements.  From  the  material  and  mechanical  properties  of  the  various  elements,  the  struc¬ 
tural  model  determines  the  eigenmodes  and  influence  matrices  (the  effect  of  a  disturbance 
at  one  point  on  others)  of  the  system  or  a  subsystem.  It  also  determines  the  deformations 
of  the  mirrors  at  selected  points,  which  are  the  nodes  of  the  finite-element  structural 
model,  hereafter  called  structural  points.  The  eigenmodes  in  conjunction  with  the  dynamics 
model  and  the  disturbances  are  used  to  carry  out  a  dynamics  analysis,  to  determine  the  move¬ 
ments  of  the  segments  of  an  optical  element,  or  the  optical  element  as  a  whole,  or  the  sys¬ 
tem  as  a  whole.  The  control  of  figure  or  rigid-body  motion  of  a  segment  or  whole  optical 
element  is  included  in  the  structural  and  dynamics  models. 


From  the  deformations  and  movements  of  the  optical  elements,  the  optics  block  determines 
the  LOS  error,  longitudinal  defocus,  the  image-plane-irradiance  distribution  and  associated 
ouantities.  A  detailed  flow  chart  of  the  optics  block  is  given  in  Figure  2.  Note  that 
from  the  structural  and  optical  points,  interpolation  coefficients  are  obtained,  which  are, 
of  course,  time  independent.  They  are  used  to  calculate  the  deformations  at  the  optical 
points  from  those  at  the  structural  points,  from  Eg.  (27) ,  wo  note  that  we  need  the  deform¬ 
ations  along  the  surface  normals,  i.e.,  n  •  ir,  and  we  interpolate  these  quantities  and  not 
the  deformations  themselves.  This  greatly  simplifies  the  calculations  because  this  way  we 
perform  a  scalar  rather  than  a  vectorial  interpolation. 


SPIB  Vol.  172  Instrumentation  in  Astronomy  III  11973)  /  4*5 


145 


MAHAJAM.  GOVIGNON 


5.  Simulation  Example  and  Results 


5.1.  Disturbance-Free  System 

A  two-mirror  optical  system  has  been  used  as  an  example  for  a  simulation.  As  shewn  in 
Figure  3,  mirrors  Mi  and  M2  have  di-neters  of  Di  and  D; ,  respectively.  The  input  to  the  sys 
tem  is  a  collimated  circular  beam  of  uniform  amplitude,  propagating  parallel  to  the  optical 
axis  of  the  system.  The  beam,  whicn  has  a  linear  obscuration  of  0.5,  is  expanded  by  the 
two-mirror  system  and  is  focused  on  the  image  plane.  The  collimated  beam  is  equivalent  to 
a  point  object  at  infinity,  with  an  additional  element  near  the  two-mirror  system  providing 
obscuration.  The  mirror  figures  and  spacing  are  such  that  they  do  not  introduce  any  aber¬ 
rations,  and  a  perfect  image  is  formed  in  the  image  plane. 

Since  the  input  beam  anc  the  two-mirror  system  are  aberration  free,  and  assuming  that  tne 
system  does  not  introduce  any  amplitude  variations,  the  image-plane  irradiar.ee  distribution 
(?SF)  can  be  calculated  analytically  using 
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where  the  irradiance  at  the  cer.  ,er  of  the 
image  is  given  by 


IiOjc) 


e2)  Pt/4X2F2 


In  Eq.  (31)  and  (32),  r  is  the  distance  of 
a  point  in  the  image  plane  from  the  image 
center,  Jt(*)  is  the  first-order  Bessel 
function  of  the  first  kind,  c,  X,  p  ,  and 
F  are  the  linear  obscuration,  wavelength, 
total  power,  and  focal  ratio  of  the  beam. 

If  D  is  the  diameter  of  the  exit  pupil  and 
R  the  radius  of  the  reference  sphere,  then  F  =  P./D. 


- I  J 

-3*  ± 


Beam  expander  and 
focusing  system. 


In  the  present  case,  D  =  D 


The  irradiance  distribution  in  the  image  plane  normalized  to  its  central  value  is  shown 
in  Figure  4.  This  distribution  is  circularly  symmetric.  In  Figure  4,  the  fraction  of  total 
power  contained  in  a  circular-  area  of  a  certain  radius  and  the  corresponding  power  in  a 
square  area  of  the  same  half-width  as  the  radius  are  also  shown.2  Note  that  the  difference 
between  the  corresponding  ensquared  and  encircled  powers  Ps  and  Pc,  respectively,  which  is 
exaggerated  in  the  figure  by  a  factor  of  10,  is  always  less  than  9%  of  the  total  power. 

The  distribution  of  rays  in  the  exit  pupil  of  the  system,  using  a  41x41  array  of  rays, 
is  shown  in  Figure  5.  Because  the  beam  is  circular  with  a  linear  obscuration  of  0.5,  only 
952  of  the  1681  rays  have  nonzero  amplitude.  The  l*s  in  the  figure  indicate  that  the  ampli¬ 
tude  across  the  beam  is  uniform.  .  The  corresponding  image-plane  irradiance  distribution 
using  a  128x128  FFT  grid  is  shown  in  Figure  6.  The  ensquared  powers  obtained  by  integrating 
such  a  distribution  using  the  trapezoidal  rule  are  listed  in  Table  1.  Semi-analytical  re¬ 
sults  obtained  by  integrating  Eq.  (31)  are  also  listed  in  this  table.  Note  that  in  going 
from  64x64  to  128x128  FFT  grid,  the  difference  between  the  semi-analytical  results  and  those 
obtained  by  use  of  the  FFT  program  reduces  to  less  than  5%. 


5.2  Thermal  Disturbance  (Effects  of  Hi-ror  Heating) 


Assume  that  the  beam  entering  the  two-mirror  system  is  of  high  power.  The  mirrors  M. 
and  h3  are  heated  by  the  beam  and  therefore  deform.  The  deformations  of  mirror  M,  along 
the  surface  normals  are  given  in  Figure  7  at  a  few  typical  time  instants.  Most  of  the 
deformation  is  a  change  in  the  radius  of  curvature  and  displacement  along  the  optical  axis. 
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Fig.  4.  Normalized  image-plane 
irradiance  and  power  distribu¬ 
tions.  Ps  and  Pe  are  fractions 
of  a  total  image  power  in  a 
square  of  half-width  r,  and  a 
circle  of  radius  rp,  respectively. 


Fig.  5.  Distribution  of  rays  in 
the  exit  pupil. 


Fig.  6.  Image-plane 
irradiance  distribu¬ 
tion  obtained  by 
using  a  128x128  FFT 
grid . 
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Table  1.  Ensquared  power  Ps(r0)  as  a  func 
tion  of  half-widts  r c  -in  units  of  AF  or  a 
square  obtained  t/  integrating  Eq.  (31) 
and  by  using  the  FFT  p.  ogram  and  then 
integrating. 


thus  producing  a  defocus.  Since  the  beam 
flux  density  is  circularly  symmetric,  the 
deformation  is  also  circularly  symmetric. 

The  mirror  M.  has  been  modeled  into  12- 
degree  segments^  The  mirror  is  supported 
along  radial  lines  every  12  degrees  by  rods. 

The  deformations  have  be  .n  determined  along 
a  supported  radial  line  and  a  free  radial 
line  between  two  supported  lines.  .Because 
of  the  circular  symmetry  of  the  beam  inci¬ 
dent  on  the  mirror,  deformations  along  all 
supported  lines  are  identical;  similarly, 
deformations  along  all  free  lines  ;ire  iden¬ 
tical.  The  deformations  along  a  supported 
and  a  free  radial  line  are  shown  in  Fig¬ 
ure  3.  The  details  of  the  deformation 
calculations  which  depend  on  the  beam 
.power  and  the  thermal  properties  of  the 
mirrors  are  not  included  here. 

The  wave  aberrations  contributed  by 
these  mirror  deformations  are  shown  in  Fig¬ 
ure  9  on  a  64*64  grid.  As  stated  earlier, 
most  of  the  error  is  defocus,  i.e.,  a 
change  in  the  radius  of  curvature  of  the 

reference  sphere.  Indeed,  as  an  example,  by  refocusing  (e.g.,  by  translating  mirror  M, 
aiong  the  optical  axis) ,  the  standard  deviation  of  the  wavefront  errors  at  23  s  is  reduced 
to  A/33. 

The  image-plane  Arradiance  distribution  as  a  function  of  time  in  the  absence  of  refo¬ 
cusing  is  shown  in  Figure  10.  At  23  s,  for  example,  the  Strehl  ratio  (central  irradiance 
ratio)  drops  to  a  value  of  0.03.  By  refocusing,  it  is  brought  back  to  a  value  o'f  0.96. 

The  variation  of  Strehl  ratio  with  time  is  shown  in  Figure  11.  Whereas  w;th  proper  focus¬ 
ing  the  Screhl  r-tio  remains  nearly  1,  without  active  focusing  it  drops  rapidly  to  very 
small  values,  when  the  aberration  is  a  pure  defocus,  the  Strehl  ratio  for  an  obscured  cir¬ 
cular  beam  is  given  by“ 
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s  =  (sin  {•’Wo  (1  -  e*)/A|/(iw,  (l  -  e:)/\)}: 


(33) 


where  W,  is  the  aberration  at  the  edge  of  the  pupil  (beam) .  Thus  the  Strehl  ratio  ap¬ 
proaches  zero  whenever  W0  =  nA/(l  -  ez)  waves,  where  n  is  an  integer.  In  Figure  11,  the 
Strehl  ratio  does  not  become  zero,  because  the  aberration  is  not  a  pure  defocus.  This  may 
be  seen  particularly  in  the  case  of  mirror  M*»  Figure  3. 
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Figure  7.  formation  of  -xrror  M  . 
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The  percentage  of  total  power  contained  in  a  square  area  as  a  function  of  the  half-width 
of  the  square  without  and  with  focus  correction  is  shown  in  Figures  12  and  13,  respectively. 
As  expected  from  the  preceding  discussion,  the  ensquared  power  decreases  considerably  unless 
the  image  is  properly  focused.  In  Figure  12,  the  case  of  100  s  is  also  shown.  In  that 
case,  the  amount  of  defocus  is  so  large  that  the  square  areas  considered  lie  practically 
inside  the  obscured  region  of  the  beam.  With  focus  correction  (not  shown  in  Figure  13) ,  a 
Strehl  ratio  of  0.78  is  obtained,  indicating  the  presence  of  higher  order  aberrations,  and 
therefore  requiring  figure  correction  of  mirrors. 

Knowing  the  longitudinal  defocus  produced  by  the  disturbance  (from  minimization  of  aber¬ 
ration  variance  and  determination  of  the  centers  of  curvature  of  the  reference  spheres  in 
the  absence  and  presence  of  the  disturbance) ,  the  movement  of  mirror  Mi  along  the  optical 
axis  can  be  determined  in  order  to  keep  the  image  focused  in  the  disturbance-free  image 
plane.  Because  of  the  symmetry  of  the  mirror  deformations,  there  is  no  tr~- averse  defocus, 
i.e.,  no  LOS  change  is  produced  by  the  disturbance  under  consideration. 
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Fig.  8.  Model  and  deformations 
of  mirror  M? .  The  distance  from 
the  renter  of  the  mirror  is  meas¬ 
ured  in  units  of  its  diameter. 
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Fig.  9.  Wavefront  errors  at  various 
time  instants.  Most  of  the  error  is 
defocus. 
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Fig.  10.  Image-plane  distribution  at  var¬ 
ious  time  instants  without  focus  correc¬ 
tion.  S  is  the  SCrenl  ratio  of  the  image. 
The  parenthetical  numbers  refer  to  Strehl 
ratio  obtained  with  focus  correction.  A 
64*64  FFT  grid  is  used. 
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Fig.  11.  Strehl  ratio  as 
a  function  of  time. 


Fig.  12.  Ensquared  power  as 
a  function  of  half-width  of 
the  square  at  various  time 
instants  without  focus 
correction. 


Fig.  13.  Ensquared  power  as 
a  function  of  half-width  of 
the  Square  at  various  time 
instants  with  active  focus 
correction. 


6.  Discussion  and  Conclusions 

Ke  have  described  a  computer  simulation  for  determining  the  optical  performance  of  imag¬ 
ing  systems  consisting  of  large,  flexible,  actively  controlled  mirrors  in  the  presence  of 
thermal  and  dynamic  disturbances.  In  the  simulation,  the  deformations  of  optical  elements 
as  determined  by  structural  analysis  are  converted  into  optical  aberrations  using  optical 
data  of  the  undeformed  system.  Using  FF?  programs,  the  simulation  determines  the  dynamic 
PSFs  of  the  system  and  real-time  or  average  quantities  such  as  the  LOS  error.  Strehl  ratio, 
and  ensquared  power. 

As  an  example,  a  two-mirror  system  that  expands  and  focuses  a  collimated  centrally  ob¬ 
scured  circular  beam  is  considered,  and  the  accuracy  of  the  ensquared  power  calculation  is 
discussed.  It  is  shown  that  the  results  obtained  are  accurate  to  within  a  few  percent. 
L'ext,  the  beam  is  assumed  to  be  of  high  power  so  that  it  heats  and  therefore  deforms  the 
mirrors.  The  optical  aberrations  produced  by  the  deformations  are  calculated  and  it  is 
shewn  that  most  of  the  error  can  be  removed  with  a  focus  correction,  e.g. ,  by  translating 
one  of  the  mirrors  along  the  optical  axis.  More  complex  situations  have  been  handled  by 
the  simulation  in  a  similar  manner. 


Apoer.dix:  Relitionshio  Between  Mirror  Deformation  and  the  Corresponding 
'.7ave  Aberration 


As  shown  in  Figure  A-l,  consider 
a  ray  incident  at  a  point  A  on  an 
intermediate  mirror  of  an  optical 
system  making  an  angle  with  the 
normal  to  the  mirror  surface.  If  u, 
u* ,  and  n  are  unit  vectors  along  the 
incident  ray,  reflected  ray,  and  the 
surface  normal,  respectively,  we 
have  according  to  Snell’s  law 
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A'h  -I  , 
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u  +  2n  cos  v 


Suppose  that  when  the  system  is  *  *<r> 

deformed  by  a  disturbance,  a  portion 

MM  of  the' mirror  under  consideration  ric.  Ge0metrv  of  mirror  deformation 

becomes  M  M  and  that  the  point  A  .  -  d  „ava  aberra'tio.,  (AA  T  AC) . 

moves  to  a  point  A  .  The  incident 
ray,  the  point  of  incidence,  and  the 

reflected  ray,  all  deviate  slightly  from  their  nominal  positions.  Assuming  that  the  mi 
is  in  free  space,  the  wave  aberratiu..  introduced  by  the  deformation  is  equal  to  the  dif 
ance  in  path' lengths  of  the  "deformed"  and  "undeformed"  rays.  For  small  deformations, 
deformed' ray  (shown  dashed  in  Figure  A-l)  lies  close  to  the  undeformed  ray  so  that  thei 


up  to  the  first  order. 
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path  lengths  at  each  deformed  mirror.  The  deformations  are  considered  small  if  the  dis¬ 
tances  between  the  points  of  incidence  of  the  corresponding  deformed  and  undeformed  rays  on 
any  mirror  is  a  srnal  1  fraction  of  the  spacing  of  points  f'f  incidence  of  two  adjacent  rays 
of  the  undoformed  system. 

From  Fermat's  principle,  the  path  length  of  a  ray  passing  through  the  point  A"  is  sta¬ 
tionary.  Therefore,  in  calculating  the  wave  aberration  introduced  by  the  deformation,  up 
to  the  first  order,  the  point  A"  on  the  deformed  surface  can  be  replaced  by  the  point  A*  on 
it.  Hence,  the  wave  aberration,  which  is  negative  in  Figure  A-l  because  the  path  length  has 
been  reduced,  can  be  written 

3W  =  -(AB  +  AC) 


where  if  ?  and  r‘  are  the  position  vectors  of  the  points  A  and  A’,  respectively 

Sr  =  ?*  -  ?  (A-3) 

Substituting  Eq.  (A-l)  into  Eq.  (A-2),  we  obtain 

6W  =  -2n  *  Sr  cos  S  (A-4) 

In  the  simulation,  the  quantity  n  ■  6r  is  calculated  by  the  structural  program  at  struc¬ 
tural  points  Pi,  P5,  ate.,  and  it  is  calculated  at  an  optical  point  such  as  A  by  inter¬ 
polating  its  values  at  four  neighboring  structural  points.  The  angle  of  incidence  S  is 
obtained  frem  the  ray-trace  program.  The  total  wave  aberration  contributed  to  a  certain 
ray  by  all  the  mirrors  of  a  deformed  system  is  obtained  by  summing  the  aberrations  contri¬ 
buted  by  each  mirror  to  that  ray. 
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